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This  work  introduces  a new  computer  simulation  method,  known  as  the 
contact-distribution  method,  to  test  mean-field,  scaling  and  density-functional  theo- 
ries of  polymer  layers.  The  method  is  illustrated  using  a number  of  model  problems 
and  is  then  applied  to  the  case  of  polymer-bridging  interactions  between  two  bare 
surfaces  and  interactions  between  undersaturated  layers.  In  the  case  of  bridging 
interactions,  it  is  shown  that  intrachain  and  bridge/loop  and  loop/loop  interactions 
can  render  scaling  predictions  based  on  bridging  forces  alone  inadequate  even  for  ex- 
tremely low  coverages,  thereby  making  simplistic  interpretations  of  experimental  data 
inaccurate.  The  results  for  interactions  between  two  physisorbed  layers  reveal  that 
the  mean-field  theories  could  fail  in  predicting  even  the  sign  of  the  forces.  The  re- 
sults further  demonstrate  that  the  forces  of  interaction  provide  a more  discriminating 
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measure  of  the  acceptability  of  the  theoretical  approximations  than  mere  compari- 
sons of  structural  properties  (as  has  been  the  practice  so  far).  It  may  be  possible 
to  improve  the  mean-field  approach  by  using  more  accurate  chain  statistics,  and  the 
simulation  technique  described  here  provides  a systematic  method  for  identifying  and 
isolating  some  of  the  critical  assumptions  in  the  theories  so  that  the  approximations 
may  be  scrutinized  systematically. 
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CHAPTER  1 
INTRODUCTION 

1.1  Polymer-Interfaces:  What  and  Why? 

Polymers  physisorbed  at  interfaces  axe  important  both  from  a fundamental 
point  of  view  and  for  their  numerous  applications  of  technological  relevance  {e.g., 
adhesion,  lubrication,  prevention  of  protein  adsorption,  self-assembled  drug-delivery 
vehicles,  and  stabilization  or  flocculation  of  colloidal  dispersions)  [1-5].  A detailed 
knowledge  of  the  layer  structure  and  a good  understanding  of  the  static  and  dynamic 
aspects  of  the  forces  induced  by  the  layers  are  essential  for  tailoring  polymer  inter- 
faces for  specific  applications.  Experimental  efforts  to  {i)  characterize  the  surfaces 
{e.g.,  characterization  of  active  sites);  (ii)  probe  the  structural  details  of  the  layer; 
and  {Hi)  measure  the  polymer-induced  forces  constitute  an  important  step  toward 
such  ‘controlled’  design  of  polymer  interfaces.  However,  the  information  one  usually 
obtains  from  such  experiments  is  incomplete  or  invariably  contains  uncontrollable 
(or  unknown)  factors  {e.g.,  surface  heterogeneity,  impurities).  For  instance,  the  only 
structural  property  that  can  be  obtained  easily  from  experiments  is  the  average  thick- 
ness of  the  layer.  The  net  force  of  interaction  between  layers  also  can  be  obtained 
using  a number  of  techniques.  However,  the  forces  measured  result  from  a plethora 
of  interactions  {e.g.,  steric,  electrostatic,  hydrodynamic,  etc.)  many  of  which  can- 
not be  isolated  easily  for  a proper  interpretation  of  their  collective  effect.  Therefore, 
although  the  experimental  efforts  mentioned  above  are  essential  for  understanding 
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structure/property  relationships,  by  themselves  they  are  insufficient.  As  a result, 
theoretical  guidelines  that  relate  the  conformations  of  the  polymers  with  the  result- 
ing forces  are  needed.  To  this  end,  two  main  groups  of  theories,  commonly  referred 
to  as  ‘scaling’  [6]  and  ‘mean-field’  [7]  theories,  have  been  developed  to  provide  such 
guidelines.  However,  to  this  day,  there  has  been  no  attempt  to  establish  what  fun- 
damental information  of  substance  can  be  obtained  reliably  from  either  the  scaling 
or  the  mean-field  approaches.  For  instance,  even  for  the  simple  case  of  adsorption 
of  homopolymers  on  flat  homogeneous  surfaces,  it  has  become  clear  recently  that  a 
number  of  subtle  issues,  even  at  the  structural  level,  remain  as  yet  unknown  or  unex- 
plored. In  the  absence  of  rigorous,  direct  tests  of  the  available  theories,  investigators 
have  resorted  to  a comparison  of  the  theoretical  predictions  directly  with  experimen- 
tal results  of  ‘model’  systems  {e.g.,  poly(styrene)  adsorbed  on  a mica  surface)  as  an 
attempt  to  assess  the  acceptability  of  the  theoretical  predictions.  These  experimental 
efforts  have  focused,  primarily,  on  force  measurements  using  the  surface  force  appara- 
tus (SFA)  [8,  9].  These  experiments,  however,  cannot  discern  the  discrepancies  among 
the  different  theoretical  approaches,  and  are  even  less  capable  of  identifying  the  type 
of  refinements  that  could  improve  the  accuracy  of  the  theories^  . The  lack  of  rigor- 
ous tests  for  the  ‘scaling’  and  ‘mean-field’  theories  therefore  constitutes  a significant 

^ In  addition  to  non-steric  contributions  to  the  net  force,  the  existence  of  unknown 
or  uncontrollable  parameters  is  unavoidable  in  experiments.  For  instance,  as  pointed 
out  recently  [10],  systematic  errors  due  to  drift  in  the  cross-cylinder  configuration  of 
the  SFA  add  considerable  uncertainty  to  the  experimental  measurements. 
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gap  in  our  current  know-how,  and  places  a limit  on  the  improvement  of  theoretical 
guidelines  in  practice. 


1.2  Scope  of  the  Present  Work 

The  above  observations,  and  related  issues,  and  a comprehensive  review  of 
the  theories  of  polymer  adsorption  can  be  used  to  develop  a general  assessment  of 
the  state  of  our  understanding  of  polymer  layers.  Here,  we  shall  restrict  ourselves  to 
the  following  points. 

• In  general,  the  mean-field  theories  simplify  the  problem  of  a multi-chain  sys- 
tem and  focus  on  the  conformations  of  a single  chain  embedded  in  a mean-field 
potential.  At  the  outset,  one  should  take  note  that  the  above  assumption  fails 
in  systems  where  there  are  large  spatial  fiuctuations  in  the  segment  concentra- 
tion {e.g.,  semi-dilute  solutions  in  good  solvent  conditions  or  starved  layers). 
Nevertheless,  comparisons  between  the  mean-field  predictions  and  exact  calcu- 
lations of  the  structure  and  interactions  are  important  to  identify  the  precise 
boundaries  of  applicability  of  mean-field  approximations  and,  furthermore,  to 
understand  the  extent  of  deviations  one  obtains  from  such  a theoretical  ap- 
proach. Moreover,  numerical  mean-field  calculations  [7]  provide  the  only  easily 
accessible  guidelines  for  ‘short’  polymers.  Therefore,  efforts  to  examine  (and 
improve)  these  predictions  have  a large  practical  significance. 

• The  scaling  theories,  based  on  results  obtained  from  highly  fluctuating  mag- 
netic systems,  have  been  used  successfully  to  describe  polymers  in  semidilute 
solutions  [6].  These  results,  however,  are  based  on  the  use  of  simple  scaling 
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laws  and  can  be  applied  only  to  a limited  number  of  systems  {e.g.,  structure 
and  interactions  of  saturated  polymer  layers  of  infinite  molecular  weight).  A 
combination  of  scaling  arguments  and  self-consistent  mean-field  theories  has 
been  used  to  explore  the  behavior  of  undersaturated  (starved)  layers  [11,  12]. 
However,  how  well  such  (somewhat  arbitrarily  formulated)  approaches  perform 
cannot  be  determined  without  detailed  comparisons  with  exact  calculations  of 
the  structure  and  interactions. 

It  is  well  known  in  statistical  mechanics  that  computer  simulations  (computer 
‘experiments’)  provide  the  means  to  test  the  theories  in  situations  similar  to  the 
above.  In  the  case  of  polymer  layers,  a number  of  computer  ‘experiments’  have  been 
attempted  over  the  years  for  generating  the  ‘experimental’  data  needed  to  assess 
theoretical  predictions  [13,  14].  However,  most  of  these  simulations  have  focused  only 
on  the  structural  properties  of  the  polymer  layers.  Interestingly,  no  attempts  have 
been  made  to  test  the  theoretical  guidelines  using  unambiguous  calculations  of  the 
interaction  force  {i.e.,  the  Helmholtz  potential  of  interaction).  The  overall  objective 
of  this  work  is  to  advance  a new  computational  method  that  allows  one  to  determine 
forces  in  polymer/colloid  systems  from  Monte  Carlo  simulations.  This  method  can 
be  used,  as  we  have  illustrated  through  a number  of  examples  (see  Section  1.3),  to 
isolate  and  test  the  most  important  of  the  assumptions  in  mean-field,  scaling  and 
density-functional  theories.  As  shown  in  this  work,  such  a test  is  more  discriminating 
than  mere  comparisons  of  structural  properties  and  provides  the  essential  first  step 
in  developing  useful,  analytical  or  numerical  theories  for  routine  guidance  in  practice. 
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1.3  Contributions  of  the  Present  Work 
The  present  work  provides  a rigorous  test  of  the  foundations  of  the  theories 
of  polymer  adsorption.  In  particular,  this  dissertation  focuses  on  the  calculation  of 
polymer-induced  forces  using  specially  developed  computer  simulations.  Such  calcu- 
lations allow  one  to  obtain  a direct  link  between  the  structure  of  the  system  (which 
is  easily  obtained  from  simulations)  and  the  induced  forces.  In  order  to  obtain  the 
interaction  forces,  we  have  developed  a new,  efficient  and  versatile  method  for  ob- 
taining free  energies  (and,  hence,  entropic  and  enthalpic  contributions)  from  lattice 
Monte  Carlo  simulations  [15].  The  results  obtained  from  such  an  approach  provide 
details  that  are  not  accessible  through  either  theory  or  experiments,  and  hence,  ad- 
dress a significant  gap  in  our  current  knowledge.  The  key  contributions  of  the  research 
presented  here  are  as  follows: 

1.  A new  method  for  determination  of  forces  and  Helmholtz  potentials 
from  lattice  Monte  Carlo  simulations  [15].  The  calculation  of  forces  from 
simulations  of  continuous  (off  lattice)  models  makes  use  of  a well-established 
method  based  on  the  calculation  of  the  pressure  from  the  virial  theorem  [16]. 
However,  the  use  of  continuous  models  for  studying  physisorbed  layers  has  been 
restricted  to  studies  of  short  polymer  chains  (containing  less  than  32  segments) 
[17].  The  study  of  long  polymers  near  adsorbing  walls  is  better  performed 
using  the  more  efficient  lattice  simulations  [14].  The  virial  theorem,  however, 
cannot  be  used  in  lattice  models,  and  the  forces  must  be  calculated  from  a 
derivative  of  the  Helmholtz  potential.  In  this  work  we  introduce  a new,  efficient 
technique  to  calculate  forces  and  differences  in  Helmholtz  potential  from  lattice 
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Monte  Carlo  simulations  of  polymer  systems.  This  technique,  which  we  refer 
to  as  the  contact- distribution  method  (CDM),  allows  one  to  obtain  energetic 
and  entropic  contributions  of  physisorbed  or  chemisorbed  polymers  interacting 
with  arbitrary  shaped  objects.  The  calculation  of  these  forces  complements  the 
structural  information  that  is  usually  obtained  in  simulations,  thus  providing  a 
more  sensitive  test  of  theories  and  a better  insight  on  the  behavior  of  polymers 
at  interfaces  or  polymers  in  confined  spaces. 

2.  Compression  of  polymer  layers  with  finite-sized  objects  [15,  18].  The 
interactions  between  polymer  layers  and  finite-sized  objects  is  important  in  a 
large  number  of  problems,  such  as  the  dispersion  of  nanoparticles,  the  preven- 
tion of  protein  adsorption,  or  the  measurement  of  forces  using  the  tip  of  an 
atomic  force  microscope  (AFM).  It  is  common  practice  in  experimental  mea- 
surements using  an  AFM  to  relate  the  interaction  between  the  polymer  layer 
and  the  AFM  tip  with  the  free  energy  of  interaction  between  two  flat  layers. 
(This  assumption  is  known  as  the  Derjaguin  approximation  [2].)  There  are 
situations,  however,  where  the  Derjaguin  approximation  cannot  be  made.  For 
instance,  the  chains  interacting  with  a finite-sized  object  may  undergo  con- 
formational transitions  and  escape  or  squeeze-out  from  under  the  compressing 
object.  Such  transitions  give  rise  to  unusual  force  profiles.  We  have  used  CDM 
to  present  the  first  computer  simulations  that  directly  link  such  a transition 
with  a force  profile.  This  study,  which  focuses  on  the  model  case  of  a single 
end-grafted  chain  compressed  by  a flat-bottom  cylinder,  shows  the  type  of  de- 
viations from  asymptotic  predictions  (valid  for  infinite  chains)  that  one  may 
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expect  in  the  case  of  ‘short’  polymers^  . We  have  also  illustrated  that  similar 
conformational  transitions  also  occur  when  compressing  polymer  brushes  with 
finite-sized  objects.  These  transitions  give  rise  to  deviations  from  the  Derjaguin 
approximation,  thereby  indicating  that  the  interaction  force  between  a polymer 
layer  and  a finite-sized  object  cannot  be  reliably  inferred  from  the  interaction 
between  two  macroscopic  bodies  (or  vice  versa).  This  segment  of  the  present 
work  is  meant  to  serve  two  complementary  purposes:  (i)  Certain  aspects  of  the 
theoretical  (scaling)  predictions  in  the  case  of  the  escape  transition  require  no 
verification.  Therefore,  they  serve  as  a measure  of  the  accuracy  of  the  simula- 
tion method,  (n)  At  the  same  time,  the  simulations  provide  new  insights  into 
aspects  of  the  problems  that  are  beyond  the  grasp  of  theories  {e.g.,  interactions 
with  finite  objects  and  polymer  layers,  finite  objects  of  arbitrary  shape,  etc.). 
We  now  discuss  the  interaction  between  two  adsorptive  surfaces.  As  a prelude  to 
starved  (undersaturated)  surfaces,  we  first  consider  bridging  interactions  caused 
by  a single  chain  or  a few  chains  confined  between  two  walls. 

3 Bridging  attraction  induced  by  polymers  [21].  We  have  presented,  for 
the  first  time,  rigorous  calculations  of  the  attractive  force  induced  by  a single 
polymer  chain  confined  by  two  adsorbing  walls.  The  theoretical  treatment  of 

^ Our  work  on  the  compression  of  a single  end-grafted  chain  has  been  complemented 
by  recent  studies.  For  instance,  Milchev  et  al.  [19]  performed  off-lattice  simulations 
for  a long  chain  (containing  more  than  1000  segments)  in  a good  solvent,  and  Ennis 
et  al.  [20]  extended  the  theoretical  calculations  to  account  for  the  effect  of  thermal 
flcutuations  on  the  escape  transition. 
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this  case  [22]  views  the  problem  as  one  of  two  walls  connected  by  elastic  ‘teth- 
ers,’ individually  governed  by  the  Pincus  law  of  elasticity.  (Notice  that  the 
number  and  size  of  the  ‘tethers’  depends  on  both  the  adsorption  energy  and  the 
separation  between  the  surfaces.)  Such  an  assumption  gives  rise  to  a force  per 
bridge  that  is  independent  of  the  separation  between  the  walls,  and  is  solely 
determined  by  the  adsorption  energy.  Our  results  from  computer  simulations 
show,  however,  that  steric  effects  (due  to  bridge-bridge  interactions  and  to  the 
exclusion  of  conformations  by  the  walls)  may  play  an  important  role  and  dimin- 
ish the  attraction  between  the  surfaces.  Interestingly,  the  ratio  of  the  net  force 
and  the  total  number  of  bridges  is  still  roughly  independent  of  the  separation 
between  the  surfaces,  but  the  resulting  force  clearly  cannot  be  related  unam- 
biguously to  the  adsorption  energy.  Here,  the  simulations  serve  to  complement 
the  theories  and  are  used  to  identify  the  limitations  of  a naive  interpretation  of 
experimental  results. 

4 Interaction  between  two  physisorbed  layers:  Competition  between 
bridging  attraction  and  steric  repulsion  [23] . Mean- field  and  Scaling  the- 
ories differ  in  their  predictions  of  the  forces  between  undersaturated  (starved) 
layers.  While  mean-field  theories  indicate  that  attraction  between  the  surfaces 
occurs  for  small  undersaturations  (i.e.,  for  a coverage  2%  lower  than  the  satu- 
ration value),  scaling  theories  indicate  that,  for  infinitely  long  chains,  the  onset 
of  attraction  occurs  for  a rather  large  undersaturation  (i.e.,  for  a coverage  60% 
lower  than  the  saturation  value).  Resorting  to  experimental  data  to  examine 
such  discrepancies  is  not  adequate.  As  mentioned  earlier,  the  experimental 


9 


measurements  of  forces  induced  by  polymer  layers  suffer  from  a large  number 
of  uncontrolled  contributions.  This  problem  is  actually  aggravated  when  the 
forces  between  the  two  surfaces  are  attractive.  For  instance,  when  the  attrac- 
tive force  between  the  two  cross  cylinders  of  the  surface  force  apparatus  (SFA) 
is  larger  than  the  spring  constant  of  the  equipment,  a mechanical  instability  oc- 
curs and  the  two  cylinders  jump  to  contact.  The  calculation  of  polymer-induced 
forces  from  computer  simulations  remains  the  only  alternative  to  elucidate  the 
controversy  alluded  to  above  and  to  examine  the  competition  between  bridging 
attraction  and  steric  repulsion  that  occurs  in  starved  layers.  We  present  here, 
for  the  first  time,  unambiguous  calculations  of  the  force  of  interaction  between 
two  starved  layers  under  good-solvent  conditions.  Although  we  are  not  able 
to  examine  the  asymptotic  prediction  of  the  scaling  theory  (because  of  the  fi- 
nite length  of  the  chain  used  in  the  simulations),  our  results  show  important 
discrepancies  in  numerical  mean-field  calculations  and  pave  the  way  for  a rig- 
orous examination  of  the  theories  of  polymer  layers.  Although  the  numerical 
lattice  theory  for  physisorbed  polymers  was  developed  more  than  twenty  years 
ago,  the  limitations  and,  more  important,  the  refinements  needed  to  improve 
the  accuracy  of  this  numerical  theory  have  remained  unknown.  The  few  ef- 
forts to  examine  these  issues  have  relied  purely  on  the  structural  information  of 
the  polymer  layer.  The  additional  study  of  interaction  forces  provides  a more 
sensitive  method  to  assess  the  validity  of  this  theory.  Therefore,  it  provides  a 
more  appropriate  way  to  identify  the  refinements  necessary  to  develop  reliable 
predictions  for  routine  guidance  in  practice. 
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1.4  Outline  of  the  Dissertation 
The  work  we  present  here  has  been  organized  as  follows: 

• In  Chapter  2 we  present  a comprehensive  review  of  the  current  status  of  the- 
oretical understanding  of  physisorbed  layers.  As  mentioned  earlier,  two  main 
groups  of  theories  have  been  used  to  provide  guidelines  on  the  structure  of  the 
layer  and  the  induced  forces.  We  first  focus  on  the  mean-field  theories  and 
describe:  (i)  the  ground  state  dominance  approximation;  (ii)  the  two-order- 
parameter  theory;  and  (Hi)  the  lattice  numerical  theory.  We  follow  these  with 
a description  of  the  scaling  approach  and  with  an  overall  assessment  of  the  as- 
sumptions and  shortcomings  of  these  theories.  The  goal  of  this  chapter  is  to 
provide  a “bird’s-eye-view”  of  the  current  status  of  understanding  of  polymer 
layers  so  that  the  contributions  of  this  work  can  be  placed  in  context. 

• For  the  same  reasons  mentioned  above,  we  provide,  in  Chapter  3,  a review  of 
the  various  computational  methods  that  have  been  used  in  the  past  for  obtain- 
ing forces  and  Helmholtz  potentials  from  computer  simulations.  We  first  focus 
on  methods  that  have  been  used  to  examine  the  force  induced  when  stretching 
a single  chain.  We  then  provide  a description  of  the  calculation  of  the  pressure 
from  the  virial  theorem.  As  mentioned  earlier,  the  calculation  of  the  pressure 
from  the  virial  is  restricted  to  continuous  (off-lattice)  models.  In  lattice  sim- 
ulations one  has  to  calculate  the  forces  as  a derivative  (a  finite  difference,  in 
this  case)  of  the  Helmholtz  potential.  In  the  last  section  of  this  chapter  we 
focus  on  two  methods  that  can  be  used  to  obtain  differences  in  the  Helmholtz 


potential  from  Monte  Carlo  simulations,  namely,  a thermodynamic  integration 
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approach  [24]  and  the  so-called  acceptance  ratio  technique  [25].  These  two  meth- 
ods inspired  the  development  of  a new,  more  efficient  technique  for  calculating 
free-energy  differences  from  lattice  Monte-Carlo  simulations,  which  forms  the 
centerpiece  of  the  present  work. 

• Chapter  4 presents  the  technical  details  of  the  new  method,  i.e.,  the  Contact- 
Distribution  Method  (CDM);  see,  also.  Item  1 in  Section  1.3.  We  first  explain 
the  implementation  of  this  technique  and  illustrate  its  efficiency  by  comparing 
CDM  with  the  thermodynamic  integration  approach.  Then  we  present  results 
on  the  interaction  between  end-grafted  chains  and  finite-sized  objects.  The 
problem  of  compressing  a polymer  layer  with  a finite-sized  object  is  relevant  to 
the  measurement  of  forces  using  ‘small’  probes  such  as  the  tip  of  an  atomic  force 
microscope  (Item  2,  Section  1.3).  As  mentioned  earlier,  it  is  commonly  assumed 
in  practice  that  such  measurements  can  provide  information  on  the  free-energy 
of  interaction  between  two  macroscopic  bodies.  We  show,  however,  that  the 
forces  in  this  case  do  not  depend  solely  on  the  characteristics  of  the  layer  but 
also  on  some  of  the  characteristics  of  the  probe,  such  as  its  size  and  shape.  These 
results  illustrate  how  the  contact  distribution  method  allows  one  to  examine  the 
relevant  parameters  that  affect  the  forces  and,  more  importantly,  to  provide  a 
direct  link  between  the  observed  force  and  the  conformations  of  the  polymer 
chains.  Finally,  we  conclude  Chapter  4 with  an  overview  of  how  CDM  can 
be  implemented  in  problems  that  contain  additional  complexity,  such  as  the 
interaction  between  two  bodies  bearing  physisorbed  polymers.  This  problem 
becomes  the  focus  of  Chapters  5 and  6. 


12 


• As  a prelude  to  the  study  of  the  interactions  between  two  physisorbed  layers, 
we  present,  in  Chapter  5,  the  attractive  force  induced  by  a single  polymer  chain 
confined  between  two  adsorbing  surfaces  (Item  3,  Section  1.3).  The  number  and 
size  of  the  bridges  between  the  two  surfaces  depends  on  both  the  strength  of  the 
adsorption  and  the  separation  between  the  surfaces.  The  contact  distribution 
method  allows  one  to  examine  these  structural  features  and,  simultaneously, 
obtain  information  on  the  force.  We  use  this  information  to  identify  the  con- 
ditions under  which  steric  (repulsive)  contributions  become  important.  This 
work  illustrates  how  the  approach  we  have  developed  can  be  used  to  comple- 
ment theories  and  to  provide  useful  information  for  the  proper  interpretation 
of  experimental  measurements. 

• In  Chapter  6 we  present,  for  the  first  time,  unambiguous  results  on  the  inter- 
action between  two  undersaturated  (starved)  layers  (see  Item  4,  Section  1.3). 
As  already  noted,  the  calculation  of  the  interaction  between  two  polymer  layers 
is  needed  to  properly  assess  the  limitations  of  the  theoretical  predictions.  We 
provide  a comparison  between  our  simulations  and  lattice  numerical  mean-field 
calculations,  and  use  this  comparison  to  set  the  stage  for  identifying  the  most 
useful  refinements  needed  for  improving  the  accuracy  of  the  theory. 

• Some  supplementary  information  is  presented  in  two  appendices.  Appendix 
A describes  the  essential  features  of  lattice  models  for  polymer  systems  and 
Appendix  B presents  the  details  of  the  numerical  mean-field  calculations  for  a 
first-order  (standard  case)  and  a second-order  Markov  chain  {i.e.,  a chain  that 
avoids  backfolding). 
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The  research  presented  here  fills  an  important  gap  in  the  field  as  well  as  in 
the  research  at  the  University  of  Florida.  We  expect  this  study  to  have  a significant 
impact  on  the  study  of  the  theoretical  foundation  of  polymer  layers  and  interactions, 
the  development  and  refinement  of  numerical  and  analytical  theories  and  on  the 
design  of  polymer  layers  with  predetermined  response  to  perturbations. 


CHAPTER  2 

AN  OVERVIEW  OF  THE  THEORIES  OF  POLYMER  ADSORPTION 

As  noted  in  Chapter  1,  physisorbed  polymers  at  interfaces  have  received 
considerable  attention  from  the  scientific  and  industrial  community  both  from  a fun- 
damental point  of  view  and  for  their  numerous  applications  of  technological  relevance 
{e.g.,  adhesion,  lubrication,  stabilization  or  controlled  flocculation  of  colloidal  disper- 
sions, etc.).  Understanding  the  static  and  dynamic  properties  of  the  polymer  layer 
offers  the  possibility  to  control  the  structure  of  the  layer  and  to  tailor  the  properties 
of  a surface  for  a given  application. 

A model  case  which  has  been  widely  investigated  is  the  adsorption  of  linear 
flexible  polymers  onto  a uniform  flat  surface.  In  this  case,  the  structure  of  the  polymer 
layer  is  described  in  terms  of  loops,  tails  and  trains  (Figure  2.1),  and  a complete 
characterization  of  the  interface  requires  the  total  amount  of  adsorbed  polymers  and 
the  density  profile  and  the  distribution  of  the  tails  and  loops.  Understanding  the 
physics  underlying  the  behavior  of  this  simple  case  is  the  essential  first  step  in  the 
analysis  of  more  complex  cases  such  as  the  adsorption  of  copolymers,  heteropolymers 
or  polyelectrolytes  on  non-flat  and  non-uniform  surfaces. 

Two  main  groups  of  theories,  commonly  labelled  as  ‘scaflng’  and  ‘mean-field,’ 
have  been  used  to  provide  guidelines  on  the  layer  structure  and  induced  forces  and 
to  analyze  experimental  measmrements.  The  basic  concepts  behind  these  two  classes 
of  theories  are  described  below. 
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Figure  2.1:  Sketch  of  a single  polymer  chain  adsorbed  on  a surface.  A loop,  a tail 
and  a train  are  indicated  in  the  figure. 

2.1  Mean-Field  Theories 

In  general  terms,  mean-field  theories  of  interfaces  extend  the  standard  Flory- 
Huggins  picture  of  polymer  solutions  to  allow  for  the  calculation  of  free  energies  in 
systems  with  concentration  gradients.  Here  we  present  three  different  approaches  in 
which  a mean-field  assumption  is  used. 

2.1.1  Free  Energy  Functional  for  Polymers  at  Interfaces:  Cahn-de  Gennes  Approach 
and  Ground  State  Dominance  Approximation 

In  deriving  a free-energy  functional  to  describe  polymer  layers,  de  Gennes 
[26,  27]  follows  the  approach  Cahn  [28]  used  to  obtain  the  interfacial  tension  and 
wetting  properties  of  fluid/solid  interfaces.  The  interfacial  tension  7 between  a solid 
surface  and  a polymer  solution  is  decomposed  into  two  terms,  as  in  Eqn.  2.1.  The  first 
term  7d(<^s)  represents  the  interactions  at  the  surface,  while  the  other  is  a free  energy 
functional,  <j>,  (j>b),  that  accounts  for  the  nonuniform  concentration  that  develops 


16 


between  the  surface  (with  concentration  (f>s)  and  the  bulk  (with  concentration 
that  is, 


^OO 

7 = 7d(<^s)+  / I{(l>s,(f>,(l>b)dz.  (2.1) 

Jo 

The  surface  interaction  parameter,  7^,  can  be  expanded  in  a virial  series  in  terms  of 
the  concentration  at  the  surface 


7d  = 70  + 71'#’.  + 72iAS  + 


(2.2) 


in  which  the  higher-order  terms  represent  segment-segment  interactions  and  can  be 
neglected  if  the  surface  is  not  completely  filled  with  segments.  One  thus  obtains  a 
simple  relation  for  7^, 

7d  = 7o  + (2-3) 

where  70  is  the  interfacial  tension  of  the  pure  solvent  and  71  is  a (negative)  binding 
energy  between  the  segments  and  the  surface.  Equation  2.3  is  rigorously  valid  for  low 
values  of  71,  i.e.,  71  < kBTa~“^  (where  a is  the  segment  size).  This  condition,  known 
as  the  weak- coupling  limit,  is  still  compatible  with  strong  adsorption  of  polymers  since 
the  energy  or  work  required  to  desorb  a chain  with  a fraction  / of  its  segments  in 
contact  with  the  wall  is  much  larger  than  the  thermal  energy,  i.e., 

fN'Tia'^  » ksT,  (2.4) 


where  N is  the  number  of  segments  in  a chain. 

As  proposed  by  Cahn  and  Hilliard  [29],  the  second  term  in  Eqn.  2.1  {i.e.,  the 
free-energy  functional  term)  contains  a ‘local’  free  energy  contribution,  Iioc{4>)i  aiid  a 
‘nonlocal’  or  ‘gradient’  contribution,  Inioc{4>)i  which  is  assumed  to  be  a function  of  the 
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square  of  the  concentration  gradient.  In  the  particular  case  of  polymer  adsorption, 
the  local  free  energy  term  is  given  by  the  virial  expansion  of  the  Flory-Huggins  free 


energy  of  interaction, 


-fioc(0)  — 


k^T 


1 ^2  1 ..3 

-v(tf  + + ... 


(2.5) 


Note  that  we  have  intentionally  dropped  the  leading  term  in  the  Flory-Huggins  free 
energy,  i.e.,  {(I>/N)  ln0,  by  assuming  N ^ <X). 

The  non-local  term  contains  both  an  entropic  and  an  energetic  contribution. 
In  the  limit  of  ->  oo,  the  entropic  contribution  was  first  derived  by  Helfand  and 
Sapse  [30]^  ; 


■fnloc(0) 


ksT  1 f d(f> 


(2.6) 


24a  4>  \dz ^ 

For  a polymer  layer  in  equilibrium  with  a bulk  polymer  solution  of  concen- 
tration (/>6,  chemical  potential  //,  and  osmotic  pressure  H,  the  overall  expression  for 
the  interfacial  tension  is  thus  given  as 


^ r°^\kBTl  (d<j)V  keT  fl  ,2  1 

7-7o  = 7A  + _(  +6 


-f- ... ) — fi(f)  -f-  n 


dz. 

(2.7) 


The  number  of  virial  coefficients  (and  their  values)  considered  in  the  Flory- 
Huggins  virial  expansion  depends  on  the  quality  of  the  solvent.  For  instance,  for 
good  solvents  it  is  sufficient  to  keep  only  the  second  virial  coefficient  [27].^  For  theta 


^ The  details  of  the  derivation  of  Eqn.  2.6  are  also  presented  in  Ref.  [7]. 

^ The  condition  where  only  the  second  virial  coefficient  is  important  (and  the  mean- 
field  assumption  is  appropriate)  is  also  referred  to  as  marginal  solvent  condition  [7]. 
This  definition  implicitly  assumes  that  the  concentration  of  segments  is  large  enough 
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conditions  u = 0 and  hence  only  the  third  virial  coefficient  is  kept  [31].  For  poor 
solvents  both  the  second  and  third  virial  coefficients  are  used  (in  this  case  v has  a 
negative  value)  [32,  33]. 

The  objective  of  the  Cahn-de  Gennes  approach  is  to  find  the  segment  concen- 
tration profile  4>{z)  and  the  surface  concentration  (fts  that  minimize  Eqn.2.7.  At  this 
point,  an  order  parameter  tp{z)  such  that  \tl^{z)f'  = <j>{z)  is  introduced  to  simplify  the 
overall  expression  for  7,  i.e., 


1/ r 

7-7o  = 7i|V's|  / 

“ JQ 


lit)  + + + n 


dz. 

(2.8) 


Minimization  with  respect  to  ipa  leads  to  an  equation  for  the  logarithmic 
derivative  of  xp  at  the  wall,  which  is  also  used  to  define  a characteristic  length  D 
associated  with  the  so-called  proximal  region,  i.e., 


1 dtp 

Xpa  dz 


2=0— 


6071 


D‘ 


(2.9) 


Prom  the  theory  of  calculus  of  variations  [34],  minimization  with  respect  to 
xp  leads  to  the  Euler  equation 


— (—']  - — - 0 
dz  \dxp' ) dxp  ’ 


(2.10) 


SO  that  spatial  fluctuations  in  the  segment  concentration  are  small.  A marginal  solvent 
regime  is  therefore  an  ill-defined  region  between  the  semi-dilute  and  the  concentrated 
regimes. 
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where  ip'  = dip/dz  and  I = I {ip,  ip')  is  the  integrand  of  the  free-energy  functional, 
i.  e., 

I{lp,  Ip')  = Inloc{lp)  + Iloc{tp')  - + n.  (2.11) 


Since  in  our  case  I {ip,  ip')  is  not  an  explicit  function  of  z,  a first  integral  of  Euler’s 
equation  is  given  by  [34] 


(2.12) 


where  C is  a constant.  Evaluation  of  the  constant  C is  done  at  — > oo,  where  ip'  = Q 
and  ip"^  = (pb. 

A particular  case  of  practical  relevance  (see  Section  2.3)  is  one  where  a polymer 
solution  is  confined  between  two  surfaces  and  the  amount  of  polymers  between  the 
walls  remains  constant,  i.e.,  f (pdz  = const.  In  this  particular  situation  the  terms 
lj,(p  and  n are  not  considered  in  Eqns.  2.7  and  2.8,  and  the  constant  C needs  to  be 
evaluated  at  the  midpoint  between  the  surfaces,  where  ip'  = 0.  Additional  details  on 
this  specific  case,  which  is  referred  to  as  the  restricted  equilibrium  case,  are  given  in 
Section  2.3. 

Let  us  now  consider  the  situation  where  we  only  keep  the  second  virial  co- 
efficient in  the  virial  expansion.  The  chemical  potential,  which  is  given  as  //  = 
dIioc{(p) / 9(p  |ft,  becomes  ^jl  = vcpb-  After  evaluation  of  the  constant  C,  we  obtain  the 
following  expression  for  the  first  integral  of  Euler’s  equation 
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which  can  be  solved  analytically  (subject  to  the  boundary  condition  at  the  surface 
given  in  Eqn.  2.9).  Before  we  present  the  solution  to  this  equation,  it  is  instructive 
to  consider  the  derivative  of  Eqn.  2.13  with  respect  to  z, 

= (2-14) 

In  Section  2.1.2  we  will  show  that  the  above  equation  is  also  obtained  when  one  retains 
only  the  leading  term  (the  “ground  state”  term)  in  the  eigenfunction  expansion  that 
solves  the  self-consistent  Edwards  equation  [35,  36].  The  Cahn-de  Gennes  approach 
is  therefore  equivalent  to  the  so-called  ground  state  dominance  approximation. 

Returning  to  Eqn.  2.13,  we  substitute  y = and  ^ = a/y/3v^  to  obtain 

- 1.  (2-15) 


which  can  be  integrated  to  get 

4>{z)  - 0(,coth^  (~^)  ■ 


The  length  ^ introduced  above  is  the  correlation  length  calculated  in  a self-consistent 
mean-field  approximation  (also  called  the  Edwards  correlation  length).  The  constant 
zo  (which  also  has  units  of  length)  is  calculated  from  the  boundary  conditions  as 


zo  = isinh->  (^^)=fcoth-'j|. 


(2.17) 


It  should  be  pointed  out  that  the  definition  of  ^ {i.e.,  ^ = a/\/3u^)  places 
a restriction  on  the  particular  solution  shown  in  Eqn.  2.16.  The  maximum  value  ^ 
can  take  is  the  size  of  the  polymer  coils  (which  is  given  by  the  unperturbed  end-to- 
end  distance  of  the  polymer,  Rq).  This  implies  that  Eqn.  2.16  is  valid  if  the  bulk 
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solution  is  in  the  semi-dilute  regime.  (The  particular  solution  for  a polymer  layer  in 


equilibrium  with  a very  dilute  solution  will  be  discussed  shortly.) 

Three  length  scales  have  been  introduced  in  Eqns.  2.16  and  2.17,  namely, 
the  monomer  size,  a,  the  proximal  region,  D,  and  the  correlation  length,  If  these 
length  scales  differ  significantly,  i.e.,  a « D « three  spatial  regions  can  be 
differentiated.  Of  particular  importance  is  the  so-called  central  region  {D  « z « 
^).  In  this  region  Eqn.  2.16  reduces  to  the  following  asymptotic  law: 


The  power  law  given  in  the  above  equation  is  characteristic  of  all  mean-field 


power  law  will  be  contrasted  in  Section  2.2  with  the  power  law  predicted  from  scaling 
theories.  Before  proceeding  to  the  next  section,  let  us  consider  briefly  the  case  of  a 


profile  as  that  described  by  the  above  asymptotic  equation.  This  equation  will  be 
valid  for  D « z « Rq. 

In  Section  2.1.2  we  will  present  a different  approach  to  examine  the  structure 
of  the  polymer  layer.  This  approach  is  based  on  the  Edwards  equation,  which  exploits 
the  analogy  between  the  conformation  of  a polymer  chain  and  the  diffusion  equation. 

2.1.2  Beyond  Ground  State  Dominance:  Two-Order-Parameter  Theory 

The  adsorption  of  polymers  at  interfaces  can  be  described,  in  a mean-field 


(2.18) 


theories  for  polymer  layers  under  marginal  solvent  conditions  (see  footnote  2).  This 


polymer  layer  in  equilibrium  with  a very  dilute  bulk  solution  {i.e.,  <j>b  — > 0).  In  this 


particular  case.  Equation  2.13  can  be  integrated  directly  to  yield  the  same  density 


treatment,  by  the  conformations  of  a single  polymer  chain  under  the  presence  of  an 
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external  field  U{r).  In  the  particular  case  of  adsorption  on  a flat  surface,  the  external 
field  U (r)  (which  in  turn  depends  on  the  polymer  conformations  or  segment  density 
distribution)  is  only  a function  of  the  distance  from  the  surface,  z.  Such  a type  of 
problems  requires  a self-consistent  approach.  The  mathematical  description  of  the 
polymer  conformation  is  based  on  the  statistical  weight,  G{z,n),  of  a chain  that 
contains  n segments,  with  the  end-segment  being  located  at  a distance  z from  the 
surface.  As  first  shown  by  Edwards  [35],  G{z,n)  satisfies  the  following  Schrodinger- 
like  equation^  : 


dG{z,n)  a^d^G{z,n) 


— U{z)G{z,  n). 


(2.19) 


dn  6 dz'^ 

The  interactions  between  the  segments  and  the  wall  can  be  introduced  in  two  different 
ways.  One  possibility  is  to  include  an  additive  delta  function  (at  2:  = 0)  in  the  external 
field  U{z).  The  other  option  is  to  account  for  the  interaction  by  using  an  effective 
boundary  condition: 


dG{z,  n) 
dn 


G{0,n) 


(2.20) 


where  6 is  a characteristic  length  associated  with  the  strength  of  the  adsorption 
energy.  If  the  latter  method  is  used,  the  external  field  U (z)  contains  only  the  mean- 
field  excluded- volume  potential  acting  on  one  segment.  For  the  particular  case  of  a 
polymer  layer  under  marginal  conditions,  U (2)  is  given  by 


U{z)  = v(j){z), 


(2.21) 


^ A lattice  version  of  Equation  2.19,  sometimes  referred  to  as  the  propagation  equa- 
tion, is  derived  in  Appendix  B.  The  analytical  solutions  to  the  Edwards  equation 
shown  in  this  section  address  the  limit  of  infinite  chain  length. 
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where  v is  the  second  virial  or  excluded-volume  coefficient  and  (j){z)  the  segment 
concentration  at  a distance  z from  the  surface. 

The  self-consistent  approach  requires  U{z)  to  be  related  to  the  statistical 
weight  G{z,n).  This  relation  is  given  by  the  so-called  composition  rule^  , 


rN 

(j>{z)  = C I G{z,n  — m)G{z,m)dm. 

Jo 


(2.22) 


The  evaluation  of  the  constant  C depends  on  whether  we  are  considering  the  case  of  a 
polymer  layer  in  equilibrium  with  a bulk  solution  [37],  or  the  case  of  a polymer  solution 
confined  between  two  adsorbing  walls  [38].  In  what  follows,  we  will  summarize  the 
procedure  used  for  the  latter  case,  i.e.,  when  a polymer  solution  is  confined  between 
two  walls  separated  by  a distance  2h  (Figure  2.2). 

At  this  point  it  is  convenient  to  distinguish  between  two  types  of  conforma- 
tions, namely,  those  that  have  at  least  one  segment  adsorbed  on  the  surface  {i.e., 
adsorbed  conformation),  and  those  that  do  not  have  any  segment  adsorbed  on  the 
surface  {i.e.,  free  conformation).  The  statistical  weight  of  the  polymer  chains  can 
therefore  we  written  as 


G{z,n)  = G^{z,n)  + G^{z,n),  (2.23) 

where  superscript  ‘a’  stands  for  ‘adsorbed’  and  ‘/’  for  ‘free’.  This  distinction  between 
‘adsorbed’  and  ‘free’  chains  allows  one  to  obtain  the  density  profile  of  the  segments 


^ At  this  time  we  will  not  attempt  to  derive  or  describe  the  physical  meaning  of  the 
propagation  equation  or  the  composition  rule.  This  can  be  done  more  easily  when 
one  considers  the  corresponding  lattice  equations,  as  shown  in  Appendix  B. 
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Figure  2.2:  Sketch  of  a dilute  polymer  solution  confined  between  two  adsorptive  walls. 
A bridge  is  indicated  in  the  figure. 


belonging  to  loops,  (pi,  and  the  corresponding  segment  density  of  tails,  <j)t, 


<pi{z)  = C f G°'{z,ri  — m)G^{z,m)dm 
Jo 

rN 

<j>i{z)  = 2C  G‘^{z,n  — m)G^{z,m)dm 
Jo 


(2.24) 

(2.25) 


Analytical  solutions  of  Eqn.  2.19  are  obtained  as  an  eigenfunction  expansion 


(2.26) 


where  {^j(2:)}  is  a complete  orthonormal  set  of  eigenfunctions,  Cj  are  the  correspond- 
ing eigenvalues,  and  K{  are  normalization  constants  associated  with  each  eigenfunc- 
tion, i.e., 


Ki=  [ <Hidz.  (2.27) 

Jo 


The  leading  term  in  the  eigenfunction  expansion  corresponds  to  chains  in 
the  lowest  energy  state  (cq),  which  correspond  to  adsorbed  chains.  Substitution  of 
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G{z,n)  = Ko^Q{z)e  in  Eqn.  2.19  therefore  results  in  the  so-called  ground  state 
dominance  approximation: 

«o(z).  (2.28) 

This  approximation,  which  in  essence  neglects  the  effect  of  tails  in  the  adsorbed 
layer,  is  equivalent  to  the  free-energy  functional  approach  presented  in  the  Section 
2.1.1,  once  the  condition  of  self-consistency  is  introduced,  i.e.,  once  U{z)  is  written 

ast/(z)~[^oWf. 

The  goal  in  this  section  is,  however,  to  go  beyond  the  ground-state  dominance 
approximation  in  order  to  characterize  the  tail  contribution  to  the  segment  density 
profile.  Semenov  et  al.  [37,  38]  lump  the  rest  of  the  terms  in  the  eigenfunction 
expansion  and  define  a second  order  parameter  as 

cp{z)=  / G^{z,m)e^^^dm=  / {Ki^i{z)e~^^^'-^°'>^  + ...)  dm.  (2.29) 

Jo  Jo 

Integration  of  the  Edwards  propagation  equation  with  the  above  extension 
over  the  chain  length  gives  rise  to  a differential  equation  for  the  second  order  param- 
eter: 

¥ " '’"‘1  V(z)  = -1,  (2.30) 

with  boundary  conditions  (^(0)  = 0 and  dg>/dz\^_^^  = 0. 

The  total  volume  fraction  can  then  be  calculated  as  a function  of  the  two  order 
parameters,  ^'o(-2)  and  g>{z).  In  the  particular  case  of  polymers  confined  between  two 
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surfaces  the  volume  fraction  is  given  as 


(j){z)  = vTi 


(2.31) 


The  first  term  in  Eqn.  2.31  represents  the  segment  density  due  to  loops  and  the 
second  term  the  corresponding  segment  density  due  to  tails.  The  proportionality 
constant  in  Eqn.  2.31  depends  on  the  loop  contribution  to  the  surface  coverage  Tj, 
which  is  defined  as  the  total  number  of  segments  per  unit  area  that  belong  to  loops. 

It  is  important  to  point  out  that  the  absence  of  a term  proportional  to 
in  Eqn.  2.31  implies  the  assumption  that  all  the  polymer  chains  are  adsorbed  to  at 
least  one  of  the  surfaces  [39] . In  case  the  two  surfaces  are  far  from  each  other  and  the 
polymer  layer  is  in  equilibrium  with  a bulk  solution,  an  additional  contribution  due 
to  free  chains  needs  to  be  introduced  in  Eqn.  2.31.  This  additional  contribution  gives 
rise  to  a term  proportional  to  <f{zy  [37].  Details  on  the  specific  solutions  for 
and  (p{z)  for  different  conditions  are  found  in  references  [37],  [38]  and  [39].  More 
important  in  this  overview  is  to  present  the  characteristic  lengths  that  result  from 
the  analysis  of  Semenov  et  al.  The  first  length  we  have  already  introduced  is  b,  which 
is  related  to  the  strength  of  the  adsorptive  potential  and  defines  the  extension  of  the 
proximal  region.  Another  characteristic  length  is  z*,  which  represents  the  crossover  of 
the  loop  and  tail  density  profiles.  This  length  appears  when  the  mean-field  equation 
is  written  in  dimensionless  form,  and  is  defined  as  z*  = (2t;r/A^i^o)"’^^^-  Finally, 
the  other  length  of  importance  is  A,  which  is  defined  as  A = 1/^/io.  This  length  is 
related  to  the  cutoff  distance  that  defines  where  the  segment  density  starts  to  decay 
exponentially,  and  therefore  is  also  referred  as  the  layer  thickness.  Here  we  make 
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use  of  these  characteristic  lengths  to  present  the  results  obtained  for  adsorption  from 
dilute  and  semi-dilute  conditions.  The  loop  density  profile  for  these  cases  is  given  by 

0/  ~ for  b « z « z*  (2.32) 

z^ 

and 

(f)i  ~ z~^  for  z*  « z « A,  (2.33) 

whereas  the  tail  density  profile  is  given  by 

(j>t  ~ T— ^ In  f— ) for  b « z « z*  (2.34) 

{z*y  \ z / 

and 

20 

(^f  ~ — for  z*  « z « A.  (2.35) 

z^ 

As  can  be  seen  from  the  above  equations,  the  overall  density  profile  presents 
the  same  scaling  law  shown  in  Section  2.1.1,  i.e., 

(ptotal  ~ (2.36) 

The  strength  of  the  approach  we  have  summarized  briefly  here  is  that  it 
allows  for  the  determination  of  the  tail-dominant  and  loop-dominant  regions  and  for 
the  determination  of  asymptotic  laws  for  the  density  of  loops,  tails  and  free  chains. 

2.1.3  A Theory  for  ‘Short’  Polymers:  Step-Weighted  Walks  in  a Mean-Field  Potential 
As  noted  in  Sections  2.1.1  and  2.1.2,  the  analytical  predictions  of  the  ground- 
state  dominance  approximation  and  of  the  two-order  parameter  theory  are  valid  only 
for  polymers  of  infinite  molecular  weight.  (This  assumption  was  necessary  in  order  to 


28 


obtain  analytical  solutions  to  the  Edwards  propagation  equation.)  In  this  section,  we 
present  a numerical  lattice  self-consistent  field  theory  developed  by  Scheutjens  and 
Fleer  [40,  41],  which  solves  the  Edwards  equation  for  any  finite  chain  length.  (For  a 
detailed  description  of  this  theory,  see  Appendix  B). 

This  theory  considers  a lattice  system  confined  by  two  walls  and  filled  with 
polymer  segments  and  solvent  molecules.  The  lattice  has  an  area  of  L sites  (in  the 
plane  parallel  to  the  surface)  and  M layers  between  the  walls.  The  number  of  segments 
in  each  chain  is  denoted  as  r.  The  surfaces  are  located  at  z = 0 and  z = M 1, 
and  adsorption  takes  place  at  2:  = 1 and  z = M.  All  inhomogeneities  parallel  to  the 
adsorbing  wall  are  neglected.  The  underlying  lattice  structure  is  determined  by  two 
parameters^  , Aq  and  Ai,  which  are  related  to  the  fraction  of  nearest  neighbors  that 
lie  in  the  same  and  adjacent  layers,  respectively.  The  dimensionless  potential  U{z) 
that  a polymer  segment  feels  in  layer  2:  is  written  as 

U{z)  = — Xsl^lz  4 Xs2^Mz  + U*"‘(2:)  -I-  u'{z),  (2.37) 

where  Xs  is  the  so-called  Silberberg  adsorption  energy  parameter  for  the  polymer/solvent 
pair,  and  the  Kronecker  deltas,  6u  and  <5 ensure  that  the  surface  adsorption  en- 
ergy is  felt  by  only  the  segments  on  the  surface  (z.e.,  the  segments  at  layers  1 and 
M).  The  term  u\z)  is  a Lagrange  multiplier  term  which  ensures  that  the  lattice  is 
completely  filled  (this  term  is  sometimes  referred  as  a ‘space-filling  potential’  and  is 

^ The  number  of  paramenters  required  depends  on  the  lattice  structure  and  other 
details  of  the  model. 
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given  by  u'{z)  = 1 - (l){z)).  The  term  u^'^*{z)  is  an  energetic  interaction  within  the 
Bragg-Williams  or  random-mixing  approximation, 


where  x is  the  Flory-Huggins  segment-solvent  interaction  parameter,  <l>s,buik  is  the 
solvent  concentration  in  the  bulk,  and  {(l>s{z))  is  the  so-called  neighbor-average  con- 
centration around  layer  z. 

The  statistical  weight  of  a first-order  Markov  chain  is  given  in  terms  of  the  end- 
segment  distribution  function  (edf),  G{z,s),  which  represents  the  statistical  weight 
of  an  s-segment  chain  with  one  of  its  ends  in  layer  2 (and  the  other  end  being  in  any 
layer).  This  weight  factor  is  also  related  to  the  edf  of  the  preceding  segment  by  a 
recurrence  relation  (For  additional  details,  see  Appendix  B).  The  use  of  this  relation 
allows  one  to  obtain  the  statistical  weight  of  a chain  of  r segments  from  the  statistical 
weight  of  a monomer,  Gz-  The  segment  density  profile  is  then  related  to  the  weight 
factors  by  the  so-called  composition  law. 


Equation  2.39  represents  the  probability  of  joining  two  different  subchains,  both  hav- 
ing the  end-segment  in  layer  z.  (Note  that  the  prefactor  1/Gz  arises  due  to  the  double 
counting  of  the  end-segment  of  each  subchain).  The  other  term  in  the  prefactor  is 
a normalization  constant  C.  (For  instance,  for  a polymer  layer  in  equilibrium  with 
a bulk  solution  C = (j>b/r.)  A self-consistent  solution  is  then  obtained  by  assuming  a 
density  profile  and  obtaining  G{z,  s)  for  s = 1 to  r.  Once  G{z,  s)  is  known,  one  can 


M*"‘(^)  = X HM^))  - (!>!>, bulk)  > 


(2.38) 


(2.39) 
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recalculate  (j){z)  by  solving  Eqn.  2.39.  The  procedure  is  repeated  until  the  solution 
becomes  self-consistent. 

Comparisons  with  the  analytical  theories  (i.e.,  ground  state  dominance  ap- 
proximation or  two  order  parameter  theory)  have  shown  that  the  (mean-field)  asymp- 
totic laws  are  obtained  only  when  the  chains  contain  several  thousand  statistical  seg- 
ments [42,  43],  which  correspond  to  a molecular  weight  that  is  well  beyond  the  molar 
mass  of  the  polymer  additives  that  are  commonly  used  in  practice. 

The  numerical  lattice  mean-field  theory  allows  one  to  obtain  density  pro- 
files for  different  bulk  concentrations,  adsorption  energies  and  different  chain  lengths. 
Furthermore,  the  theory  can  be  modified  to  examine  a large  number  of  cases  {e.g., 
block  copolymers,  branched  polymers,  semiflexible  polymers  [7,  44]).  A careful  anal- 
ysis of  the  conditions  for  which  the  numerical  predictions  are  valid  and,  additionally, 
an  understanding  of  the  deviations  of  these  predictions  (in  those  situations  where 
mean-field  is  not  appropriate)  are,  therefore,  of  utmost  importance. 

2.2  Scaling  Theories 

The  description  of  the  spatial  fluctuations  in  polymers  under  good  solvent 
conditions  has  similarities  with  the  description  of  highly  fluctuating  magnetic  systems 
near  the  critical  temperature  [6] . Several  approaches,  which  are  based  on  the  so-called 
‘magnet  analogy’,  share  the  name  ‘scaling  theory’  in  the  scientific  literature.  Here  we 
briefly  describe  the  basis  of  these  approaches  and  point  out  some  of  the  differences 


between  them. 
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2.2.1  Self-Similar  Structure:  A Simplified  View  Based  on  Scaling  Arguments 

The  existence  of  simple  scaling  laws  that  describe  the  structure  of  bulk  poly- 
mer solutions  was  an  important  result  from  the  magnet  analogy.  These  scaling  laws, 
first  realized  by  de  Gennes  [6],  give  simple  relations  for  quantities  such  as  the  end- 
to-end  distance,  Rq,  the  characteristic  length  of  the  spatial  fluctuations  (also  called 
correlation  length  or  mesh-size),  or  the  osmotic  pressure,  11,  in  terms  of  the  length, 
N,  or  the  concentration  of  the  polymer  chains,  (j).  For  example,  for  single  unperturbed 
chains  in  good  solvents,  Ro  scales  as 

Ro  ~ (2.40) 

and,  for  semi-dilute  solutions,  the  correlation  length  and  osmotic  pressure  scale  as 

e ~ (2.41) 

n ~ (2.42) 

Equation  2.42  is  known  as  the  des  Cloizeaux  law.  The  semidilute  regime  occurs  for 
(j)h>  (f)*  fn  N/Rl,  i.e.,  when  the  polymer  coils  begin  to  overlap. 

In  order  to  examine  the  structure  of  physisorbed  layers,  de  Gennes  [45]  first 
divides  the  layer  into  different  zones: 

1.  A proximal  region  (characterized  by  a length  D),  where  the  density  of  the 
polymer  segments  depends  on  the  characteristics  of  both  the  polymer  and  the 
surface, 

2.  A central  region,  where  universal  scaling  laws  relate  the  segmental  density,  (j), 
to  the  distance  from  the  surface,  z,  and 
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Figure  2.3:  Self-similar  structure  of  a polymer  layer. 

3.  A distal  region,  where  only  a few  long  loops  and  tails  influence  the  structure  of 
the  layer. 

De  Gennes  assumes  the  central  layer  to  be  ‘self-similar’  {i.e.,  ^(^)  ~ z)  and 
uses  the  scaling  laws  developed  for  polymer  solutions  to  derive  scaling  relations  for 
the  central  region  in  the  polymer  layer  (see  Figure  2.3).  The  relations  for  the  segment 
density  obtained  in  this  manner  show  a decay  that  is  less  pronounced  than  the  one 
predicted  from  mean-field  arguments: 


(j){z)  ~ for  D«z«  ^b,  (2.43) 


where  ^6  is  the  characteristic  length  of  a bulk  solution  of  density  (f)b-  Note  that  for 
dilute  solutions,  i.e.,  when  (j)b  < </>*,  Eqn.  2.43  is  valid  for  D « z « Rq. 
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2.2.2  Free-Enersy-Functional  (FEF) /Scaling  Approach 

As  pointed  out  by  de  Gennes,  one  deficiency  of  the  pure  scaling  analysis  is 
its  lack  of  precise  formulas.  In  order  to  circumvent  this  deficiency,  and  to  obtain 
information  on  the  density  profile  for  the  entire  range  of  2:,  i.e.,  for  D < z < 


functional  used  in  the  Cahn-Hilliard  approach  (see  Section  2.1.1,  Eqn.  2.7),  which, 
for  the  case  of  segments  in  good  (marginal)  solvent,  is 


As  we  mentioned  earlier,  the  first  term  in  the  integral  represents  the  non-local  con- 
tribution and  the  second  term  the  local  contribution  based  on  the  Flory-Huggins  free 
energy  of  interaction.  In  order  to  obtain  a ‘scaling-based’  free-energy  functional,  the 
Flory-Huggins  term  is  replaced  by  the  relation  for  the  osmotic  pressure  obtained  from 
scaling  arguments  (z.e.,  the  so-called  des  Cloizeaux  law,  shown  in  Eqn.  2.42),  which 
is  written  in  terms  of  the  correlation  length  as 


where  a and  mo  are  numerical  constants.  This  equation  can  also  be  rewritten  in  terms 


de  Gennes  proposed  a ‘scaling-based’  free-energy  functional  to  replace  the  mean-field 


7-70  = 


dz.  (2.44) 


(2.45) 


The  non-local  (mean-field)  term  is  replaced  by  The  resulting  ‘scaling- 


based’  free-energy  functional  becomes 


7o 


of  the  segment  concentration  using  the  scaling  relations  between  ^ and  (j)  shown  in 
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Eqn.  2.42, 


7 - 7o  = li<f>s  + 


akeT  T 

Jo 


ml  ( d<j) 


- 


dz. 


(2.47) 


0^/“^  \dz  ^ 

The  concentration  profile  obtained  from  the  minimization  of  this  ‘scaling- 
based’  functional®  is  in  agreement  with  the  scaling  laws  [26,  27].  The  characteristic 
length  of  the  proximal  region  is  given  as 


D = a 


keT 

,l7i|a^ 


3/2 


(2.48) 


where  71  and  a are  the  binding  energy  and  monomer  size,  respectively,  as  mentioned 
in  Section  2.1.  The  asymptotic  relation  for  (j){z)  in  the  central  or  self-similar  region 


IS 


4/3 


“ (^zT'ib)  ’ D «z«  Cft,  (2.49) 

where  is  the  characteristic  or  mesh  size  of  a semi-dilute  solution.  (Note  that  for 
a very  dilute  solution  ~ i?o,  where  Rq  is  the  average  end-to-end  distance  of  the 
polymer  chains.) 

The  FEF-Scaling  approach  has  since  been  extended  by  Rossi  and  Pincus 
[11,  12]  to  the  case  of  undersaturated  polymer  layers.  It  is  our  opinion,  however,  that 
the  free-energy  functional  in  Equation  2.47  was  defined  in  an  arbitrary  manner.  The 
limitations  of  such  an  approach  need  to  be  examined  carefully. 


® This  approach  has  also  been  called  ‘renormalized  mean-field  theory’  [38]. 
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2.3  Scaling  or  Asymptotic  Laws  for  the  Interaction  Force 
In  the  previous  sections  we  have  focused  on  the  structure  of  polymer  layers. 
The  relation  between  such  structure  and  the  induced  forces  is  essential  to  tailor  poly- 
mer layers  for  specific  applications.  It  has  been  observed,  in  experimental  measure- 
ments, that  the  interaction  energy  between  surfaces  with  physisorbed  chains  depends 
strongly  on  the  ‘history’  of  the  sample  {i.e.,  incubation  time,  washing  steps,  etc.)  and 
the  characteristic  timescale  of  compression  [9,  26].  This  observation  has  led  to  the 
proposition  of  two  basic  cases  (for  theoretical  analysis),  namely,  full  and  restricted 
equilibrium.  These  cases  are  described  in  the  following  subsections. 

2.3.1  Full  Equilibrium 

The  full-equilibrium  situation  rests  on  the  assumption  that  the  polymer  chains 
between  the  two  surfaces  are  in  equilibrium  with  a bulk  solution  of  concentration 
(j)b.  Under  this  condition  it  has  been  shown  by  both  scaling  and  mean-field  theories 
(ground  state  approximation)  that  there  is  an  attraction  between  the  two  surfaces 

/ ~ -h~^.  (2.50) 

Calculations  with  the  two-order-parameter  theory  show  a small  repulsion  at  long  dis- 
tances (due  to  the  interaction  between  the  tails  of  the  adsorbed  polymers)  followed 
by  attraction  [46].  This  attraction  occurs  because  several  chains  desorb  from  the  sur- 
faces and  diffuse  to  the  bulk  as  the  surfaces  come  closer  to  one  another.  Experiments, 
however,  show  repulsion  between  the  plates,  since  the  chains  are  so  strongly  attached 
to  the  surface  that  they  cannot  desorb  or  diffuse  within  the  timescale  of  compres- 
sion. Furthermore,  it  is  believed  that  most  of  the  experimental  me£isurements  have 
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been  taken  under  conditions  where  the  total  amount  of  polymers  between  the  inter- 
acting surfaces  remains  constant.  The  above  condition  has  been  called  irreversible 
adsorption  by  some  researchers  [37],  although  restricted  equilibrium^  may  be  a better 
description  of  the  situation. 

2.3.2  Restricted  Equilibrium 

The  forces  obtained  under  restricted  equilibrium  depend  on  whether  the  sur- 
face coverage,  T,  defined  as  the  total  number  of  segments  per  unit  area  belonging  to 
adsorbed  chains,  is  larger  or  smaller  than  the  saturation  coverage,  Tq.  The  saturation 
coverage  is  defined  as  the  coverage  obtained  when  the  polymer  layer  is  in  equilibrium 
with  a very  dilute  polymer  solution. 

Forces  between  saturated  surfaces 

For  good  solvent  conditions,  the  scaling/ free-energy-functional  approach  pre- 
dicts a repulsive  force  at  all  distances.  The  scaling  of  the  force  with  the  separation 
between  surfaces  is  [26] 

/ ~ for  h»D  (2.51) 

and,  for  a very  small  separation  or  gap  between  the  surfaces, 

/ ~ for  h « D,  (2.52) 

^ The  term  ‘irreversible  adsorption’  is  not  appropriate  since  one  assumes  that  the 
chains  are  allowed  to  equilibrate  within  the  gap,  z.e.,  the  timescale  for  desorption 
and  diffusion  within  the  gap  between  the  plates  is  smaller  than  the  timescale  of 
compression. 
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where  D is  the  characteristic  dimension  of  the  proximal  region,  as  defined  in  Eqn. 
2.48. 

Attempts  to  obtain  the  force  relations  for  a saturated  layer  under  good 
(marginal)  solvent  conditions  using  the  ground  state  dominance  approximation  (Cahn- 
de  Gennes  approach)  lead  to  / ~ 0.  However,  when  the  calculations  are  performed 
using  the  two-order-parameter  mean-field  theory  [38],  a steric  repulsion  is  obtained. 
This  difference  is  partly  due  to  the  fact  that  the  ground  state  dominance  approxi- 
mation underestimates  the  saturation  coverage  and  neglects  the  effect  of  tails.  The 
asymptotic  relations  obtained  with  the  two-order-parameter  theory,  for  the  proximal 
region,  are: 

/■proximal  ~ h for  h <<  b,  (2.53) 

whereas  in  the  loop-dominated  and  tail-dominated  regimes  one  has,  respectively, 

/ ~ h~^  for  b«h«  z*,  (2.54) 

and 

/ ~ for  z*  <h<  Rend-  (2.55) 

Forces  between  undersaturated  surfaces 

Rossi  and  Pincus  [9,  12]  considered  the  effect  of  surface  coverage  on  the  in- 
teraction force  under  good  solvent  conditions.  Their  results,  which  are  an  extension 
of  the  scaling/free-energy-functional  approach,  show  that  for  moderate  undersatura- 
tions, i.e.,  0.5  < r/Fo  < 1.0,  the  forces  become  less  repulsive  (as  compared  to  the 
saturated  case).  For  large  undersaturations,  i.e.,  F/Fq  < 0.4,  an  attractive  force 
occurs,  presumably  due  to  the  presence  of  bridges. (see  Figure  2.4). 
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Scaling  / Free-Energy  Functional  Approach 
• Attraction  occurs  for  n r,  < 0.4 


Figure  2.4:  Sketch  of  the  interaction  between  two  undersaturated  layers  as  predicted 
by  the  scaling/free-energy-functional  approach  [9,  12]. 

It  is  interesting  to  contrast  the  ‘scaling’  results  with  the  analytical  predic- 
tions based  on  the  Cahn-de  Gennes  free-energy  functional  approach  (ground  state 
dominance  approximation)  and  the  two-order-parameter  theory.  The  ground  state 
dominance  approximation  predicts  attractive  forces  for  any  F/Fo  < 1/0  (and  / = 0 
for  F/Fq  = 1.0),  while  the  two-order-parameter  theory  predicts  a decrease  in  repul- 
sion for  very  small  undersaturations,  and  attraction  for  any  F/Fq  < 0.98  (see  Figure 
2.5). 

One  should  keep  in  mind  that  the  scaling/free-energy-functional  approach  was 
not  derived  from  first  principles,  and  hence,  the  validity  of  the  extensions  by  Rossi 
and  Pincus  is  unknown.  In  addition,  the  mean-field  predictions  are  supposed  to  be 
valid  for  marginal  solvent  conditions,  which  implies  that  the  segment  concentration 
is  between  the  semi-dilute  and  the  concentrated  regimes.  Whether  that  assumption 
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Figure  2.5:  Sketch  of  the  interaction  between  two  undersaturated  layers  as  predicted 
by  mean-field  theories  [39]. 

is  satisfied  in  the  case  of  undersaturated  layers  is  uncertain.  The  generation  of  unam- 
biguous force  profiles  for  undersaturated  layers  is,  therefore,  of  utmost  importance. 

Force  profiles,  for  theta  and  poor  solvent  conditions,  have  been  obtained  us- 
ing the  mean-field  Cahn-de  Gennes  approach  (for  both  saturated  and  undersaturated 
surfaces).  It  would  be  interesting  to  compare  these  predictions  with  numerical  cal- 
culations based  on  the  Scheutjens-Fleer  theory  (which  incorporates  a solvent  quality 
parameter). 

Non-equilibrium  forces  between  polymer  layers 

Besides  the  two  cases  we  have  considered  above,  namely,  reversible  adsorp- 
tion and  restricted  equilibrium,  one  can  envision  different  situations  where  a higher 
degree  of  irreversibility  may  occur.  For  example,  a case  named  ‘double  restricted 
equilibrium’  (DRE),  in  which  not  only  the  surface  coverage  but  also  the  concentra- 
tion of  monomers  at  the  surface  is  kept  constant,  has  been  considered  by  Brooks  and 
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Figure  2.6:  Sketch  of  the  interactions  between  two  saturated  surfaces  under  restricted 
equilibrium  conditions.  The  forces  predicted  by  the  double-restricted  equilibrium 
condition  [47]  are  larger  than  the  ones  predicted  by  (standard)  restricted  equilibrium. 
The  scaling  laws,  however,  remain  the  same. 

Cates  [47].  This  modification  was  implemented  within  the  ‘scaling’/FEF  approach. 
The  numerical  results  indicate  the  presence  of  larger  forces  when  compared  with  the 
restricted  equilibrium  case  (see  Figure  2.6).  The  additional  restriction,  however,  does 
not  affect  the  scaling  law  in  the  central  regime  {h  » D)  but  only  the  one  in  the 
proximal  regime  [h  « D),  i.e., 


foRE  ~ h for  h « D. 


(2.56) 


The  limiting  case  of  permanent  adsorption  (PA),  e.g.,  complete  ‘irreversibil- 
ity’, was  considered  by  Aubouy  et  al.  [48].  Their  results,  based  on  simple  scaling 
arguments,  show  that  a permanently  adsorbed  layer  exhibits  the  same  behavior  as  a 
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polymer  brush  under  compression, 

fpA  ~ h-ol*.  (2.67) 

The  cases  shown  here  are,  however,  just  meant  to  provide  an  insight  into  the 
forces  induced  under  non-equilibrium  situations.  The  predictions  of  these  forces  still 
remains  an  almost  unexplored  area. 

2.4  Experimental  Studies  on  the  Layer  Structure  and  Induced  Forces 
In  the  previous  sections  we  have  summarized  some  of  the  theoretical  pre- 
dictions for  the  density  profile  and  the  force  of  compression  of  physisorbed  polymer 
layers.  The  structure  of  such  layers  has  also  been  studied  by  different  experimental 
techniques  (in  addition  to  computer  simulations).  The  experimental  efforts  that  have 
been  made  so  far  can  be  classified  broadly  into  three  categories,  namely,  those  which 
focus  on  the  amount  of  adsorbed  polymers  {global  features),  those  which  give  infor- 
mation on  the  microstructure  of  the  layer  and  those  focusing  on  force  measurements. 
We  present  a brief  discussion  of  some  of  these  here. 

1.  Hydrodynamic  measurements,  ellipsometry  [45]  and  evanescent  wave  induced 
fiuorescence  [50]  are  some  of  the  techniques  that  have  been  used  to  measure  the 
layer  thickness  or  the  total  amount  of  adsorbed  polymers.  These  experiments 
have  confirmed  that  the  thickness  of  the  adsorbed  layer  is  of  the  order  of  the  coil 
size  and  are  used  extensively  in  studies  of  adsorption  kinetics  or  competitive 
adsorption. 

2.  Neutron  scattering  techniques,  based  on  either  small-angle  scattering  or  reflec- 
tivity experiments,  have  been  used  to  probe  the  density  profile  of  adsorbed  and 
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end-grafted  layers  [51,  52].  The  results  obtained  with  these  techniques  are  in 
qualitative  agreement  with  theoretical  predictions,  but  cannot  be  used  to  make 
definite  (quantitative)  comparisons  with  theories. 

3.  The  third  group  of  experiments  focuses  on  the  interaction  force  between  ph- 
ysisorbed  layers.  Most  of  such  force  experiments  have  made  use  of  the  sur- 
face force  apparatus  (SFA)  and  have  measured  the  forces  under  good,  theta 
and  below-theta  conditions  [8,  9].  Again,  analysis  of  the  experimental  data 
shows  that  the  results  are  in  qualitative  agreement  with  theories,  but  no  def- 
inite comparisons  with  the  theories  are  possible  due  to  the  large  number  of 
uncontrolled  factors  that  exist  in  the  experimental  realization.  For  example,  it 
has  been  observed  experimentally  that  polymer  chains  do  not  adsorb  uniformly 
on  surfaces  but  form  ‘islands’  [53,  54].  Therefore,  it  is  difficult  to  use  ‘average’ 
measurements  such  as  the  ones  obtained  from  SFA  and  compare  those  forces 
with  theoretical  predictions.  Atomic  force  microscopy  (AFM)  can  be  used  as  a 
‘local’  probe  in  this  respect,  but  one  must  pay  attention  to  the  influence  of  the 
shape  and  size  of  the  AFM  tip  and  its  affinity  (or  lack  thereof)  to  the  polymer 
chains  in  order  to  interpret  the  measurements  meaningfully.  Moreover,  current 
theories  ignore  lateral  inhomogeneities  and  are  inadequate  for  predicting  the 
force  of  compression  by  finite-sized  objects. 

Due  to  the  severe  limitations  found  in  experiments,  computer  simulations 
serve  as  an  ideal  tool  to  examine  the  theoretical  predictions  on  the  structure  and 
forces.  In  the  next  section  we  comment  briefly  on  the  type  of  results  that  have  been 
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obtained  in  previous  simulations  studies  and  the  particular  issues  that  need  to  be 
addressed. 

2.5  Outstanding  Issues:  What  is  the  Role  of  Computer  Simulations? 

The  role  of  computer  simulations  is,  actually,  twofold:  (i)  as  a test  of  scaling 
and  theoretical  predictions  and  {ii)  as  a predictive  tool  to  get  insight  in  cases  that 
cannot  be  easily  analyzed  theoretically. 

Several  simulations  have  been  reported  in  the  literature  for  the  structure  and 
dynamics  of  physisorbed  polymers  [55-61].  While  the  results  available  confirm,  at 
least  in  qualitative  terms,  the  picture  of  polymer  layers  developed  theoretically  {e.g., 
qualitative  distribution  of  loops,  trains  and  tails)  and  are  in  agreement  with  the  over- 
all segment  density  profile  expected  from  the  theories,  no  simulations  are  available 
on  detailed  quantitative  predictions  of  structural  features  and  segment  distributions. 
Moreover,  the  range  of  adsorption  energies  and  molecular  weights  for  which  the  avail- 
able theories  are  valid  have  not  been  examined  rigorously.  (In  general,  the  theoretical 
approaches  presented  in  this  chapter,  with  the  exception  of  the  Scheutjens-Fleer  the- 
ory, are  strictly  valid  for  very  high  molecular  weights  and  relatively  low  or  moderate 
adsorption  energies  - constraints  that  are  not  necessarily  satisfied  in  practical  sit- 
uations, e.g.,  in  the  case  of  polymers  used  in  colloid  stabilization.)  A first  step  in 
this  direction  has  been  taken  recently  by  de  Joannis  and  Bitsanis  [14],  who  have  ini- 
tiated a careful  examination  of  the  structure  of  physisorbed  layers  for  a wide  range 
of  adsorption  energies  and  molecular  weights  and  have  developed  an  efficient  Monte 
Carlo  algorithm  to  sample  the  conformations  of  long  polymer  chains  [62] . This  work. 
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however,  leaves  an  important  question  unanswered;  What  is  the  relation  between 
the  structure  of  the  polymer  layer  and  the  induced  forces,  and  how  does  a change 
in  the  density  profile  affect  the  forces?  It  is  therefore  evident  that  the  structural 
information  on  the  polymer  layers  should  be  complemented  with  information  on  the 
polymer-induced  forces.  The  calculation  of  the  induced  forces  is,  however,  not  an 
easy  task,  and  no  attempt  has  been  made  so  far  to  examine  the  forces  due  to  the 
compression  of  physisorbed  layers  and  the  relation  between  the  structure  of  the  layer 
and  the  resulting  forces.  In  the  next  chapter  we  will  review  some  of  the  methods  that 
can  be  used  to  obtain  forces  from  computer  simulations. 

Finally,  it  is  important  to  mention  that,  in  addition  to  being  an  ideal  tool  to 
examine  the  theoretical  predictions,  computer  simulations  are  also  very  important  for 
probing  the  structure  and  forces  in  cases  that  cannot  be  easily  examined  by  purely 
theoretical  means.  Examples  of  such  cases  are  the  compression  of  polymer  chains  or 
layers  by  finite-sized  objects,  adsorption  of  copolymers  and  heteropolymers  [63-65], 
adsorption  on  non-uniform  surfaces  [53,  54],  competitive  adsorption  [66-69],  etc. 


CHAPTER  3 

POLYMER-INDUCED  FORCES  FROM  COMPUTER  SIMULATIONS 
Computer  simulations  are  useful  both  in  testing  the  validity  of  theoretical 
models  as  well  as  in  helping  in  the  interpretation  of  the  experiments.  They  provide 
details  that  are  inaccessible  to  either  theory  or  experiments.  For  instance,  one  can 
use  simulations  to  examine  physisorbed  polymers  at  interfaces  in  order  to  obtain 
structural  features  such  as  the  fraction  of  tails,  trains,  loops  and  bridges.  One  can  go 
a step  further  and  get  details  on  the  number  and  size  distributions  of  those  entities, 
for  a wide  range  of  adsorption  energies  and  chain  lengths.  Additional  information 
on  the  induced  force  allows  one  to  make  a direct  link  between  the  forces  and  the 
conformations.  Although  this  link  is  strictly  valid  only  for  the  particular  model  used 
in  the  simulations,  the  information  gained  from  these  types  of  studies  gives  a better 
insight  into  structure-property  relationships  in  practical  situations. 

The  choice  of  a particular  model  used  in  computer  simulation  studies  depends 
on  the  specific  objectives  of  the  simulation.  For  instance,  it  depends  on  whether  one 
is  interested  in  the  static  or  dynamic  properties.  It  is  also  a compromise  between 
computer  time  and  the  level  of  realism  (or  coarse-graining)  required  in  the  model. 
Our  goal  here  is  to  summarize  some  of  the  methods  one  can  use  to  obtain  polymer- 
induced  forces  from  either  lattice  or  continuum  models.  We  first  focus,  in  Section 
3.1,  on  the  elastic  forces  that  arise  when  stretching  single  polymer  chains  (as  in  the 
case  of  a stretched  polymer  bridge  between  two  surfaces).  In  Section  3.2  we  discuss 
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the  calculation  of  the  pressure  tensor  from  off-lattice  simulations.  Calculation  of 


the  pressure  from  the  virial  equation  provides  a robust  method  to  estimate  forces 
induced  by  polymers  that  are  confined  by  either  repulsive  or  attractive  walls.  Finally, 
in  Section  3.3  we  turn  our  attention  to  more  general  methods  that  can  be  used  to 
estimate  free-energy  (i.e.,  Helmholtz  potential)  differences.  These  methods  are  very 
useful  to  estimate  forces  in  lattice  simulations. 

3.1  Stretching  Single  Polymer  Chains:  Universal  Force-Length  Relations 
The  elastic  behavior  of  single  polymer  chains  is  of  significant  importance  for 
understanding  the  behavior  of  such  macromolecules  in  solution  or  when  confined  in 
small  spaces.  In  the  case  of  ideal  chains,  or  for  chains  under  theta  conditions,  the 
force,  /,  required  to  stretch  a polymer  coil  contains  only  an  entropic  contribution. 


(3.1) 


where  T is  the  temperature,  S the  entropy  and  R (=  |R|)  the  end-to-end  distance  of 


the  polymer  chain.  For  Gaussian  chains,  Eqn.  3.1  can  be  rewritten  (in  vector  form) 


as 


(3.2) 


where  ks  is  the  Boltzmann  constant,  N the  number  of  segments  in  the  chain  and  a 


the  segment  size. 


For  polymers  under  good-solvent  conditions,  the  excluded-volume  interactions 
modify  the  force-length  behavior.  In  the  case  of  weak  forces  (or  weak  stretching),  the 
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force  is  given  as 


/ ~ 


keT 

Ne/5a2 


R, 


for 


/< 


keT 

N^/^a 


(3.3) 


where  = Rq  is  the  unperturbed  end-to-end  distance  of  the  polymer  coil.  For 

the  strong  stretching  regime,  i.e.,  ksT/Ro  < f < kgT/a,  the  polymer  coil  does  not 
deform  uniformly  (see  Figure  3.1),  but  breaks  up  into  a series  of  “tensile  blobs”  of 
size  [70].  (A  blob  is  a region  in  space  where  the  chain  retains  the  local  correlation 
of  a Flory  chain,  i.e.,  a region  in  space  where  the  chains  remain  unperturbed.)  The 
size  of  the  tensile  blobs  depends  on  the  force  intensity  as 

~ (3.4) 

and  is  also  related  to  the  number  of  segments  it  contains,  g,  by 


(3.5) 


The  chain  elongation  is  then  proportional  to  the  total  number  of  blobs  as 


R"  ” 7^"  ~ (k^) 


(3.6) 


which  can  also  be  rewritten  as 


r ^ ^bT  ^3/2 


(3.7) 


Notice  that  Rp  in  the  above  equation  is  the  projection  of  the  end-to-end  distance, 
R,  on  the  direction  of  the  force,  f,  and  Cp  is  a proportionality  constant  that  is 
model-dependent.  Equation  3.7  is  commonly  known  as  the  Pincus  law. 
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Figure  3.1:  Entropic  elasticity  of  single  polymer  chains  under  good  solvent  conditions. 

The  concept  of  ‘tensile  blobs’  has  applications  in  other  problems  in  polymer 
physics  such  as  polymer  bridging  and  polymers  under  shear.  The  use  of  computer  sim- 
ulations to  examine  these  force-length  relations  and  the  influence  of  additional  factors 
such  as  chain  flexibility,  solvent-segment  interactions  or  the  presence  of  obstacles  {e.g., 
a nearby  surface)  is  important  to  get  a better  insight  into  the  deformation  behavior 
of  polymer  chains.  In  what  follows,  we  present  a description  of  some  of  the  methods 
that  can  be  used  to  obtain  ‘universal’  force-length  relations  using  coarse-grained  or 
mesoscopic  models  (z.e.,  models  in  which  the  chemical  details  of  the  macromolecule 
are  not  important). 

3.1.1  Strain  Ensemble 

In  the  strain  ensemble  the  end  atoms  of  the  polymer  chain  are  fixed  at  a 
distance  Rend  from  each  other,  and  one  monitors  the  force  f needed  to  maintain  this 
distance.  This  force  fluctuates  as  the  rest  of  the  chain  undergoes  thermal  motion. 
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In  off-lattice  simulations,  the  contributions  to  this  force  consist  of  the  non-bonded 
interactions  between  the  end  segments  and  the  rest  of  the  segments  in  the  chain 
and  the  ‘internal’  force  arising  from  the  fluctuations  in  the  bonds  between  the  end 
segments  and  their  immediate  neighbors.  The  average  force  is  thus  obtained  as 


Ir  — 


(f)  • Rend 

Rend 


(3.8) 


where  /r  is  the  component  of  the  force  in  the  direction  of  the  end-to-end  vector. 

An  interesting  variation  of  this  method,  especially  suited  for  Monte  Carlo 
simulations,  was  recently  introduced  by  Wittkop  et  al.  [71,  72].  In  their  procedure, 
the  end-to-end  distance  in  the  z direction  was  initially  fixed  at  and  the  force 
was  measured  by  counting  the  attempted  jumps  with  + a and  n_  with 

Rz  Rz  — a,  where  a is  the  lattice  spacing.  (Note  that  once  an  end-segment  is 
selected,  the  attempted  move  is  used  to  obtain  n+  or  n_,  but  the  position  of  the 
end-segments  is  not  updated.)  Then,  the  force  in  the  z direction  was  computed  as^ 


ksT  n. 

fz  = In  — 


(3.9) 


a n+ 

Wittkop  et  al.  [71]  used  this  method  to  study  the  deformation  behavior  at 
temperatures  below  theta  conditions  using  lattice  simulations  of  the  bond  fluctuation 
model.  Their  results  confirmed  that,  below  the  theta  temperature,  a phase  transition 
from  a collapsed  globule  to  a swollen  coil  occurs  and  that  the  transition  region  is 
characterized  by  a coexistence  between  a globule  and  a stretched  strand.  A similar 


^ A derivation  of  this  equation  was  recently  presented  by  Kreitmeier  et  al.  [73]. 
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Figure  3.2:  Competition  between  entropic  and  enthalpic  forces  when  stretching  single 
polymer  chains. 

balance  between  the  entropic  and  enthalpic  contributions  to  the  free  energy  of  a 
polymer  is  also  expected  when  stretching  (and  detaching)  a chain  that  is  strongly 
adsorbed  to  a surface  (see  Figure  3.2)  or,  as  pointed  recently  by  Kreitmeier  et  al. 
[73],  when  stretching  a folded  protein. 

3.1.2  Stress  Ensemble 

In  the  stress  ensemble  one  end  of  the  chain  is  fixed  at  the  origin,  and  a fixed 
force  F is  applied  to  the  other  end.  The  distance  Rp,  the  projection  of  the  end-to-end 
distance  Re„d  along  the  direction  of  the  force,  now  fluctuates.  The  average  value  of 
Rp  is  then  obtained  as. 


Rp  = (Rend)  • F/F. 


(3.10) 
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In  Monte  Carlo  simulations  the  force  F is  introduced  in  the  Boltzmann  factor, 
w.  In  the  usual  Metropolis  sampling  scheme,  w is  written  for  the  end  segments  as 

Wendsegment  = exp  [-/S  (AC/  - ~ Rfm))]  (^-H) 

where  AC/  is  the  change  in  the  configurational  energy.  (Note  that  for  all  other  seg- 
ments w is  given  in  terms  of  the  change  in  the  configurational  energy  alone.) 

It  is  important  to  note  that  the  force-length  relationship  obtained  in  the 
stress  ensemble,  i.e.,  F vs.  Rp,  is  not  equivalent  to  the  one  obtained  in  the  strain 
ensemble,  fp  vs.  Rend-  This  non-equivalence  has  been  shown  rigorously  [74,  75,  76], 
and  can  be  easily  understood  once  the  difference  in  the  definition  of  ‘length’  in  each 
ensemble  is  stated.  In  particular,  a difference  between  Rend,  the  end-to-end  distance, 
and  Rp,  the  projection  of  Rend  over  F,  is  expected  to  occur.  Moreover,  this  difference 
is  a strong  function  of  the  quality  of  the  solvent  and  is  especially  significant  for  low 
deformations.  However,  for  moderate  to  large  deformations,  the  difference  decreases 
and  both  ensembles  give  similar  results  [71,  77]. 

Simulations  based  on  the  stress  ensemble  have  been  used  to  study  force-length 
relations  using  the  athermal  bond-fluctuation  model  [78].  The  results  obtained  in  this 
study  were  in  good  agreement  with  the  scaling  laws  that  predict  a linear  deformation 
for  small  forces  and  a nonlinear  relation,  known  as  the  Pincus  scaling  law,  for  larger 
forces.  The  proportionality  constant  obtained  for  the  Pincus  law  was  Cp  « 3.25.^ 

^ More  recently,  off-lattice  simulations  of  very  long  polymer  chains  {N  « 6000) 
[79]  also  confirmed  the  predictions  of  the  Pincus  law.  The  proportionality  constant 
obtained  in  this  study  was  Cp  « 3.2. 
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It  was  also  shown  that  for  strong  stretching  (ie.,  when  R > 2Rq)  the  results  deviate 
from  the  Pincus  law  and  become  model-dependent.  In  this  regime  the  energetic 
characteristics  of  the  model  {i.e.,  stretching  and  twisting  of  the  bonds)  dominate  over 
the  entropic  forces. 

3.1.3  Non-Universal  Force  -vs-  Length  Relationships 

The  recent  development  of  experimental  techniques  to  manipulate  single 
macromolecules,  e.g.,  using  optical  traps  [80,  81],  magnetic  tweezers  [82]  or  an  atomic 
force  microscope  [81,  83,  84],  allows  one  to  measure  intrachain  forces  more  directly 
than  has  been  possible  so  far.  These  forces,  which  fall  in  the  strong  stretching  regime, 
lead  to  non-universal  force-length  relationships  that  have  provided  information  on  the 
secondary  and  tertiary  structures  of  several  biopolymers.  For  instance,  the  double- 
stranded  DNA  helix  exhibits  a reversible,  conformational  transition  at  a force  of  0.06 
nN  to  a more  tightly  coiled,  extended  structure  [85,  86];  dextran  [87]  and  xanthan 
also  undergo  conformational  transitions  at  0.25  and  0.4  nN,  respectively;  and  titin 
exhibits  a highly  inhomogeneous  deformation  {i.e.,  a sawtooth  pattern)  [88,  89,  90]. 
Similarly,  the  deformation  of  synthetic  polymers  [91]  allows  one  to  obtain  the  elastic- 
ity of  the  macromolecule  and  to  examine  the  conformation  of  adsorbed  or  tethered 
chains  [92,  93,  94].  Analysis  of  the  experimental  data  from  such  experiments  is  very 
often  based  on  simple  theoretical  models  such  as  the  freely  jointed  chain^  and/or 

^ This  model  consists  of  n segments  of  length  6 joined  by  freely  rotating  pivots.  The 
classical  model  considers  that  the  length  of  each  segment,  which  may  be  interpreted 
as  the  Kuhn  statistical  segment  length,  remains  fixed.  There  are,  however,  some 
extensions  where  the  segments  have  a finite  elasticity. 
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the  wormlike  chain  model^  . Although  such  ideal  models  are,  in  some  cases,  useful 
to  estimate  the  length  of  the  statistical  segment  of  the  polymer  chain  or  to  estimate 
the  entropic  and  enthalpic  contributions  to  the  total  force,  they  place  a limit  on  the 
information  and  understanding  one  could  possibly  gather  from  such  micromanipula- 
tion experiments.  In  this  context,  atomistic  molecular  dynamics  simulations  might 
be  helpful  to  interpret  the  experimental  data.  Examples  of  studies  where  experiments 
and  simulations  have  been  combined  are  given  by  Grubmuller  et  al.  [95],  who  ob- 
tained the  streptavidin-biotin  bond  strength  from  computer  simulations,  and  by  Rief 
et  al.  [87],  who  used  an  AFM  tip  to  stretch  a polysaccharide  that  was  chemisorbed 
to  a surface.  The  method  used  in  these  two  cases  belong  to  neither  the  strain  en- 
semble nor  the  stress  ensemble,  i.e.,  it  represents  a variable  stress  and  variable  strain 
simulation. 

A common  implementation  of  a variable  stress-variable  strain  simulation  (to 
obtain  force-length  relationships)  resorts  to  a non-equilibrium  molecular  dynamics 
simulation.  Here  one  considers  that  a segment  or  atom  of  the  molecule  is  subject  to  a 
harmonic  potential,  V (which  basically  implies  that  the  segment  is  attached  to  either 
an  Atomic  Force  Microscope  (AFM)  cantilever  or  to  a spherical  particle  in  an  optical 
trap).  This  potential,  which  is  shifted  at  a constant  speed  in  the  z direction,  is  given 
by 

V = kc  [Zattachedit)  ~ ^o(^)]^  /2,  (3.12) 

^ The  wormlike  chain  describes  a semiflexible  polymer,  i.e.,  one  that  possesses  local 
stiffness  and  long-range  flexibility. 
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where  is  the  spring  constant  of  the  potential,  Zattached{t)  is  the  position  of  the  atom 
connected  to  the  cantilever  or  the  particle  at  any  time  t and  Zo{t)  is  the  position  of 
the  cantilever  or  the  particle,  which  are  pulled  at  a constant  velocity,  Uq-  (Note  that 
this  expression  neglects  the  thermal  fluctuations  of  the  cantilever  or  trapped  particle.) 
One  can  rewrite  zo{t)  as 

Zo{t)  = 2!o(^  = 0)  + Vot,  (3.13) 

where  Zo{t  — 0)  actually  coincides  with  the  position  of  the  attached  atom  at  t = 0. 
The  exerted  force  follows  from  Hooke’s  law, 

f ~ kc[ZaUached{^)  2^0 (t)].  (^’l^) 

Therefore,  at  any  time  t during  the  simulation,  one  can  calculate  the  end-to-end 
distance  of  the  molecule  and  the  external  force  exerted  on  it,  /. 

It  is  important  to  point  out  that,  although  this  method  has  been  used  primar- 
ily with  atomistic  models,  it  is  not  exclusive  to  molecular  dynamics  simulations.  Hozl 
et  al.  [96,  97],  for  example,  used  a similar  approach  to  calculate  the  forces  required 
to  stretch  chains  grafted  to  two  different  surfaces  (see  Figure  3.3)  from  lattice  Monte 
Carlo  simulations.  The  goal  of  their  study  was  to  examine  the  influence  of  the  wall 
and  the  influence  of  inter-chain  interactions  on  the  elasticity  of  the  tethers.  Hozl  et 
al.  fixed  the  position  of  the  external  potential,  zq  (which  is  equivalent  to  setting  the 
pulling  velocity,  v,  to  zero),  and  used  the  average  position  of  the  ‘attached’  (end) 
segment  to  obtain  the  force  from  Eqn.  3.14,  which  in  this  case  is  written  as 


f — kc[{Zgjici)  Zq\. 


(3.15) 
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Figure  3.3:  The  elasticity  of  single  chains  is  modified  by  wall/segment  interactions  and 
interchain  interactions.  For  instance,  the  ‘elasticity’  of  a collection  of  chains  grafted 
to  two  different  surfaces  is  a function  of  the  grafting  density.  For  large  separations 
{H/Rq  > 1)  the  forces  are  attractive  (except  for  very  large  grafting  densities),  whereas 
for  moderate/short  separations  {H/Rq  > 1),  the  forces  are  repulsive.  The  walls  shown 
in  this  figure  are  modelled  as  regions  in  space  where  the  energy  of  the  segments  is 
large,  i.e.,  they  are  ‘soft’  walls  [96,  97]. 
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In  order  to  model  the  upper  surface,  Hozl  et  al.  applied  a repulsive  potential  at 


position  Ze.nd  + 1.  Note  that  the  purpose  of  this  repulsive  potential  (which  moves 
according  to  the  position  of  the  end-segment)  is  to  confine  the  majority  of  the  seg- 
ments to  positions  2 < Zgnd,  and  that  the  surface  modelled  in  such  a manner  is  not 
a hard  wall.  The  results  obtained  from  this  study  show  that,  at  large  separations 
{i.e.,  H > Ro,  where  H is  the  separation  between  the  walls),  the  chains  exerted  an 
attractive  force  on  the  walls  (their  collective  elasticity  being  dependent  of  the  grafting 
density  or  number  of  chains  per  area),  whereas  for  H < Rq  the  force  was  repulsive 
(due  to  the  number  of  conformations  excluded  by  the  ‘wall’). 

The  repulsive  or  steric  forces  induced  by  confined  polymers  can  be  calculated 
by  different  methods.  The  calculation  of  a pressure  tensor  in  off-lattice  simulations 
is  a robust  technique  that  can  be  applied  to  a large  number  of  problems.  In  the  next 
section  we  discuss  the  details  of  such  a calculation. 

3.2  Confined  Polymers:  Calculation  of  the  Pressure  from  the  Virial 

The  calculation  of  polymer-induced  forces  in  off-lattice  models  is  generally 
based  on  the  calculation  of  the  pressure[16,  17,  98,  99].  The  element  of  the 
pressure  tensor  is  obtained  from  the  virial  equation,  which  relates  the  pressure  to  the 
ensemble  average  of  the  microscopic  stress  tensor: 


n 


n— 1 n 


(3.16) 


Here  n is  the  total  number  of  polymer  segments  in  the  system  (i.e.,  n = Nm,  where 
m is  the  number  of  polymer  chains  and  N is  the  chain  length),  V is  the  volume  of 
the  system,  mj  and  Vj  are  the  molecular  mass  and  velocity  of  segment  i,  respectively. 
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and  Tij  is  the  vector  between  segments  i and  j.  The  force,  fjj,  between  segments  i 
and  j is  given  by 


where  U is  the  total  interaction  potential  between  segments  i and  j.  It  includes  both 
bonded  interactions  and  non-bonded  interactions.  Bonded  interactions  are  usually 
described  by  a harmonic  potential  or  by  the  finite  extensible  nonlinear  elastic  (FENE) 
potential,  whereas  nonbonded  interactions  are  generally  modeled  by  a Lennard  Jones 
or  a Morse  potential.  The  only  restriction  that  is  implicit  to  the  above  equation  is 
that  all  the  potentials  U{rij)  should  be  continuous  at  all  separations  Uj. 

The  systems  we  usually  consider  consist  of  polymers  confined  by  two  surfaces. 
These  confining  surfaces  are  generally  located  at  2 = 0 and  dX  z — H.  The  force 
induced  by  the  polymers  is  thus  related  to  the  ^-component  of  the  pressure  tensor, 
Pzz-  For  instance,  the  interaction  between  two  polymer  brushes  is  obtained  as  [100] 


where  Pzz{H)  is  the  pressure  at  a given  separation  H,  and  Pzz{H  00)  is  the 
pressure  when  the  two  brushes  are  not  yet  in  contact  with  each  other  (see  Figure 
3.4).  In  these  cases,  one  can  drop  the  thermal  contribution  {i.e.,  the  first  term  in  the 
virial  equation)  and  get  Pzz  as 


-r.j  dUjrij) 


(3.17) 


f{H)  = [P„{H)  - P„{H  -+  °o)l , 


(3.18) 


n— 1 n 


(3.19) 


t=l  j>l 


where  L is  the  side  of  the  lattice  and  H the  separation  between  walls. 
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Figure  3.4;  The  interaction  between  two  polymer  brushes  is  determined  by  the  com- 
ponent of  the  pressure  tensor  normal  to  the  surface  {Pzz  in  this  case).  Note  that  at 
equilibrium  the  lateral  components  of  the  tensor  should  be  identical  (i.e.,  Pxx  = Pyy)- 


In  some  circumstances,  it  is  advantageous  to  calculate  a local  pressure  tensor 
in  the  z-direction  [17,  98,  99].  This  local  component  of  the  pressure  is  given  by 


Pap{z)  = p{z)kBTSap  - 


\ (ro)/9  dUjrjj) 

e{n,) 


[5{zi  -z)  + 5{zj  - z)] . 


(3.20) 


Here  is  the  component  of  the  local  pressure  tensor,  p{z)  is  the  density  of  the 
polymer  segments,  is  the  Kronecker  delta  (i.e.,  Sij  = 1 if  i = j and  % = 0 if 
i ^ j)  and  5(x)  is  the  step  function  (5(a:)  = 0 if  x < 0,  and  <5(x)  = 1 if  x > 0).? 
The  use  of  the  local  pressure  defined  in  Eqn.  3.20  presents  certain  advantages  over 
the  calculation  of  since  it  allows  one  to  examine  whether  equilibrium  has  been 
reached.  For  instance,  in  most  of  the  confined  systems,  one  expects  symmetry  parallel 
to  the  surface,  which  implies  that  P^x  = Pyy  Furthermore,  in  some  cases,  as  in 
the  study  of  polymer  layers  in  equilibrium  with  a bulk  solution,  the  local  pressure 
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Figure  3.5:  The  calculation  of  the  local  pressure  can  be  used  to  verify  whether  equi- 
librium conditions  have  been  reached.  This  figure  shows  a polymer  solution  confined 
by  two  adsorbing  walls.  If  the  conditions  in  the  middle  of  the  gap  are  similar  to  those 
of  a bulk  solution,  the  normal  and  lateral  components  of  the  pressure  tensor  should 
be  identical. 

tensor  can  be  used  to  corroborate  the  ‘bulk’  behavior  at  the  middle  of  the  gap  {i.e., 
Pxx{H/2)  = Pyy{H/2)  = Pzz{H/2))  (see  Figure  3.5).  It  is  important  to  point  out  that 
integration  of  the  local  pressure  tensor,  Pa^{z),  over  the  2— direction  recovers  the 
component  of  the  pressure  defined  in  Eqn.  3.16. 

Calculations  of  forces  from  the  pressure  tensor  have  been  used  in  a wide  variety 
of  systems,  e.g.,  single  confined  chains  [101],  confined  polymer  solutions  [17,  102], 
interaction  between  polymer  brushes  [100, 103],  etc.  The  use  of  this  approach  is  highly 
recommended,  provided  an  efficient  algorithm  is  available  to  equilibrate  the  system  in 
a reasonable  time.  There  are,  however,  situations  where  a lattice  model  is  preferred 
over  an  off-lattice  model.  For  instance,  examination  of  the  validity  of  numerical  mean- 
field  lattice  calculations  evidently  requires  the  use  of  lattice  simulations.  In  addition. 
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the  simulation  of  long  polymer  chains  in  the  presence  of  attractive  surfaces  requires 
lattice  simulations  unless  access  to  a massively  parallel  supercomputer  is  available. 

The  calculation  of  forces  from  lattice  simulations  requires  a diflferent  approach. 
Instead  of  direct  calculation  of  the  pressure  or  forces,  in  lattice  simulations  one  needs 
to  calculate  differences  in  the  free  energy  {i.e.,  Helmholtz  potential).  In  the  next 
section  we  describe  one  method,  based  on  thermodynamic  integration,  that  has  been 
used  to  estimate  the  pressure  in  a lattice.  We  then  follow  with  the  details  of  an  efficient 
technique  to  estimate  Helmholtz  potential  differences  from  Monte  Carlo  simulations. 
(This  technique,  which  is  called  the  acceptance  ratio  method,  was  actually  developed 
for  off-lattice  simulations  of  simple  fluids.) 


3.3  Estimation  of  Free  Energy  Differences  from  Lattice  Monte  Carlo  Simulations 
The  methods  described  above  to  obtain  the  pressure  from  the  virial  cannot 
be  implemented  in  a lattice,  where  the  space  is  discretized  and  the  potential  between 
segments  is  usually  given  by  a step  function  {e.g.,  a square- well  potential).  The  forces 
in  this  case  need  to  be  estimated  from  the  derivative  of  the  Helmholtz  energy  with 
respect  to  a relevant  coordinate.  Consider,  for  example,  the  case  of  polymer  chains 
conflned  by  two  walls  located  at  z = 0 and  z — H.  The  force  of  interaction  between 
the  two  surfaces,  /,  arising  from  the  presence  of  the  polymer  chains,  is  simply  the 
change  in  free  energy,  T,  as  the  surfaces  are  moved  from  d = H to  d = H — where 
d is  the  distance  between  the  surfaces  in  lattice  units. 


/ ld=:H-l/2  — 


T{H) 


(3.21) 


61 


Figure  3.6:  Compression  of  polymer  chains  by  a wall.  Introduction  of  a repulsive 
potential  and  definition  of  contact. 

The  problem  of  moving  a ‘hard’  wall  (or  in  general,  the  boundary  of  a solid  object) 
is  equivalent  to  considering  two  systems  with  identical  potential  fields  except  in  the 
layer  immediately  adjacent  to  the  moving  wall  (say,  the  ‘surface  layer’).  At  the  surface 
layer  an  external  potential,  Ugxu  ranging  from  0 to  co  is  applied,  with  the  application 
of  the  infinite  potential  signifying  the  relocation  of  the  wall  to  the  location  of  the 
original  surface  layer  (see  Figure  3.6). 

In  order  to  estimate  the  change  in  the  Helmholtz  potential  between  the  initial 
state  {d  = H)  and  the  final  state  {d  — H — 1),  Dickman  [24,  104,  105]  proposed  the 
use  of  a thermodynamic  integration  procedure,  which  we  describe  below. 

3.3.1  Repulsive  Wall  Technique:  A Thermodynamic  Integration  for  Lattice  Systems 
The  first  step  in  Dickman’s  repulsive  wall  method  is  the  definition  of  a param- 
eter A = , which  is  a parameter  that  changes  smoothly  from  the  initial  state 
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{d  = H)  to  the  final  state  {d  = H - 1).  Note  that  A = 1 corresponds  to  Uext  = 0, 
that  is,  the  upper  surface  is  at  z = if;  in  contrast,  A = 0 (i.e.,  Uext  = oo)  implies 
that  the  surface  is  at  2:  = if  - 1 (since  the  excluded  volume  criterion  prevents  the 
polymer  segments  from  occupying  z = H — 1).  The  change  in  Helmholtz  potential 
between  the  initial  and  final  states  can  be  written  as 

= Q U=i  -Q  U.o=  r , (3.22) 


where  Q is  the  partition  function  of  the  system. 

For  this  particular  application,  the  derivative  in  Eqn.  3.22  can  be  written 
in  terms  of  the  number  of  polymer  segments  in  layer  if  — 1 {i.e.,  the  number  of 
contacts,  see  also  Figure  3.6)  as  follows: 


51nQ 

dX 


(3.23) 


where  the  configurational  energy  U does  not  include  the  contribution  NcUext,  which 
is  written  separately  as  shown.  Equation  3.22  then  becomes 

= [ ^dX  . (3.24) 

Jo  X 

The  above  integral  is  evaluated  by  performing  simulations  for  different  values  of  Uextj 
and  hence,  for  different  values  of  A.  One  calculates  the  average  number  of  contacts, 
(Nc),  from  each  simulation  and  obtains  the  value  of  {Nc)/X  at  A = 0 {i.e.,  Uext  = 00) 
by  extrapolation  (see  Figure  3.7).  Equation  3.24  is  then  evaluated  by  numerical 


integration. 


63 


AF  = 

0 A 

0 A 1 

Dickman’s  Method 


Figure  3.7:  Estimation  of  free-energy  differences  based  on  a thermodynamic  integra- 
tion approach. 

One  drawback  of  Dickman’s  method  is  that  it  requires  a large  number  of 
intermediate  steps,  i.e.,  several  values  of  A,  in  order  to  obtain  a smooth  set  of  points 
which  can  be  easily  extrapolated  to  A = 0.  Dickman  suggests  that  at  least  six  values 
of  A be  used  [104].  In  the  next  section  we  present  a more  general  and  efficient  method 
to  estimate  free-energy  differences  from  Monte  Carlo  simulations.  This  method  will 
be  the  starting  point  for  a new  technique  (the  Contact- Distribution  method)  to  obtain 
forces  from  lattice  simulations,  which  is  described  in  Chapter  4. 

3.3.2  Acceptance  Ratio  Method 

The  other  technique  that  is  of  interest  here,  developed  by  Bennet  [25],  is 
called  the  acceptance-ratio  method  and  was  introduced  as  an  efficient  way  to  estimate 
free  energy  differences  from  Monte  Carlo  data  in  any  kind  of  system.  It  considers  two 
systems,  denoted  as  0 and  1,  with  configurational  energies  Uq  and  U\  obtained  from 


64 


respective  pair-potentials  Uq  and  Ui.  The  ratio  of  the  canonical  partition  functions  of 
these  systems,  which  defines  the  difference  in  Helmholtz  energy  = ln(<3o/Qi)) 
can  be  written  as 


0„ 

0i  Q,/dr''m(r'')e-«"«+“0  (to(r'‘)e-«'>)„  ’ ''  ’ 

where  w{r^)  is  a weight  function  for  the  configuration  specified  by  r^,  the  positions 
of  the  N particles  in  the  system.  The  optimum  choice  for  w(r^)  is  obtained  by 
minimizing  the  statistical  error  in  /3AT  with  respect  to  w{r^)  [25,  106].  Such  a 
procedure  leads  to 


where  C is  a constant.  Rewriting  Eqn.  3.26  in  terms  of  the  Fermi-Dirac  function 


(3.26) 


9{x) 


1 

1 -f  ’ 


(3.27) 


one  obtains 


Qo  _ (g(^o  — Ui+  C))i 


(3.28) 


Qi  {g{U,-Uo-C))o  ■ 

Bennet  [25]  has  shown  that  the  optimal  choice  of  C,  denoted  henceforth  as  C*,  cor- 


responds to 


Qo 

Qi  no’ 


(3.29) 


where  no  and  ni  are  the  numbers  of  statistically  independent  samples  taken  from 
Systems  0 and  1,  respectively.  Substituting  for  {Qo/Qi)  in  Eqn.  3.28  and  taking 
C — C*,  one  gets 


n\  no 

i:  - us + C-) = E - <’  - c')' 


Z=1 


t=l 


(3.30) 
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where  is  the  energy  obtained  from  the  ith  sample  of  System  j using  the  pair- 
potential  Ml  and,  similarly,  UqJ  is  the  energy  calculated  with  the  pair-potential  Mq. 
(Note  that  one  needs  to  run  simulations  of  both  systems  and  that  the  Boltzmann 
sampling  in  each  system  is  based  on  its  own  configurational  energy.  For  example, 
u[°i  is  the  configurational  energy  calculated  using  the  pair-potential  Ui  for  a sample 
which  is  chosen  with  a probability  proportional  to  exp  {—^Uq)  from  System  0.)  Once 
C*  is  obtained  from  Eqn.  3.30,  it  can  be  substituted  into  Eqn.  3.29  for  (Qo/Qi)  so 
that  A.F  can  be  obtained  from 

In  f = ^AJF  = ^C*  - In  — . (3.31) 

V Q 1 / ^0 

In  the  next  chapter,  we  extend  Bennet’s  method  to  obtain  forces  from  polymer 
systems  and  lattice  Monte  Carlo  simulations.  The  new  technique  is  computationally 
efficient  and  requires  less  number  of  simulations  as  compared  to  Dickman’s  method. 


CHAPTER  4 

A NEW  SIMULATION  METHOD  FOR  POLYMER-INDUCED  FORCES 


In  this  chapter^  we  introduce  a new  simulation  method,  which  we  call  the 
contact-distribution  method  (CDM),  for  the  determination  of  the  Helmholtz  potential 
for  polymer /colloid  systems  from  lattice  Monte  Carlo  simulations.  The  method  is 
general  enough  to  allow  examination  of  systems  with  chemisorbed  (Le.,  grafted)  or 
physisorbed  polymer  chains  on  surfaces  of  arbitrary  shape  and  dimensions  and  is 
computationally  very  efficient.  Our  objective  here  is  to  present  the  details  of  the 
simulation  method  and  to  illustrate  its  application  using  three  examples:  (i)  the 
compression  of  a physisorbed  layer  by  a bare  surface,  {ii)  the  interaction  between  a 
finite-sized  object  and  a single,  end-grafted  polymer  chain  {i.e.,  a polymer  mushroom) 
and  {in)  the  interaction  between  a finite-sized  object  and  a collection  of  end-grafted 
chains  (also  called  a polymer  brush).  Numerical  results  for  the  first  case  are  used  to 
discuss  certain  details  on  the  implementation  of  CDM.  The  results  on  the  compression 
of  polymer  mushrooms  and  polymer  brushes  illustrate  how  the  method  can  be  used 
to  confirm  and  extend  scaling  laws  for  forces  and  Helmholtz  potentials.  It  also  allows 
one  to  examine  conformational  transitions  and  their  relation  with  the  induced  forces. 

^ This  chapter  is  a modified  and  updated  version  of  Jimenez  and  Rajagopalan 
[15,  18]. 
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(a)  Moving  a Surface 
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(b)  Applying  a Repulsive  Potential 
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Figure  4.1:  Compression  of  polymer  chains  by  a bare  surface:  Two  equivalent  views. 

As  mentioned  in  Chapter  3,  the  methods  used  to  calculate  the  pressure  from 
the  virial  cannot  be  implemented  in  a lattice,  and  hence,  the  polymer-induced  forces 
need  to  be  estimated  from  the  derivative  of  the  Helmholtz  energy  with  respect  to  the 
relevant  spatial  variable.  Consider,  for  example,  the  case  of  polymer  chains  confined 
by  two  walls  located  at  2 = 0 and  z — H,  as  shown  in  Figure  4.1a.  The  force  of 
interaction  between  the  two  surfaces,  /,  arising  from  the  presence  of  the  polymer 
chains  is  simply  the  change  in  free  energy,  as  the  surfaces  are  moved  from  d — H 
to  d = H — 1,  where  d is  the  distance  between  the  surfaces  in  lattice  units. 

/ lfcH-i/2  = -%  = na  - 1)  - HH)  (4.1) 

The  difference  in  Helmholtz  potential  in  Eqn.  4.1  is  equivalent  to  that  between  two 
systems  with  identical  potential  fields  except  in  the  layer  immediately  adjacent  to 
the  moving  wall  (say,  the  ‘surface  layer’).  At  the  surface  layer  an  external  potential. 
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Uexu  ranging  from  0 to  oo  is  applied,  with  the  application  of  the  infinite  potential 
signifying  the  relocation  of  the  wall  to  the  location  of  the  original  surface  layer  (see 
Figure  4.1b),  i.e., 

dT 

f lfc«-l/2  = - grf  = = °°)  - = “)  (‘‘•2) 

In  Chapter  3,  Section  3.3  we  showed  that  the  change  in  the  Helmholtz  po- 
tential in  Eqn.  4.2  could  be  obtained  using  a thermodynamic  integration  procedure 
[24,  104,  105].  In  what  follows  we  show  that  an  alternative  procedure,  based  on  the 
acceptance-ratio  technique  [25]  (see  also  Chapter  3,  Section  3.3),  results  in  a more 
efficient  method  to  estimate  Helmholtz  energy  differences.  This  procedure  is  what  we 
call  the  contact-distribution  method,  CDM. 

4.1  The  New  Technique:  Contact-Distribution  Method 
We  begin  the  description  of  CDM  with  the  definition  of  two  systems  labelled 
0 and  1.  System  0 is  a lattice  of  sides  x Ly  and  height  that  contains  a certain 
number  of  polymer  chains  and  a ‘hard-body’  {i.e.,  an  object  that  occupies  a certain 
number  of  lattice  sites).  Our  goal  is  to  obtain  the  difference  in  Helmholtz  energy 
that  occurs  when  the  coordinates  of  the  ‘hard-body’  are  moved  one  lattice  unit  in 
any  given  direction.  In  this  context.  System  1 is  identical  in  every  respect  to  System 
0 but  has  a repulsive  potential  Uext  = V (with  V > 0)  at  the  moving  boundary 
of  the  ‘hard-body’.  All  other  interaction  energies  {e.g.,  polymer  segment/segment 
interaction  energy)  remain  the  same  in  both  systems.  The  ‘hard  body’  alluded  to 
above  could  be  a solid  surface,  as  in  Figure  4.1,  or  any  arbitrary  object  as  shown  in 
Figure  4.2.  With  this  definition,  the  difference  between  the  configurational  energies 
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Figure  4.2:  Compression  of  polymer  chains  by  a finite-sized,  arbitrarily-shaped  object, 
of  the  two  systems  is  given  by 

Ui-Uo=^N,V,  (4.3) 

where  Nc  is  the  number  of  contacts  {i.e.,  the  number  of  segments  of  the  polymer 
chains  in  the  layer  next  to  the  moving  boundary  of  the  object). 

Following  Bennet’s  approach  [25,  106],  the  ratio  of  the  canonical  partition 
functions  of  Systems  0 and  1,  which  defines  the  difference  in  Helmholtz  energy  /3AT  = 
\n{Qo/Qi),  can  be  written  as 


where  (7  is  a constant.  Rewriting  Eqn.  4.4  in  terms  of  the  Fermi-Dirac  function 

1 


(4.4) 


g{^)  = 


1 e^x  ’ 


(4.5) 
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and  substituting  Eqn.  4.3  in  Eqn.  4.4,  one  gets 


where  Nc^  is  the  number  of  contacts  obtained  in  a simulation  of  System  j. 

Bennet  [25]  has  shown  that  the  optimal  choice  of  C,  denoted  henceforth  as 
C*,  corresponds  to 

Qo  ^1 

Qi  no’ 


where  no  and  ni  are  the  numbers  of  statistically  independent  samples  taken  from 
Systems  0 and  1,  respectively.  Substituting  for  {Qo/Qi)  in  Eqn.  4.6  and  taking 
C = C*,  one  gets 

ni  no 

'£g{C-  - - C*).  (4.8) 

t=l  i=l 

The  term  in  Eqn.  4.8  indicates  the  number  of  contacts  drawn  from  sample  i in 
System  j.  The  only  unknown  in  the  previous  equation  is  C*.  The  number  of  terms 
in  this  equation  is,  however,  rii  + no,  which  is  typically  of  the  order  of  10^  or  10^.  It 
is  therefore  useful  to  rewrite  Eqn.  4.8  in  a more  convenient  way,  i.e.,  not  in  terms 
of  the  number  of  contacts  in  each  sample  but  in  terms  of  the  probability  of  having  a 
given  number  of  contacts.  Let  be  the  probability  of  the  polymer  chains  having 
k contacts  with  the  hard  body  in  System  1 (in  the  system  with  repulsive  potential 
Uext  = at  the  boundary)  and  the  corresponding  probability  in  System  0 (in 
the  system  with  no  repulsive  potential).  Assuming  no  = ni,  one  can  rewrite  Eqn.  4.8 
as 

^max  ^max 

£ P<‘>9(C<  - AV)  = £ - C-)  . 

fe=0  fc=0 


(4.9) 
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where  k,  the  number  of  contacts  with  the  boundary,  goes  from  0 to  kmax,  with  Aimax 
being  the  largest  possible  number  of  contacts  (determined  by  the  number  of  lattice 
sites  available  at  the  boundary).  Using  the  Fermi-Dirac  equation  for  the  function 
g{C*  — kV)  one  gets 


^(1)  _ pW^0{C*—kV) 


= 0. 


(4.10) 


lpe^{c*-kv) 

fc=0 

Equation  4.10  is  a nonlinear  equation  with  one  unknown,  namely,  C*,  which  can  be 
obtained  once  and  are  determined  using  simulations.  From  Eqn.  4.7  it  then 
follows  that 


AP  = P{Ue,t  = V)-  P{Ue.t  = 0)=)3C'*.  (4.11) 


It  is  important  to  point  out  that,  following  Bennet’s  arguments  [25],  the  statistical 
error  in  in  Eqn.  4.11  depends  on  the  extent  of  overlap  of  the  configurational 
spaces  sampled  in  the  two  systems,  which  in  turn  depends  on  the  extent  of  overlap 
of  the  distributions  of  contacts  and  We  will  examine  this  issue  further  in 
Section  4.2. 

At  this  point  we  should  also  remember  that  we  are  interested  in  the  change 
in  the  Helmholtz  potential  that  occurs  when  Uext  varies  from  0 to  oo  (see  Eqn.  4.2). 
It  is  therefore  useful  to  examine  what  happens  to  Eqn.  4.10  in  the  limit  of  Uext  oo. 
In  this  case,  the  polymer  segments  are  excluded  from  the  layer  next  to  the  moving 
object,  i.e.,  = 1,  P^^^  = 0 and  = 0 for  any  A;  > 1.  Equation  4.10 

therefore  reduces  to 


1 - 
1 + 


= 0 


(4.12) 
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which  implies  that 


PC’  = = oo)  - = 0)  = - lnP(i">. 


(4.13) 


As  already  noted  by  Dickman  [24],  Eqn.  4.13  is  well  known  and  shows  that  if 
the  probability  of  zero  contacts  at  the  boundary  of  the  hard-body  is  known,  one  can 
get  the  interaction  force  (as  the  solid  object  is  moved  from  z — Htoz  — H — \)diS 


Although  very  simple,  this  equation  is  not  useful  for  cases  where  the  density  of  the 
segments  is  high  (i.e.,  for  high  compressions)  since  the  probability  of  zero  contacts 
is  very  small  and,  therefore,  hard  to  determine  accurately.  It  is  in  such  cases  where 


needed.  The  idea  is  to  use  a potential  Vi  that  guarantees  an  overlap  of  the  resulting 
contact  distribution  with  the  distribution  corresponding  to  System  0 (a  system  with 
no  repulsive  potential),  and  at  the  same  time,  that  allows  for  an  accurate  determina- 
tion of  the  probability  of  having  zero  contacts.  There  are  situations,  however,  where 
the  two  conditions  described  above  cannot  be  met  simultaneously  {e.g.,  in  very  dense 
systems  or  at  large  compressions).  It  is  in  these  situations  where  it  may  be  necessary 
to  introduce  additional  intermediate  steps  with  repulsive  potentials  V2,  V3,  etc.,  and 
decompose  the  total  change  in  the  Helmholtz  potential  into  the  appropriate  number 
of  contributions  (see  Section  4.2  for  an  example). 


(4.14) 


the  introduction  of  a finite  repulsive  potential  at  the  boundary  of  the  hard-body  is 


For  instance,  let  us  consider  that  in  a particular  case  three  simulations  are 
performed  for  Uext  = {0,  Vi,  V2},  which  denote  Systems  0,  1,  and  2,  respectively.  Let 
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us  also  consider  that  the  probability  of  having  zero  contacts  is  essentially  zero  for 
Systems  0 and  1 and  is  determined  accurately  for  System  2,  i.e.,  = 0 and 

= p,  where  p > 0.  The  change  in  the  Helmholtz  potential  that  occurs  as  Uext  is 
varied  from  0 to  oo  is  then  given  by  the  following  sum; 

P{Uext  = oo)  - P{Uext  = 0)  = APq^i  + A ~ In  pS"\  (4.15) 


where  APq-H)  the  Helmholtz  potential  difference  between  Systems  1 and  0,  is  ob- 
tained from 


fcmax  p(l)  p(0)p(A^O->l-fcVi) 

E^k  ^ _ n 

1J_  p(A:Fo_n-fcVi) 

fc=0 


(4.16) 


and  AP’i_).2,  the  corresponding  free  energy  difference  between  Systems  2 and  1,  is 
obtained  from 

femax  p(2)  _ p(l)  [A.Fi_>2-fc(V2-Vi)] 

ii-_ = 0 (4.17) 

1 + e[A:Fi_2-A;(V2-Vi)]  ''  > 

k=0 

The  last  term  in  Eqn.  4.15  represents  the  difference  in  the  Helmholtz  potential 
between  System  2 and  a system  where  Uext  — oo. 

The  need  for  additional  steps,  however,  has  to  be  tested  for  each  particu- 
lar case.  In  the  next  section  we  present  numerical  results  for  the  compression  of  a 
physisorbed  polymer  layer  by  a bare  surface.  Our  goal  here  is  not  the  physics  of 
the  problem  (which  is  discussed  in  length  by  de  Joannis  et  al.  [107]),  but  to  sketch 
a strategy  for  the  implementation  of  CDM  {e.g.,  for  choosing  the  number  and  val- 
ues of  the  repulsive  potentials),  and  to  compare  CDM  with  a method  based  on  a 
thermodynamic  integration,  the  so-called  Dickman’s  method. 
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4.2  Implementation  of  CDM:  Compression  of  a Polymer  Layer  with  a Bare  Surface 

The  example  we  consider  here  is  the  compression  of  physisorbed  polymers  by 
a bare  wall,  as  shown  schematically  in  Figure  4.3.  The  lattice  simulations  used  for 
this  purpose  employ  the  simple  cubic  lattice  and  the  configurational  bias  method, 
both  of  which  are  described  in  Appendix  A.  A lattice  size  of  40  x 40  x is  used  with 
impenetrable  walls  at  2 = 0 and  z = H and  with  periodic  boundary  conditions  in 
the  X-  and  y-directions.  The  simulations  consider  m chains  with  N segments  each. 
The  segments  interact  only  with  the  surface  located  at  2:  = 0.  The  upper  surface, 
located  a,t  z = H,  is  said  to  be  ‘athermal’.  The  discrete  potential,  Uads,  between  the 
segments  and  the  bottom  surface  is  effective  at  layer  1.  The  strength  of  this  potential 
is  Uads  = —1.0  ksT,  where  ks  is  the  Boltzmann  constant  and  T the  temperature  of  the 
system.  The  problem  of  confining  a physisorbed  layer  with  a bare  surface  is  an  ideal 
system  to  examine  the  interplay  between  the  loops,  tails  and  trains  that  constitute  a 
physisorbed  layer  (see  Chapter  2).  An  extensive  analysis  of  this  problem  was  recently 
presented  by  de  Joannis  et  al.  [107].  In  what  follows  we  focus  on  numerical  results 
obtained  for  m = 8 and  N = 200.  Our  goal  here  is  to  illustrate  the  implementation 
of  CDM  as  a method  to  calculate  differences  in  the  Helmholtz  potential  and  forces 
induced  by  polymers  under  compression. 

4.2.1  Implementation  of  CDM 

In  Table  4.1  below  we  show  the  probability,  Pq,  of  having  zero  contacts  be- 
tween the  polymer  chains  and  the  surface  located  at  z = H.  The  probability  Pq  is 
related  to  the  number  of  samples  for  which  the  surface  layer,  located  at  2 = if  — 1,  is 
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‘Athermal’ 

wall 


Figure  4.3:  Compression  of  physisorbed  polymers  by  a bare,  ‘athermal’  surface. 

empty  (we  denote  this  number  of  samples  as  Uempty)-  The  probability  of  zero  contacts 
is  then  given  as 

Po  = (4.18) 

n 

where  n is  the  total  number  of  samples  drawn  from  the  simulations.  If  Pq  can  be 
determined  accurately,  the  force  of  compression  is  obtained  directly  from  Eqn.  4.14. 
The  accuracy  of  Pq  is  strongly  dependent  on  the  total  number  of  samples,  n,  used  in 
the  simulation.  For  instance,  the  typical  value  of  n used  for  the  simulations  reported 
in  Table  4.1  below  was  n « 10^-10^,  which  implies  that  any  value  of  Pq  ^ 10“^  may 
have  a considerable  uncertainty. 

Based  on  our  previous  comments  on  the  uncertainty  of  Po,  one  can  argue 
that,  for  i7  > 7 (see  Table  4.1),  one  can  determine  the  force  of  compression  directly 
from  Eqn.  4.14,  i.e.,  from  a single  simulation.  For  H = 6 and  Ff  = 5,  it  is  advisable 
to  perform  an  additional  simulation  with  a repulsive  potential  Vi  at  layer  z = H — 1. 


76 


Table  4.1:  Compression  of  a physisorbed  layer  by  a bare  wall. 


Separation  between  Walls  Probability  of  Zero  Contacts 
H/a  Po 

9 

0.535 

8 

0.450 

7 

0.187 

6 

0.049 

5 

0.0012 

4 

— 

This  additional  simulation  would  allow  one  to  determine  a more  accurate  value  of  Pq 
in  this  case).  The  difference  in  Helmholtz  potential  between  the  system  where 
U^xt  = Vi  and  that  where  Uext  = 0 (i.e.,  the  ‘athermal’  case)  would  be  obtained  using 
Eqn.  4.10  or  its  equivalent,  Eqn.  4.16. 

The  more  interesting  case,  which  we  analyze  in  detail  here,  is  that  of  = 4, 
for  which  the  value  of  Pq  obtained  from  the  simulations  was  zero.  The  need  to 
perform  additional  simulations  in  this  case  is  evident.  We  first  present,  in  Figure  4.4, 
the  distribution  of  contacts  obtained  for  the  athermal  wall  (denoted  here  as  System 

0)  and  for  a wall  with  a repulsive  potential  Uext  = 0.2  ksT  (denoted  here  as  System 

1) .  As  shown  in  the  figure,  there  is  an  evident  overlap  between  the  two  distributions. 
The  probability  of  having  zero  contacts  in  System  1,  Pq^\  is  however,  very  small. 
Nevertheless,  it  is  now  possible  to  make  an  estimate  of  the  force  of  compression.  In 
order  to  illustrate  this  calculation,  we  define  here  E(C’)  as  the  sum  in  Eqn.  4.10 
which,  for  this  particular  case,  is  given  as 

p(l)  p(0)  (C-0.2A:) 

^ 

1 I g(C-0.2fc) 

fc=0 


(4.19) 
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Figure  4.4:  Distribution  of  contacts  between  polymer  segments  and  a wall.  The  filled 
diamonds  correspond  to  the  system  with  no  repulsive  potential  {i.e.,  System  0)  and 
the  open  circles  correspond  to  a system  with  a repulsive  potential  Uext  = 0.2  ksT. 

The  value  of  C for  which  S = 0 (i.e.,  the  so-called  C*)  determines  the  differ- 
ence in  the  Helmholtz  potential  between  Systems  0 and  1.  It  is  evident,  from  Figure 
4.5,  that  there  is  an  easily  obtainable  single  solution  for  Eqn.  4.19.  The  value  of  C* 
obtained  in  this  manner  and  the  value  of  Pq  obtained  in  System  1 are  then  used  to 
estimate  the  force  of  compression  when  the  wall  is  moved  from  H = 4 to  H = 3 (see 
the  inset  of  Figure  4.5).  The  result  obtained  in  this  case  is 

f{H  = 3.5)  = 16.89  ksT/a, 
where  a represents  the  lattice  unit. 

It  is  interesting  to  compare  what  would  happen  if  we  had  not  chosen  a re- 
pulsive potential  Uext  = 0.2  ksT  but  a larger  one.  In  Figure  4.6  we  present  the 
distribution  of  contacts  obtained  for  the  athermal  system  (where  Uext  = 0)  and  for  a 
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Figure  4.5:  Calculation  of  the  difference  in  Helmholtz  potential  between  System  1 
and  System  0.  The  value  of  is  taken  from  Figure  4.4. 

system  with  a repulsive  potential  Uext  ~ 0-8  ksT  (denoted  here  as  System  2).  The  use 
of  a larger  repulsive  potential  allows  one  to  estimate  more  accurately  the  probability 
of  having  zero  contacts,  but,  as  shown  in  the  figure,  the  larger  potential  decreases  the 
overlap  between  the  two  distributions. 

The  effect  of  the  low  degree  of  overlapping  can  be  seen  in  Figure  4.7.  In  this 
case,  the  crossover  of  H(C)  from  positive  to  negative  values  is  not  as  clear  as  the  one 
shown  in  Figure  4.5.  The  increased  difficulty  in  obtaining  the  value  of  C*  offers  a 
clear  evidence  of  how  the  degree  of  overlap  of  the  distributions  affects  the  uncertainty 
of  the  solution.  The  force  of  compression  obtained  in  this  case  is 

f{H  = 3.5)  = 16.8  ksT/a. 

The  most  appropriate  procedure  for  this  particular  case  is  to  use  both  the 
distribution  of  contacts  obtained  for  U^xt  = 0.2  ksT  and  the  one  obtained  for  Uext  = 
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Figure  4.6:  Distribution  of  contacts  between  polymer  segments  and  a wall.  The  filled 
diamonds  correspond  to  the  system  with  no  repulsive  potential  (i.e.,  System  0)  and 
the  squares  correspond  to  a system  with  a repulsive  potential  Uext  = 0.8  ksT. 
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Figure  4.7:  Calculation  of  the  difference  in  Helmholtz  potential  between  System  2 
and  System  0. 
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Figure  4.8:  Distribution  of  contacts  between  the  polymer  segments  and  the  upper 
surface  obtained  from  simulations  of  System  0 (with  no  repulsive  potential),  System 
l(with  repulsive  potential  C/gif  = 0.2  ksT)  and  System  2 (with  repulsive  potential 
Uext  = 0.8  ksT). 

0.8  ksT,  as  shown  in  Figure  4.8.  Such  a procedure  provides  a good  overlapping  of 
the  distribution  of  contacts  and  an  accurate  value  of  the  probability  of  zero  contacts. 
The  force  of  compression  in  this  case  is  estimated  as 


f{H  = 3.5)  = 17.1  keT/a. 


4.2.2  Comparison  with  Thermodynamic  Integration 

It  is  now  illustrative  to  examine  the  implementation  of  the  method  developed 
by  Dickman  (see  Chapter  3,  Section  3.3)  and  to  compare  the  results  obtained  from 
CDM  to  those  obtained  using  the  thermodynamic  integration  approach.  Dickman’s 
procedure  uses  only  the  average  number  of  contacts,  (iVc),  as  opposed  to  the  entire 
distribution  used  in  CDM.  The  values  of  (Nc)  for  Systems  0,  1 and  2 are  obtained  from 
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Table  4.2:  Compression  of  a physisorbed  layer  by  a bare  wall:  Thermodynamic  Inte- 
gration. 


System 

Ue.t(i/) 

A 

{Nc) 

0 

0 

1 

50.00 

1 

0.2 

0.819 

28.10 

2 

0.8 

0.449 

4.63 

3* 

1.2 

0.301 

1.62 

4* 

1.4 

0.246 

0.95 

5* 

2.0 

0.135 

0.32 

* The  simulations  corresponding  to  Systems  3 -5  were  needed  to  properly 

extrapolate  (Nc)  /A  to  A — > 0. 

the  distributions  shown  in  Figure  4.8.  These  values,  and  additional  results  obtained 
for  other  repulsive  potentials,  are  shown  in  Table  4.2. 

Following  Dickman’s  procedure,  we  present,  in  Figure  4.9,  a plot  of  (Nc)  / A 
as  a function  of  A,  using  the  information  obtained  for  Systems  0,  1 and  2.  The 
integral  f (Nc)  /\  dX  shown  in  this  figure  is  related  to  the  difference  in  Helmholtz 
potential.  As  we  mentioned  in  Chapter  3,  one  of  the  main  disadvantages  with  this 
procedure  is  the  need  to  use  several  simulations  in  order  to:  (i)  fit  the  points  to 
an  empirical  equation  (which  may  be  needed  to  perform  the  integration),  and  (ii) 
make  an  adequate  extrapolation  of  (Nc)  /A  to  A = 0.  In  this  particular  example,  it  is 
evident  that  with  only  the  first  two  points  («.e.,  the  simulations  of  Systems  0 and  1) 
it  would  be  impossible  to  make  an  estimate  of  the  force  of  compression.  The  addition 
of  a third  point  (the  one  corresponding  to  the  simulation  of  System  2)  allows  one  to 
make  an  extrapolation  to  A = 0,  but  the  uncertainty  of  the  procedure  is  quite  large. 
The  result  obtained  from  the  information  contained  in  the  figure  leads  to 


f{H  — 3.5)  « 18  ksT/a. 
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Figure  4.9:  Compression  of  polymer  chains  by  a bare  surface:  A thermodynamic 
integration  approach.  The  calculation  shown  in  this  figure  is  an  estimate  based  on 
three  simulations.  The  open  circle  corresponds  to  the  extrapolated  value  at  A = 0. 
Once  additional  simulations  are  included,  the  thermodynamic  integration  approach 
yields  the  same  results  as  CDM. 

In  order  to  make  a better  estimate  of  the  force  of  compression  one  would 
need  additional  simulations.  For  instance,  Dickman  suggests  the  use  of  at  least  six 
simulations  [24,  104,  105].  The  result  obtained  when  one  uses  all  the  values  shown  in 
Table  4.2  is  f{H  = 3.5)  w 17.2  hgT/a.  This  effort  is,  however,  both  time  consuming 
and  unnecessary.  Even  if  one  uses  a thermodynamic  integration  approach,  there  is 
no  need  to  make  any  extrapolations  since  the  area  between  A = 0 and  the  lowest 
value  of  A (i.e.,  the  difference  in  Helmholtz  potential  between  the  system  with  largest 
repulsive  potential  used  in  the  simulations  and  a system  with  Uext  — oo)  can  be 
obtained  from  the  corresponding  probability  of  zero  contacts. 
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Figure  4.10:  Difference  in  Helmholtz  potential  between  a system  confined  by  a bare, 
‘athermal’  surface  and  a system  confined  by  an  adsorptive  wall. 

4.2.3  Adsorptive  Walls 

To  further  illustrate  the  implementation  of  CDM  and  the  comparison  with 
the  thermodynamic  integration  approach,  we  now  present  the  Helmholtz  potential 
difference  between  a system  where  the  upper  wall  is  athermal  (i.e.,  our  so-called 
System  0)  and  another  system  where  the  upper  wall  is  adsorptive,  i.e.,  Uext  = — 1-0 
fcgT.  These  systems  are  shown  schematically  in  Figure  4.10.  The  reason  we  are 
interested  in  the  difference  in  Helmholtz  potential  between  such  systems  is  made 
clear  in  Section  4.4.3,  where  we  explain  how  to  use  CDM  to  estimate  the  force  of 
interaction  between  two  physisorbed  layers. 

In  Figure  4.11  we  present  the  distribution  of  contacts  obtained  for  both  the 
athermal  and  the  adsorptive  walls.  As  evident  from  the  figure,  the  two  systems  are 
very  different  and  no  overlap  is  obtained  in  the  distribution  of  contacts. 
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Figure  4.11:  Distribution  of  contacts  between  polymer  segments  and  a wall.  The 
filled  symbols  correspond  to  the  system  with  no  repulsive  potential  {i.e.,  System  0) 
and  the  open  circles  correspond  to  a system  with  an  adsorptive  potential  Uext  = ~ 1-0 
ksT. 


It  is  illustrative  to  present  a plot  of  E!  as  a function  of  C for  this  case.  This 
plot  is  shown  in  Figure  4.12.  Although  there  is  only  one  single  value  of  C for  which 
E = 0,  the  uncertainty  in  this  value  is  very  large. 

Additional  simulations  at  intermediate  values  of  Uext  between  0 and  —1.0 
ksT  are  therefore  required  to  make  a good  estimate  of  the  difference  in  the  Helmholtz 
potential  between  the  adsorptive  and  the  athermal  wall,  AT  — TiJJext  = — 1.0  ksT)  — 
T{Uext  = 0).  Figure  4.13  shows  the  distribution  of  contacts  for  a large  number  of  these 
intermediate  potentials.  The  calculation  of  AT  in  this  case  can  be  made  using  either 
CDM  or  using  a thermodynamic  integration  approach.  The  results  obtained  by  these 
two  methods  are  identical,  CDM  yielding  AT  = 274.8  ksT  and  the  thermodynamic 
integration  yielding  AT  = 274.9  ksT. 
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Figure  4.12:  Calculation  of  the  difference  in  Helmholtz  potential  between  a system 
with  an  adsorptive  wall  and  a system  with  an  athermal  wall.  The  inset  in  the  figure 
shows  the  value  of  C at  which  H = 0. 
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Figure  4.13:  Distribution  of  contacts  between  polymer  segments  and  a wall  obtained 
for  a large  number  of  intermediate  potentials  between  Uext  = 0 and  Uext  = — 1-0  ksT. 
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The  numerical  calculations  presented  in  this  section  were  intended  to  give  a 
better  insight  on:  {i)  how  to  decide  whether  additional  simulations  are  necessary,  and 
(a)  how  to  choose  the  number  and  values  for  the  additional  external  potentials  that 
one  should  use.  At  the  same  time,  we  have  illustrated  that  CDM  is  more  efficient 
than  Dickman’s  technique,  which,  to  our  knowledge,  constitutes  the  only  method 
that  was  available  until  CDM  was  introduced  for  calculation  of  polymer  induced 
forces  from  lattice  Monte  Carlo  simulations.  In  what  follows,  our  task  is  to  show 
that  the  calculation  of  forces  from  lattice  simulations  gives  important  information 
that  can  be  used  to  examine  the  behavior  of  polymers  at  interfaces  or  the  behavior 
of  polymers  in  confined  spaces. 

4.3  Illustrative  Applications:  Compression  of  End-Grafted  Chains  by  Finite-Sized 

Objects 

The  deformation  of  polymer  layers  by  finite-sized  objects  (such  as  the  tip 
of  an  Atomic  Force  Microscope,  AFM)  has  received  considerable  attention  in  recent 
years  [108-114].  In  this  section  we  illustrate  the  use  of  CDM  in  two  different  prob- 
lems, namely,  the  interaction  of  an  AFM  tip  with  a single  end-grafted  chain  (the 
so-called  polymer  mushroom),  and  the  interaction  of  an  AFM  tip  with  a collection 
of  end-grafted  chains  (also  called  polymer  brush).  These  two  problems  are  shown 
schematically  in  Figure  4.14.  The  results  on  the  compression  of  polymer  mushrooms 
and  polymer  brushes  illustrate  how  CDM  can  be  used  to  confirm  and  extend  scaling 
laws  for  forces  and  Helmholtz  potentials.  It  also  allows  one  to  examine  the  effects 
of  the  shapes  and  sizes  of  the  compressing  objects  and  to  examine  conformational 


transitions  in  the  polymer  chains. 
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Figure  4.14:  Compression  of  end-grafted  chains  by  finite-sized  objects. 

4.3.1  Single  End-Grafted  Chain  Under  Compression 

We  first  focus  on  the  compression  of  grafted  polymers  in  the  so-called  mush- 
room regime.  While  from  a practical  perspective  this  case  is  of  importance  in  some 
biological  systems  [115,  116],  it  is  also  of  great  interest  from  a more  fundamental 
standpoint  due  to  the  novel  and  interesting  conformational  transitions^  that  could 
occur  in  a polymer  chain  under  compression  [108,  109,  110].  There  are  a number 
of  issues  that  influence  the  force  of  compression  and  the  conformations  available  to 
the  polymer  chain;  these  include,  for  example,  the  solvent  quality  and  the  shape 
and  alignment  of  the  tip  (relative  to  the  grafting  point).  Theoretical  predictions  can 
be  developed  in  the  simpler  cases  {e.g.,  compression  by  flat-bottom  cylindrical  tips) 


^ As  it  will  be  discussed  shortly,  end-grafted  chains  undergo  a conformational  tran- 
sition, from  a ‘confined’  state  to  an  ‘escaped’  state,  under  compression.  Although  this 
transition  has  been  referred  to  as  ‘phase’  transition  in  earlier  studies  [108,  109,  110], 
we  prefer  to  use  the  term  conformational  transitions  to  avoid  confusion. 
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Confined  Chain  Escaped  Chain 


Figure  4.15:  Conformational  transition  from  a ‘confined’  state  to  an  ‘escaped’  state. 

but  cannot  be  extended  to  more  complicated  ones.  The  simulation  method  we  have 
developed  offers  a means  to  study  the  latter  cases.  Our  objectives  here  are  to  ex- 
amine the  accuracy  of  the  scaling  relations  and  the  nature  of  the  above-mentioned 
conformational  transition. 

The  model  case  of  a cylinder  with  a fiat  bottom  compressing  a polymer 
mushroom  has  been  studied  using  scaling  arguments  in  both  good  [108,  109]  and 
theta  solvents  [110].  As  mentioned  earlier,  these  studies  identify  a conformational 
transition  upon  compression  (see  Figure  4.15).  In  particular,  the  chain  enters  a 
metastable  state  as  it  is  compressed  by  the  fiat  bottom  of  a cylinder  of  radius  Rt  and 
will  undergo  a conformational  transition  from  a “confined”  state  to  an  “escaped”  state 
as  the  distance  H between  the  bottom  of  the  cylinder  and  the  surface  is  decreased. 
The  use  of  scaling  arguments  allows  one  to  obtain  simple  relations  for  the  Helmholtz 
potential  T and  the  force  of  compression  /;  for  example,  for  theta  conditions,  Guffond 
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Figure  4.16:  Scaling  predictions  for  the  Helmholtz  potential  and  the  induced-force  of 
a single  end-grafted  chain  under  compression  by  a finite-sized  object. 


et  al.  [110]  have  shown  that 


^ :onflned  — ^1' 


H2  ' 


/confined  — ‘lC\ 


Na^ 

m ’ 


(4.20) 


and 


^ sscaped  — 2c2 


H' 


(4.21) 


where  a is  the  monomer  size,  N the  length  or  degree  of  polymerization  of  the  chain, 
and  Cl  and  C2  are  proportionality  constants. 

The  transition  from  the  confined  to  the  escaped  state  occurs  at  a distance 
H* , the  crossover  point  between  the  free  energies  of  the  two  states, 

* _ Cl  Na^ 

02.  2i?t 


H 


(4.22) 
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leading  to  a sudden  change  in  the  end-to-end  distance,  and  hence,  in  the  force  profile 
(see  Figure  4.16).  In  reality,  however,  the  transition  will  be  less  drastic  as  the  confor- 
mation of  the  chains  fluctuates  between  the  two  states.  One  would  expect  the  extent 
of  the  ‘noise’  in  the  force  in  the  transition  region  to  depend  on  the  free  energy  barrier 
between  the  two  states,  the  length  of  the  chain,  and  the  temperature  (or  the  qual- 
ity) of  the  solvent.  Scaling  relations  cannot,  of  course,  examine  the  effect  of  thermal 
fluctuations  on  this  transition,  but  simulations  can. 

The  simulations  presented  here  employ  the  three-dimensional  bond-fluctuation 
model  described  in  Appendix  A.  A lattice  size  of  200  x 200  x 200  is  used  with  im- 
penetrable walls  at  ^ = 0 and  z = 200  and  periodic  boundary  conditions  in  the  x- 
and  ^/-directions.  We  consider  a single  chain  grafted  to  the  center  of  the  lower  wall 
(i.e.,  at  X = 100,  y = 100,  z = 1)  and  allow  the  non-bonded  segments  to  interact 
with  each  other  through  a square-well  potential  with  an  interaction  parameter  given 
by  XcI^bT  = —0.53,  where  fcjg  is  the  Boltzmann  constant.  This  interaction  gives  rise 
to  a behavior  characteristic  of  theta  conditions  [117].  The  shape  and  the  size  of  the 
tip  and  other  conditions  used  in  the  simulations  are  summarized  in  Table  4.3.  For 
each  set  of  parameters,  several  simulations  are  run  for  different  values  of  H (with 
at  least  one  run  for  each  H).  In  each  case,  the  system  is  first  equilibrated  and  then 
sampled  over  1 - 3 x 10^  Monte  Carlo  steps  to  obtain  the  probability  of  the  chain 
having  contacts  with  the  tip,  from  which  the  Helmholtz  potentials  and  the  forces  can 
be  obtained. 

We  first  examine  the  accuracy  of  the  scaling  laws  for  the  Helmholtz  potential 
and  the  nature  of  the  conformational  transition.  Figure  4.17  shows  good  agreement 


91 


Table  4.3:  Compression  of  a polymer  mushroom:  Simulation  parameters. 


a 


Set 

Chain 

Characteristics  of  the  AFM  Tip 

Quality 

Number 

Length 

Shape 

Bottom 

Size^ 

Solvent 

I 

100 

Monolith 

Flat 

Lt  = 120 

Theta 

II 

60 

Monolith 

Flat 

Lt  = 60 

Theta 

III 

60 

Monolith*^ 

Flat 

Lt  = 60 

Theta 

IV 

100 

Monolith 

Flat 

Lt  = 100 

Theta 

V 

100 

Cylinder 

Flat 

Rt  = 50 

Theta 

VI 

100 

Cylinder 

Hemispherical 

Rt,  Rc  = 50 

Theta 

VII 

100 

Monolith 

Flat 

Lt  = 120 

Good 

Size  given  in  lattice  units.  Lt  is  the  side  of  the  monolith  tip,  Rt  is  the  radius 
the  cylindrical  tip  and  Rc  is  the  curvature  radius  at  the  bottom  of  the  tip. 

**  The  monolith  in  case  III  is  not  centered  over  the  grafting  point. 


of 


between  our  simulations  and  the  scaling  relations,  i.e.,  Eqns.  4.20  and  4.21.  (The 
proportionality  constants  Ci  and  02  in  the  equations  are  treated  as  adjustable  pa- 
rameters to  obtain  the  dashed  lines.)  Moreover,  the  value  of  H*  predicted  by  Eqn. 
4.22  from  ci,  C2  and  the  value  of  Na^  obtained  from  the  simulations  (see  the  figure 
caption)  is  roughly  about  8.2  and  falls  within  the  transition  region,  thus  confirming 
the  consistency  of  the  scaling  relations. 

Figure  4.18  shows  that  the  slope  of  the  2:-component  of  the  radius  of  gyration, 
Rgz,  with  respect  to  H changes  sign  from  positive  to  negative  at  the  onset  of  escape 
of  the  chain.  Previously,  the  variation  of  the  end-to-end  distance  has  been  used 
[108,  109,  110]  as  an  indication  of  the  escape  of  the  chain  from  under  the  tip.  We 
find,  however,  that  Rgz,  rather  than  the  end-to-end  distance,  is  a better  indicator  of 
the  point  at  which  the  escape  probability  becomes  significant. 

One  would  expect  the  size,  shape  and  alignment  of  the  tip  (relative  to  the 
grafting  point)  to  influence  the  chain  conformation  and  the  transition,  and  this  is 
evident  from  Figures  4.19  and  4.20.  These  figures  shows  that  the  distance  at  which 
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Figure  4.17:  Force  profile  when  compressing  a 100-segment  chain  with  a flat  bottom 
monolith  {Lt  = 120).  To  fit  Eqns.  4.20  and  4.21,  we  use  ~ 27  (average 

end-to-end  distance  for  the  unperturbed  chain  obtained  from  the  simulations)  and 
z=  Lt/2  = 60.  The  values  of  Ciand  C2  obtained  from  fitting  the  simulation  data  are 
1.9  and  1.4,  respectively. 


Tip/Surface  Distance  ( H/a ) 


Figure  4.18:  The  z-component  of  the  radius  of  gyration  for  the  end-grafted  chain. 
The  bars  above  and  below  the  points  indicate  the  statistical  errors  obtained  in  the 
calculation. 
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Figure  4.19:  Force  profiles  when  compressing  a 60-segment  chain  with  an  aligned  and 
misaligned  fiat  bottom  monolith  {Lt  = 60).  In  the  former  case,  the  monolith,  with 
its  center  at  (x,  y)  = (100, 100),  and  the  grafted  end  of  the  polymer  chain  are  aligned 
one  over  the  other.  In  the  latter,  the  center  of  the  tip  is  located  at  (x,  y)  = (90, 90). 

the  chain  escapes  and  the  occurrence  of  an  ‘inflection’  point  in  the  force  profile  depend 
on  the  actual  size  and  shape  of  the  tip  and  the  alignment  of  the  tip  with  respect  to 
the  polymer  chain.  The  transition  region  can  be  either  small  or  not  clearly  defined 
depending  on  the  above  considerations  (see,  for  example,  the  case  of  a tip  with  a 
hemispherical  bottom  in  Figure  4.20).  These  observations  imply  that  it  may  be  very 
difficult  to  detect  this  transition  experimentally.  However,  as  shown  in  the  above 
figures,  we  expect  the  fluctuations  in  the  force  to  increase  markedly  within  this  region, 
which  suggests  that  recent  methods  based  on  the  analysis  of  the  noise  of  the  AFM 
tip  [81,  118,  119,  120]  might  serve  as  better  probes  to  detect  the  transition  region 
from  the  force  profile. 

Figure  4.21  illustrates  the  influence  of  the  temperature  on  the  force  profile. 
The  results  were  obtained  for  the  same  system  shown  in  Fig.  4.17,  but  at  a higher 
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Figure  4.20:  This  figure  summarizes  the  results  obtained  for  Sets  IV,  V and  VI  in 
Table  4.3.  The  broken  lines  correspond  to  the  statistical  errors  in  the  determination 
of  the  force. 

temperature,  i.e.,  at  good-solvent  conditions.  The  functional  form  of  the  force  for 
a confined  chain  in  good  solvents  differs  noticeably  from  what  is  obtained  for  theta 
solvents.  (Scaling  arguments  for  good  solvents  suggest  that  / ~ [108].)  It  is 

also  evident  that,  for  the  particular  case  shown  in  Figure  4.21,  an  increase  in  the  tem- 
perature makes  the  transition  region  smaller  and  almost  imperceptible.  This  suggests 
that  there  is  minimum  length  below  which  thermal  fluctuations  blur  the  character- 
ization of  the  states  as  confined  or  escaped.  More  recent  off-lattice  simulations  [19] 
have  indeed  confirmed  that,  in  good  solvents,  the  transition  is  observed  only  for  large 
chains  {N  > 500). 

Although  presented  in  the  context  of  the  interaction  between  a polymer  mush- 
room and  an  AFM  tip,  the  phenomenon  described  here  has  broader  implications,  both 
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Figure  4.21;  Force  profile  for  a 100-segment  chain  and  a flat  bottom  monolith  (Lj  = 
120)  in  a good  and  6 solvents. 

to  other  force  measurement  techniques  and  particle/substrate  interactions  in  the  pres- 
ence of  polymers.  For  example,  instead  of  an  AFM,  one  can  use  evanescent  wave  light 
scattering  from  a particle  held  by  an  optical  trap  to  measure  particle/surface  force 
as  a function  of  H [81].  This  technique  is  much  more  sensitive  than  the  AFM  and 
is  capable  of  detecting  forces  of  the  order  of  picoNewtons.  The  phenomenon  de- 
scribed in  this  section  is  of  importance  in  interpreting  such  measurements,  as  it  is  for 
other  applications  such  as  colloid  stability  and  biological  cellular  interactions.  Some 
additional  details  on  these  issues  may  be  found  in  Refs.  [15]  and  [18]. 

4.3.2  Compression  of  a Polymer  Brush 

In  this  section  we  consider  an  extension  of  the  above  single  chain  case  to  that 
of  a multi-chain  system,  namely,  the  force  of  compression  of  a polymer  brush  in  a 
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Table  4.4:  Compression  of  a polymer  brush:  Simulation  parameters. 

Case  Tip  Size  Tip  Shape 

I Lt  = 20  Monolith,  Flat  Bottom 

II  = 10  Cylinder,  Hemispherical  Bottom 

III  i?t  = 5 Cylinder,  Hemispherical  Bottom 

Note:  The  grafting  density  is  held  constant  at  0.1.  The  number  of  chains,  m,  is  40, 

each  with  40  segments. 

theta  solvent.  The  lattice  simulations  used  for  this  purpose  also  employ  the  three- 
dimensional  bond-fluctuation  model  (see  Appendix  A).  A lattice  size  of  40  x 40  x 100 
is  used  with  impenetrable  walls  at  z = 0 and  2:  = 100  and  with  periodic  boundary 
conditions  in  the  x and  y directions.  The  simulations  consider  m chains  with  N 
segments  each,  and  each  chain  is  end-grafted  to  the  lower  wall  in  a periodic  array. 
The  shape  and  the  size  of  the  tip  used  to  obtain  the  results  are  listed  in  Table  4.4. 

To  the  best  of  our  knowledge,  only  one  previous  attempt  has  been  made  on  the 
simulation  of  the  interaction  between  an  AFM  tip  and  a polymer  brush.  This  work, 
by  Murat  and  Crest  [121],  uses  molecular  dynamics  simulations  to  obtain  the  effect  of 
the  shape  and  size  of  the  tip  on  the  force  profile  under  good-solvent  conditions.  The 
results  show  that  the  polymer  chains  splay-out  upon  compression,  thereby  giving  rise 
to  a deviation  of  the  force  from  theoretical  predictions  based  on  the  assumption  that 
the  chains  cannot  escape  from  under  the  tip  in  the  brush  regime  [109]. 

In  contrast  to  the  work  of  Murat  and  Crest  [100,  121],  the  present  work 
focuses  on  the  AFM  tip/polymer  brush  interactions  under  theta  conditions.  We  have 
selected  the  theta  condition  here  in  view  of  the  interesting  conformational  transition 
shown  in  Figure  4.17  for  single  chains.  It  is,  therefore,  interesting  to  examine  if  such  a 
transition  exists  in  the  case  of  a polymer  brush.  (The  calculations  of  Murat  and  Crest 
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[121],  which  are  based  on  simulations  of  100-segment  polymer  chains,  do  not  indicate 
such  a transition.  However,  our  results  for  good  solvent  conditions  reported  in  the 
previous  section  suggest  that  the  demarcation  of  the  states  as  confined  or  escaped 
and,  therefore,  the  ‘transition  region’  are  only  observed  for  very  large  chains.) 

We  first  present,  in  Figure  4.22,  the  interaction  between  a cylindrical  tip  with 
a hemispherical  bottom  and  a polymer  brush.  We  have  chosen  to  report  (in  Figure 
4.22)  the  ratio  of  the  force  exerted  on  the  tip  and  the  radius  of  the  tip,  Le.,  f /Rt- 
According  to  theoretical  calculations  based  on  the  Derjaguin  approximation  [2],  this 
ratio  is  related  to  the  Helmholtz  energy  of  interaction  between  two  flat  plates.  This 
ratio  is,  therefore,  supposed  to  depend  only  on  the  characteristics  of  the  polymer 
layer  and  not  on  the  characteristics  of  the  probe.  Figure  4.22  illustrates  however, 
the  qualitatively  different  behavior  that  one  can  obtain  when  probing  a layer  with 
finite-sized  objects. 

The  differences  in  Figure  4.22  are  caused  by  the  splay-out  or  ‘escaping’  of  the 
chains  from  under  the  AFM  tip.  In  Figure  4.23  we  replot  the  force  of  compression  due 
to  a cylindrical  tip  (this  time  we  present  only  the  force  exerted  on  the  AFM  tip)  along 
with  the  result  for  compression  by  a larger  object,  a monolith  with  a flat  bottom. 
We  have  also  defined  = (Rt/Rgxy),  the  relative  size  of  the  tip  Rt  with  respect 
to  the  radius  of  gyration  R^xy  of  the  grafted  chains  in  the  x and  y directions.  (The 
magnitude  of  R^xy  is  related  to  the  separation  between  the  grafting  points  on  the 
surface.  For  the  examples  shown  here,  R^xy  is  roughly  5.6  lattice  units.).  Figure  4.23 
illustrates  that  for  < 1,  the  induced  forces  are  very  weak  since  the  tip  penetrates 
the  brush  easily  (and,  as  a result,  the  force  is  essentially  constant),  whereas  for  > 1 
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Figure  4.22:  Force  between  a cylindrical  AFM  tip  with  a hemispherical  bottom  and 
a polymer  brush  in  a theta  solvent.  The  simulation  parameters  are  specified  in  Table 
4.4.  The  plot  illustrates  that  the  interaction  energy  measured  with  finite-sized  objects 
depends  not  only  on  the  characteristics  of  the  layer,  but  also  on  the  characteristics  of 
the  probe. 

conformational  rearrangements  in  the  chains  can  be  seen  for  the  conditions  illustrated 
in  the  figure.  These  rearrangements  manifest  themselves  in  the  force  profile  through  a 
‘transition’  region.  Details  of  this  transition  region  are,  however,  not  easily  accessible 
from  lattice  simulations.  Off  lattice  simulations,  which  permit  finer  intervals  in  H, 
are  required  for  this  particular  task. 

The  conformational  rearrangements  can  also  be  observed  from  the  segment 
density  profiles  shown  in  Figure  4.24  for  three  different  positions  of  the  AFM  tip  (the 
monolith,  in  this  case).  These  positions  are  marked  by  Points  A,  B and  C in  Figure 
4.23  and  correspond  to  i/t  = 50,20  and  12,  respectively.  Curve  A in  Figure  4.24 
illustrates  the  density  profile  for  an  unperturbed  brush  (f.e.,  the  tip  is  quite  far  from 
the  brush).  The  segment  density  profile  corresponding  to  Position  B (Curve  B in 
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Figure  4.23;  Force  between  an  AFM  tip  and  a polymer  brush  in  a theta  solvent.  The 
simulation  parameters  are  specified  in  Table  4.4.  The  radius  of  gyration,  Rgxy,  in  the 
x-y  plane  under  these  conditions  is  5.6  lattice  units. 

Figure  4.24)  shows  that  the  chains  under  the  tip  are  in  a ‘confined’  state,  indicated 
by  the  increased  density  for  H ~ 16.  (The  decrease  in  density  for  16  ~ < i/t  = 20 

is  the  result  of  packing  restrictions.)  For  H > 20,  the  density  is  essentially  the  same 
as  that  for  a free  brush.  The  situation  in  Curve  B thus  corresponds  essentially  to 
pure  compression  of  the  chains  beneath  the  tip.  In  contrast,  a significant  level  of 
‘escape’  can  be  seen  for  Curve  C (/ft  = 12),  as  evident  from  the  increase  in  the 
density  for  H > 12.  We  should  point  that  the  densities  reported  in  Figure  4.24  are 
‘local’  densities,  computed  over  a small  area  in  the  vicinity  of  the  tip.  For  instance, 
in  the  particular  example  shown  here,  we  use  an  area  four  times  the  cross-sectional 
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Figure  4.24:  Segment  density  profiles  of  the  brush  for  three  different  positions  of  the 
AFM  tip  as  indicated  in  Figure  20  for  the  monolith.  The  parameters  used  correspond 
to  Case  I in  Table  4.4.  Densities  are  in  number  of  segments  per  unit  area  parallel  to 
the  surface. 

area  of  the  tip  for  calculations  beneath  the  tip  {H  < Ht)  and  three  times  the  tip  area 
for  H > Ht  {i.e.,  the  area  occupied  by  the  tip  is  not  available  for  the  chains).^ 

Understanding  the  relation  between  the  conformational  rearrangements  of 
polymers  and  the  induced  forces  is  one  of  the  key-issues  needed  for  engineering  poly- 
mer layers  for  different  applications.  (For  instance,  understanding  the  interactions 
between  finite-sized  particles  and  polymer  layers  has  important  implications  for  the 
dispersion  of  nanoparticles,  drug  delivery,  and  (biological)  cellular  interactions.)  Our 


^ The  density  profiles  shown  in  Figure  4.24  are  meant  to  merely  illustrate  the 
conformational  changes  arising  from  the  compression.  In  general,  for  detailed  quan- 
titative analyses,  the  area  used  for  calculating  the  densities  must  be  small  enough  to 
reflect  the  density  changes  in  the  neighborhood  of  the  tip  accurately.  In  addition,  the 
size  of  the  simulation  box  must  be  large  enough  to  avoid  the  effect  of  the  periodic 
boundary  condition  on  the  calculations. 
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goal  in  this  section  was  to  illustrate  how  CDM  can  be  used  for  such  a purpose,  i.e., 
to  provide  a direct  link  between  the  structure  and  forces,  and  to  provide  guidelines 
to  complement  both  theories  and  experiments. 

Finally,  we  should  point  out  that  in  the  examples  presented  above,  the  AFM 
tip  is  not  fully  contained  within  the  simulation  box  but  ‘grows’  from  the  upper  surface. 
Furthermore,  the  tip  does  not  have  any  physisorbed  or  grafted  chains.  The  contact- 
distribution  method  can  be  applied  to  situations  involving  objects  that  are  fully 
contained  in  the  simulation  box  or  objects  that  have  some  affinity  for  the  polymer 
chains.  The  estimation  of  the  changes  in  free  energy  in  these  cases  require  some 
additional  steps.  In  the  next  section  we  briefly  describe  how  to  implement  CDM  to 
other,  more  complex  problems  of  relevance  to  polymer  and  colloid  physics. 

4.4  CDM  as  a Tool  for  Examining  Polymer-Induced  Forces  in  Colloidal  Systems 
In  this  section  we  provide  some  guidelines  on  how  to  implement  CDM  in 
problems  that  involve  additional  complexity,  such  as  (i)  the  interaction  between  two 
surfaces  bearing  end-grafted  polymers,  (n)  the  interaction  with  finite-sized  objects, 
fully  contained  within  the  simulation  box,  and  {in)  the  interaction  between  two  sur- 
faces bearing  physisorbed  polymers.  These  three  cases  are  depicted  in  Figure  4.25.^ 

^ Our  goal  here  is  not  to  review  or  specify  the  issues  that  should  be  examined 
in  the  type  of  problems  shown  in  Figure  4.25,  but  to  illustrate  how  CDM  can  be 
implemented  in  such  systems. 
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Figure  4.25:  Examples  of  polymer/colloid  interactions  of  practical  relevance. 

4.4.1  Interaction  between  Two  Surfaces  Bearing  End-Grafted  Polymers 

Here  we  shall  look  at  the  interaction  between  two  surfaces  bearing  end-grafted 
polymers,  such  as  the  interaction  between  two  polymer  brushes  or  the  interaction 
between  a polymer  brush  and  a polymer  mushroom.  The  application  of  CDM  to 
examine  this  class  of  problems  requires  some  further  modifications  on  how  Nc,  the 
number  of  contacts,  is  defined.  In  all  the  work  we  have  presented  so  far,  Nc  was 
defined  as  the  number  of  segments  occupying  the  layer  next  to  the  moving  boundary 
of  the  object.  In  the  present  case,  however,  contacts  occur  not  only  between  the 
segments  and  the  object  but  also  between  the  polymer  segments  of  the  two  objects^ 


® A contact  between  two  segments  is  defined  as  follows.  Consider  two  surfaces,  say 
T’  and  ‘2’,  with  Surface  1 at  z = 0 and  the  other  at  z = H.  Let  the  coordinates  of 
a segment  of  a chain  attached  to  Surface  1 be  {x,  y,  z}.  This  segment  is  said  to  be 
in  contact  with  a segment  of  a chain  from  Surface  2 if  the  coordinates  of  the  second 
segment  are  {x,y,z  -\-\} 
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° Contacts 
° ‘Entanglements’ 


Figure  4.26;  Interaction  between  end-grafted  polymers.  Definition  of  contact  and 
entanglement. 

(see,  for  example,  Figure  4.26).  Therefore,  the  repulsive  potential  that  is  applied  in 
CDM  at  the  moving  boundary  is  not  applied  in  the  present  case  at  a fixed  layer  next 
to  the  object  but  at  any  site  where  a contact  occurs. 

Moreover,  the  use  of  CDM  alone  is  not  sufficient  for  obtaining  the  change  in 
the  Helmholtz  potential  in  the  present  case,  as  there  are  a number  of  conformations 
d = H — 1 which  cannot  be  generated  just  from  the  compression  of  a conformation 
aX  d = H.  An  example  of  such  type  of  conformations  are  those  with  entanglements 
between  the  chains  belonging  to  the  two  different  surfaces®  , as  shown  in  Figure  4.26. 
A correction  to  the  Helmholtz  potential  determined  from  either  CDM  or  a thermody- 
namic approach,  say  A.Fcompression)  is  needed  and  is  obtained  by  running  simulations 


® An  entanglement  is  said  to  occur  whenever  a segment  of  a chain  from  Surface  2 is 
in  contact  with  a segment  (at  position  {x,  y,  z})  from  Surface  1,  but  with  coordinates 
{x,y,z-l}. 
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at  d = H — 1 and  getting  the  probability  Pg  of  the  absence  of  entanglements.  The 
net  change  in  the  Helmholtz  potential  is  then  given  as 

AP  = AP jompression  T InPg,  (4.23) 

Additional  details  on  this  type  of  calculations  have  been  presented  by  Dickman  [105], 
who  examined  the  interaction  between  two  polymer  brushes  under  good  solvent  con- 
ditions. 

4.4.2  Interactions  with  Finite-Sized  Objects,  Fully  Contained  within  the  Simulation 

Box 

In  the  simulations  presented  earlier  in  this  chapter,  the  ‘hard  body’  interacting 
with  the  polymer  chains  was  either  a semi-infinite  object  (such  as  a bare  surface)  or 
a finite-sized  object  that  ‘grew’  from  one  of  the  walls  (such  as  the  AFM  tip).  There 
are,  however,  a large  class  of  problems  in  which  the  object  is  fully  contained  within 
the  simulation  box.  The  calculation  of  the  force  for  this  case  requires  some  additional 
steps,  as  described  below. 

Consider  for  instance,  the  case  of  a colloidal  or  nanoparticle  immersed  in  a 
polymer  solution,  and  let  us  assume  that  we  are  interested  in  the  (polymer-mediated) 
interaction  between  the  particle  and  a wall.  One  can  obtain  the  interaction  force  from 
the  change  in  the  Helmholtz  potential  as  the  distance  between  the  particle  and  the 
surface  is  changed  from  D lattice  units  (denoted  here  as  System  a)  to  D — 1 lattice 
units  (denoted  as  System  ^),  as  shown  in  Figure  4.27. 

In  order  to  obtain  the  difference  in  Helmholtz  potential  between  Systems  a 
and  ,5,  we  first  apply  a repulsive  potential  V at  the  moving  boundary  of  the  particle. 
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Figure  4.27:  Polymer  mediated  interaction  between  a colloidal  particle  and  a wall. 

However,  in  this  case  oo  corresponds  to  a distortion,  rather  than  a displacement, 
of  the  particle,  as  shown  in  Figure  4.27.  The  use  of  the  procedure  described  in 
Section  4.1  leads  to  the  change  in  free  energy  between  System  a and  the  distorted 
case  (System  7).  In  order  to  obtain  the  change  in  free  energy  between  Systems  a and 
/?,  an  additional  step  is  needed,  namely,  the  free  energy  change  between  Systems  7 
and  j3.  One  can  determine  the  free  energy  change  in  going  from  System  ^ to  System 
7 by  applying  the  potential  V on  the  other  side  of  the  particle.  The  required  result, 
is  then  obtained  from  the  difference  between  the  two  results  AJ-a-y  and 

i.e., 

ATap  = AJ^ay  ~ ATp^y  , (4.24) 


as  shown  in  the  figure. 


106 


4.4.3  Interactions  between  Two  Physisorbed  Layers 

As  illustrated  in  the  previous  example,  the  actual  path  used  to  calculate 
Helmholtz  energy  differences  is  immaterial  and  can,  in  fact,  be  ‘unphysical’.  In  order 
to  examine  the  interaction  between  two  surfaces  bearing  physisorbed  polymers  we 
also  need  to  follow  a somewhat  circuitous  path.  This  path,  which  is  illustrated  in 
Figure  4.28,  consists  of  three  steps:  [i]  desorption  of  the  chains  from  one  of  the  walls 
(taken  to  be  the  upper  one  here),  (ii)  moving  that  wall  to  the  new  position,  and  (in) 
re-adsorption  of  the  chains  back  on  that  wall.  To  achieve  this,  the  calculation  requires 
the  use  of  a sequence  of  fictitious  potentials  at  the  upper  surface,  Uext{z  = H).  During 
step  (i),  Uext{H)  is  gradually  increased  (in  increments  of  Se  > 0)  starting  from  — e, 
until  the  probability  of  zero  surface  contacts  is  obtained  accurately.  In  particular,  this 
proceeds  as  follows.  One  changes  the  external  potential  in  steps,  i.e.,  Uext{H)  = — e, 
Uext{H)  = — e+(5e,  Uext{H)  = —e+2Se, ...,  and  determines  the  distribution  of  segments 
in  the  surface  layer  (and,  therefore,  the  probability  of  contacts)  at  each  stage.  The 
number  of  such  increments  needed  depends  on  the  conditions  of  the  simulation  (such 
as  the  number  of  polymer  chains,  the  separation  between  the  plates  and  the  original 
magnitude  of  the  adsorption  energy). 

Once  the  probability  of  having  zero  contacts  is  determined  with  sufficient  ac- 
curacy, the  energy  change  corresponding  to  the  compression  is  obtained  as  described 
in  Section  4.2.  The  energy  change  determined  in  this  manner  characterizes  the  differ- 
ence in  Helmholtz  potential  that  occurs  when  the  adsorptive  wall  moves  by  one  lattice 
unit  and  simultaneously  becomes  non-adsorptive.  Step  (ii)  therefore  corresponds  to 
a polymer  system  of  width  i?  — 1,  in  equilibrium  with  one  adsorptive  wall  (the  lower 
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Figure  4.28:  Interaction  between  two  physisorbed  layers.  Illustration  of  the  desorp- 
tion/readsorption path. 

one)  and  one  non-adsorptive  wall  (the  upper  one).  Finally,  in  step  (Hi)  the  surface 
potential  at  the  upper  wall,  Uext{z  = H — 1)  is  gradually  decreased  from  0 to  — e,  so 
that  the  chains  are  re-adsorbed  at  the  upper  wall.  (This  calculation  is  illustrated  in 
Section  4.2.)  The  net  effect  of  these  three  steps^  is  to  go  from  two  adsorptive  walls 
at  a distance  H to  two  adsorptive  walls  separated  by  a distance  H — 1.  The  sum  of 
the  free  energy  changes  obtained  in  these  three  steps  is,  therefore,  the  desired  total 
change  in  the  Helmholtz  energy. 

In  the  next  two  chapters  we  will  use  the  methodology  described  here  to 
examine  the  bridging  force  induced  by  a single  polymer  chain  confined  between  two 

^ It  is  acually  more  convenient  to  proceed  in  the  reverse  order,  from  step  (Hi)  to  {i), 
when  calculating  the  interaction  between  two  physisorbed  layers.  When  proceeding 
in  the  reverse  order,  one  can  use  the  final  configuration  from  step  {Hi)  as  the  initial 
configuration  for  (i)  thus  reducing  the  time  required  to  generate  a new  system. 
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adsorptive  walls  and  the  interaction  force  between  two  undersaturated  polymer  layers. 
The  results  obtained  from  such  analysis  addresses  some  of  the  gaps  and  key  issues  we 
have  identified  in  Chapter  2,  and  provide  important  guidelines  to  better  understand 
the  forces  induced  by  physisorbed  polymers. 


CHAPTER  5 

ADSORBED  SINGLE  CHAINS:  PRELUDE  TO  STARVED  LAYERS 


In  this  chapter^  we  focus  on  the  conformations  and  forces  induced  by  single 
polymer  chains  confined  between  adsorptive  surfaces.  The  correct  description  of 
single  chains  at  interfaces  enhances  one’s  understanding  of  more  realistic  many-chain 
problems  of  direct  technological  relevance  [122].  Also  important  in  this  context  are 
interactions  with  adsorbing  walls  or  with  undersaturated  (starved)  layers,  where  the 
so-called  bridging  interactions  play  a notable  role.  A bridge  is  a section  of  a polymer 
chain  that  is  attached  (either  physically  adsorbed,  or  chemically  grafted)  to  two 
different  bodies.  Understanding  the  mechanism  for  bridging  is  crucial  for  controlling 
the  quality  of  adhesives.  Bridge  formation  is  also  important  in  colloid  stabilization, 
or  rather  destabilization  by  flocculation,  and  crystallization  of  molten  polymer,  where 
bridging  between  crystalline  seeds  affects  the  morphology  of  the  crystalline  structure. 

Models  of  bridging  in  a variety  of  multi-chain  systems  have  used  mean-field 
[12,  22,  26,  31,  32,  33,  38,  39,  41,  47,  123],  scaling  [9,  12,  22,  26,  47,  45],  and  Monte 
Carlo  [97,  124,  125]  methods.  Evidence  for  the  attractive  forces  caused  by  bridging 
has  been  obtained  by  measuring  directly  the  forces  between  two  surfaces  using  the 
surface  force  apparatus  (SEA)  [126,  127]  and  atomic  force  microscopy  (AFM)  [112], 
and  indirectly  from  neutron  diffraction  experiments  on  suspensions  [128,  129].  In  the 

^ This  chapter  is  a slightly  modified  version  of  Jimenez  et  al.  [21]. 
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measurements  reported  recently  by  Swenson  et  al.  [128,  129],  the  effect  of  adding  poly 
(ethylene  oxide)  (PEO)  to  an  aqueous/ clay  suspension  is  analyzed.  The  magnitude  of 
the  polymer-related  forces  is  inferred  from  accurately  measured  decreases  in  colloidal 
spacing  after  addition  of  PEO.  The  flat  geometry  of  the  clay  particles  makes  it  an 
ideal  system  for  a bridging  analysis.  At  each  separation  between  the  clay  particles, 
Swenson  et  al.  estimate  the  number  of  bridges  per  unit  area  and  suggest  that  a 
constant  force  per  bridge  of  roughly  1.4  pN  develops  in  the  system  under  equilibrium 
conditions. 

There  are  fewer  theoretical  [22,  130,  131]  and  experimental  [84,  132,  133] 
studies  on  the  simpler  case  of  the  bridging  interactions  of  an  isolated  chain.  Single- 
chain bridging  may  be  considered  insignificant  under  conditions  typically  realized 
in  technological  applications.  However,  even  in  such  cases  single-chain  events  could 
play  a role  in  the  early  stages  of  a process.  For  instance,  measurements  of  the  total 
adsorbed  amount  [7]  of  polymers  from  dilute  solutions  are  high  enough  to  indicate  that 
the  surface  is  still  characterized  by  a high  degree  of  chain  interactions,  although  the 
polymer  chains  in  the  bulk  solution  do  not  overlap.  Nevertheless,  it  is  possible  to  have 
isolated  chains  at  the  surface  for  an  appreciable  amount  of  time  before  equilibrium  is 
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reached  in  many  systems^  . Biological  systems  also  exhibit  more  frequently  bridging 
interactions  due  to  single  biopolymer  chains. 

DiMarzio  and  Rubin  [130]  were  the  first  to  study  bridge  formation,  using 
an  isolated  ideal  chain.  Their  matrix  formalism,  which  was  extended  and  used  to 
treat  multi-chain  systems  by  Scheutjens  and  Fleer  {e.g.,  see  Ref.  [41]),  provided  early 
numerical  results  for  the  distribution  of  bridges  and  force  as  a function  of  separation 
and  adsorption  energy.  More  relevant  for  good  solvents  is  the  scaling  analysis  of  Ji 
et  al.  [22],  are  obtained  for  the  average  size  and  number  of  bridges.  A combination 
of  these  results  with  the  assumption  that  each  bridge  is  individually  governed  by  the 
Pincus  law  of  elasticity  [70]  yields  a force  per  bridge  that  is  independent  of  separation 
between  the  surfaces  - much  like  the  neutron  diffraction  results  mentioned  above.  The 
scaling  analysis  of  Ji  et  al.  also  predicts  that  the  force  per  bridge,  ft^,  scales  with 
adsorption  energy  as 

(5.1) 

where  eo  is  the  effective  theoretical  adsorption  energy  and  a is  the  size  of  a segment 
or  a repeating  unit. 


^ It  is  believed  that  polymer  physisorption  is  characterized  by  extremely  long  ad- 
sorption times  (of  the  order  of  hours  or  days)  as  shown,  for  example,  in  the  experi- 
ments of  Luckham  and  Klein  on  PEO  [126].  A few  hours  after  adding  polymer  to  a 
good  solvent,  an  overall  attractive  interaction  between  surfaces  is  observed.  A steady 
state  is  found  many  hours  later,  in  which  an  overall  repulsion  is  consistently  observed. 
There  is  little  knowledge  beyond  this  on  the  kinetics  of  physisorption,  especially  on 
the  amount  and  structure  of  polymer  chains  at  the  surface  during  the  incubation 
phase. 
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Illustrations  (meant  as  caricatures)  of  a single  adsorbed  chain  are  shown  in 
Figure  5.1a  &:  b.  The  first  case  (Figure  5.1a)  corresponds  to  a system  with  strongly 
adsorbing  surfaces  such  that  there  are  relatively  few  bridges  and  these  are,  in  general, 
highly  stretched.  Most  of  the  chain  is  tightly  adsorbed  to  the  surface  in  short  loop 
conformations^  . This  is  usually  referred  to  as  the  strong-stretching  bridge  regime. 
Strong  stretching  is  easier  to  model  theoretically  since  the  role  of  loop  conformations 
may  safely  be  ignored,  and  only  the  elasticity  of  non-interacting  bridges  needs  to  be 
accounted  for,  as  suggested  by  Ji  et  al.  [22]. 

A weak-stretching  regime  (Figure  5.1b)  occurs  when  the  surfaces  are  less 
adsorptive.  This  is  characterized  by  a large  number  of  slightly  stretched  bridges  and 
fewer,  but  larger,  loop  conformations.  Finally  in  the  limit  of  large  separation  and 
adsorption  energy,  the  chain  adsorbs  only  onto  one  wall. 

A full  understanding  of  the  bridging  force  and  interactions  requires  additional 
details  such  as  the  number  and  size  distribution  of  bridges  and  possible  interaction 
of  the  bridges  with  loops  and  tails  of  the  chain.  To  our  knowledge  there  have  been 
no  computer  simulations  that  attempt  to  study  directly  the  force  induced  by  a single 
polymer  chain  that  adsorbs  onto  two  parallel  surfaces.  Here,  we  use  CDM  to  examine 
the  force  induced  by  bridges.  One  of  our  objectives  is  to  test  the  predictions  of 
the  scaling  analysis  presented  by  Ji  et  al.  [22].  In  addition,  we  examine  in  detail 

^ Tail  conformations  (unadsorbed  ends  of  the  chain)  become  irrelevant  as  the  chain 
length  is  sufficiently  large.  The  numbers  of  bridges  and  loops  increase  with  chain 
length  while  the  number  of  tails  is  always  two.  Furthermore,  the  tail-size  cannot 
grow  indefinitely  since  tails  will  form  bridges  once  their  size  is  comparable  to  the  gap 
width. 
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Figure  5.1:  Single  chain  confined  between  two  adsorbing  walls,  (a)  Large  e,  (b)  small 
€. 

whether  there  is  any  effect  due  to  interactions  between  bridges  of  a single  chain,  the 
distribution  of  bridge  sizes,  and  possible  interference  from  loop  conformations.  We 
also  provide,  in  Section  5.3,  a scaling  analysis  for  a related  problem,  the  retraction  of 
a strongly  adsorbing  single  chain  in  a good  solvent. 

5.1  Details  of  the  Simulation  Procedure 
We  model  the  polymer  chain  as  a confined  self-avoiding  walk  (SAW)  inscribed 
on  a simple  cubic  lattice  (see  also  Appendix  A).  Steps  into  either  surface  layer  are 
favored  by  a Boltzmann  factor  exp{—e/kBT).  The  number  of  segments  in  the  chain, 
N,  is  chosen  as  1000,  which  corresponds  to  a very  long  polymer  chain,  for  which 
scaling  laws  should  hold.  Based  on  our  experience  from  neutral-wall  simulations,  this 
chain  length  proved  to  be  sufficient  to  test  the  scaling  limit  [122].  If  the  direction,  z, 
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is  normal  to  both  surfaces,  the  discrete  surface-segment  potential  is 


Uext  — ^ 


oo  z < 1,  z > H 
—e  z — z = H 
0 I <z  < H 


(5.2) 


The  method  used  for  sampling  the  conformations  of  the  polymer  chain  is 
based  on  the  ‘configurational  bias’  algorithm  of  Siepmann  and  Frenkel  [134]  which, 
as  mentioned  in  Appendix  A,  consists  of  breaking  an  existing  chain  configuration 
into  two  parts,  erasing  one  of  these  subchains,  and  randomly  regrowing  it  using  an 
algorithm  based  on  the  well-known  Rosenbluth-Rosenbluth  method  [135].  The  force 
exerted  by  the  polymer  chain,  /,  is  estimated  from  the  change  in  Helmholtz  potential, 
T,  that  occurs  when  the  separation  between  the  walls  is  decreased  from  H to  H — 1, 


i.e., 


m - 1/2)  = 


dh 


H-l/2 


T{H  - 1)  - 


(5.3) 


defined  so  that  a negative  force  implies  an  attraction  between  the  walls.  Here,  h is 
used  to  represent  the  wall-separation  as  a free  variable. 

The  Helmholtz  potential  difference  in  Equation  5.3  is  obtained  using  the 
recently  developed  contact-distribution  method  (CDM)  [15].  As  described  in  Chapter 
4,  the  path  used  to  calculate  the  free  energy  difference  is  immaterial  and  can,  in  fact, 
be  ‘unphysicaP.  The  path  required  in  CDM  when  the  moving  wall  is  adsorptive 
consists  of  three  steps:  (^)  desorption  of  the  chain,  {ii)  moving  the  surface  to  the 
new  position,  and  {in)  re-adsorption  of  the  chain.  Specific  details  of  the  simulation 


procedure  are  presented  in  Chapter  4,  Section  4.4.3. 
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A snapshot  of  the  bridges  in  a typical  case  is  shown  in  Figures  5.2a  & b.  The 
forces  arising  from  bridging  have  two  competing  factors.  Attractive  forces  arise  from 
the  tension  in  bridge-structures,  and  repulsive  forces  arise  from  excluded-volume  in- 
teractions and  the  entropic  cost  of  confining  the  chains.  In  Table  5.1  the  forces  are 
listed  for  adsorption  energies  ranging  from  0.3  to  1.0  ksT  for  separations  between 
the  surfaces,  H/a  (with  a being  the  size  of  the  lattice  unit),  ranging  from  4 up  to  the 
separation  at  which  bridging  ceases  (between  7 and  11).  Separations  below  approx- 
imately 4 lattice  units  exhibit  a deviation  from  universal  (be.,  lattice  independent) 
behavior  in  some^  but  not  all  of  the  system  properties. 

5.2  Polymer  Conformations  and  Induced  Forces 
In  what  follows,  we  first  focus  on  the  number  of  bridges  as  a function  of 
separation  between  the  walls  and  the  force  due  to  the  bridges.  The  latter  will  also 
be  compared  with  results  based  on  scaling  arguments.  Following  this  we  shall  ex- 
amine the  distribution  of  bridges  in  terms  of  size  and  under  what  conditions  steric 
interactions  between  the  bridges  and  the  rest  of  the  chain  and  the  effects  of  exclusion 
of  conformations  by  the  wall  set  in.  Finally,  in  Section  5.3,  we  present  a scaling 
argument  for  the  force  of  retraction  of  an  adsorbed  chain  and  compare  it  with  the 
bridging  forces  observed  in  the  simulation. 

^ The  concentration  profile  will  be  certainly  affected,  but  properties  such  as  the 
lateral  component  of  the  radius  of  gyration  would  not. 
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Figure  5.2:  Snapshots  from  a simulation  showing  a side  view  (2-D  projection)  of 
a chain  and  a 3-D  perspective  of  a section  of  a chain.  The  black  segments  denote 
adsorbed  segments  or  trains;  the  gray  segments  are  the  loops,  and  the  light-gray  ones 
belong  to  bridges.  Case  (a)  corresponds  to  H = 5 and  e = 1.0  ksT,  and  case  (b)  to 
H = 6 and  e = 0.5  fcsT. 
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Table  5.1:  Force  induced  by  a single  polymer  chain  at  different  adsorption  energies. 


H/a 

0.3 

0.5 

e/kjT 

0.6 

0.75 

0.85 

1.0 

4.5 

-4.65 

-16.69 

-21.24 

-23.84 

-23.65 

-22.19 

5.5 

-0.55 

-11.08 

-12.21 

-12.63 

-10.84 

-7.79 

6.5 

-0.21 

-6.80 

-7.58 

-6.22 

-4.86 

-3.36 

7.5 

-0.37 

-4.60 

-5.30 

-3.15 

-2.54 

-2.03 

8.5 

-0.66 

-3.28 

-2.48 

-1.13 

— 

— 

9.5 

-0.12 

-2.82 

-0.95 

— 

— 

— 

10.5 

-0.58 

— 

— 

— 

— 

— 

5.2.1  Number  of  Bridges 

In  Figure  5.3,  the  number  of  bridges,  Ubr,  is  plotted  as  a function  of  H/a 
for  a sequence  of  adsorption  energies.  The  trend  seen  is  consistent  with  the  numer- 
ical results  of  DiMarzio  and  Rubin  [130],  which  predict  that  the  number  of  bridges 
vanishes  when  the  adsorption  energy  becomes  overwhelming  and  the  chain  adopts 
completely  flattened  conformations.  At  each  separation,  a maximum  in  the  number 
of  bridges  must  occur  at  a certain  e that  is  somewhere  between  strongly  adsorptive 
and  repulsive.  Figure  5.3  suggests  that  this  maximum  occurs  at  a magnitude  of  e 
below  0.5  ksT  for  all  of  the  separations  shown  since  none  of  the  curves  overlap.  The 
decrease  of  the  number  of  bridges  with  separation,  shown  in  Figure  5.3,  is  consistent 
with  the  exponential  decay  of  the  scaling  prediction  [22] : 

/£)\  -5/3 

~ 77  ( — j exp{-H/2D),  (5.4) 

where  H is  a length  comparable  to  the  average  loop  extension  and  depends  only  on 


the  adsorption  energy. 
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Figure  5.3;  Number  of  bridges,  Ubr,  as  a function  of  the  separation  between  surfaces 
H/a,  where  a is  the  lattice  spacing.  Exponential  trendlines  are  shown  for  e > Q.bksT. 

5.2.2  Bridging  Forces 

For  all  e > 0.3  fcflT,  the  magnitude  of  the  observed  force  (see  Figure  5.4) 
decreases  monotonically  with  increasing  separation  between  surfaces.  For  the  lowest 
energy,  e = 0.3  the  force  fluctuates  near  zero  for  large  separations  and  becomes 
slightly  attractive  for  small  separations.  This  indicates  that  the  system  is  close  to 
the  ‘threshold’  adsorption  energy,  deflned  as  the  point  at  which  the  volume  fraction 
profile  becomes  uniform  (the  adsorption  energy  at  which  the  energetic  profit  is  exactly 
balanced  by  the  entropic  cost  for  the  chain  to  assume  an  adsorbed  conformation^  ). 
In  what  follows  we  examine  adsorption  energies  in  the  range  between  0.5  and  1.0 
keT. 


^ This  observation  is  actually  confirmed  by  the  analysis  of  the  segment- 
concentration  profile  between  the  walls. 
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Figure  5.4:  Force,  /,  as  a function  of  the  separation  between  surfaces  H/a  , for 
different  values  of  the  adsorption  energy. 

As  we  have  noted  earlier,  one  might  idealize  the  system  under  study  as  two 
planes  connected  by  several  elastic  tethers  individually  governed  by  the  Pincus  law 
(described  in  Chapter  3,  Section  3.1).  This  model  ignores  possible  repulsive  contri- 
butions arising  from  loop  conformations,  and  the  interaction  between  bridges  that 
might  alter  their  collective  elcisticity.  With  this  in  mind,  it  is  interesting  to  study 
the  force  per  bridge,  i.e.,  the  ratio  of  the  total  force  and  the  number  of  bridges 
obtained  from  the  simulation: 


1/2)  = 


f{H  - 1/2) 

nbr{H  - 1/2)' 


(5.5) 


Note  that  this  calculation  requires  interpolation  of  the  number  of  bridges  since  Ubr  is 
obtained  at  integer- valued  separations:  H and  H — 1.  The  uncertainty  introduced  in 
the  values  of  fbr  by  the  interpolation  of  ntr  is  not  expected  to  be  of  major  concern 
since  the  ni^  data  are  relatively  smooth. 
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Figure  5.5:  Force  per  bridge,  fbr,  vs.  H/a  for  a range  of  adsorption  energy. 

The  force  per  bridge,  fbr,  is  plotted  as  a function  of  H in  Figure  5.5,  and 
as  a function  of  e in  Figure  5.6.  These  results  yield  direct  evidence  that  there  is  a 
constant  force  per  bridge.  The  observed  relation  between  the  force  per  bridge  and 
the  adsorption  energy  is  in  very  good  agreement  with  Eqn.  5.1.  This  result  insinuates 
that  the  interaction  between  different  bridges  does  not  play  a major  role  and  that  the 
forces  are  simply  additive.  A more  careful  examination  of  this  statement  is  presented 
below,  where  we  analyze  the  distribution  of  the  bridge  sizes  for  a number  of  adsorption 
energies. 

5.2.3  Distribution  of  Bridges 

The  interpretation  of  the  results  presented  in  the  previous  section  is  facilitated 
by  the  size  distribution  of  the  bridges  obtained  for  various  adsorption  energies.  Figure 
5.7  shows  the  probability  that  a bridge  has  a certain  degree  of  stretching,  H/Rbr- 
Here  Rbr  is  the  root-mean-squared  end-to-end  length  of  an  unperturbed  chain  with 
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Figure  5.6:  Force  per  bridge  vs.  adsorption  energy  e.  The  solid  line  corresponds  to 
fbr  = 1.42ei-48. 

as  many  segments  as  are  in  each  bridge  (known  as  the  natural  length  of  the  bridge). 
This  length  is  characterized  by  the  well  known  Flory  exponent,  Ri,r  ~ where 

the  prefactor  depends  on  the  model  [97].  In  other  words,  what  is  plotted  is  the 
probability  that  a bridge  with  s segments  is  stretched  by  a factor  of  approximately 
H/as^!^  times  its  natural  length.  It  is  evident,  from  Eqn.  5.1,  that  for  e = 1.0  ksT 
most  of  the  bridges  are  stretched  beyond  their  natural  length,  whereas  for  e = 0.5 
HbT  there  is  a considerable  amount  (~  20%)  of  large,  loosely  stretched  bridges. 


122 


Figure  5.7:  Discrete  probability,  p{x),  that  a bridge  is  stretched  by  a factor  x,  where 
X = H/Rhr  ~ s being  the  number  of  segments  in  a bridge  (see  text).  Sep- 

aration is  at  H/a  = 5.  For  c = 0.5  ksT,  80%  of  the  bridges  are  strongly  stretched 
{H/Rbr  > !)•  For  e = 1.0  ksT,  96%  of  the  bridges  are  strongly  stretched. 

The  tension  in  each  strongly  stretched  bridge  can  be  described  by  a law 
proposed  by  Pincus  [70],  which  in  our  context,  is® 

^3/2 

/ ~ — 5/2  , H > Rbr.  (5.6) 

The  restoring  force  of  those  bridges  that  are  not  strongly  stretched  is  described  by 
the  spring-like  linear  relation  between  force  and  separation  [6], 

H<Rhr.  (5.7) 

^br 


® The  exact  relation  we  use  for  the  Flory  radius  is  Rbr/ a = 0.99A^®  ®^®  which  is 
suitable  for  chain  lengths  above  20  on  a simple  cubic  lattice. 
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The  strong-  and  weak-stretching  regimes  described  by  Eqns.  5.6  and  5.7  have 
been  observed  in  lattice  Monte  Carlo  simulations  on  single  chains  stretched  by  a 


assume  that  the  chain  can  take  any  conformation  which  ends  exactly  H layers  from 
its  starting  point,  and  hence  do  not  account  for  the  presence  of  the  walls.  It  appears 
that  Eqn.  5.7  will  not  be  entirely  valid  as  a description  of  weakly  stretched  bridges 
since,  for  this  particular  case,  there  would  be  a significant  number  of  disallowed, 
wall-penetrating  conformations. 

The  force  profiles  corresponding  to  the  simulations  in  Figure  5.7  are  shown 
in  Figures  5.8  and  5.9.  The  dashed  lines  represent  a theoretical  comparison  based 
on  the  distribution  of  bridges  found  in  the  simulation.  If  the  polymer  chain  has  n{s) 
bridges  with  s segments  each,  the  total  force  can  be  estimated  as 


with  the  force  due  to  each  bridge  being  calculated  using  either  Eqn.  5.6  or  Eqn.  5.7, 
as  appropriate. 

As  is  evident  from  Figure  5.9,  the  assumption  that  a system  is  governed 
completely  by  the  elasticity  of  its  bridges  does  not  hold  in  cases  where  there  is  a sig- 
nificant number  of  weakly  stretched  bridges.  In  such  cases,  other  effects  such  as  steric 
forces  (bridge-bridge  interactions  and  interaction  of  bridges  with  the  rest  of  the  chain) 

^ The  upper  limit  for  the  validity  of  this  equation  has  been  established  through 
simulations  [78]  to  be  approximately  2Rbr.  When  the  chain  is  stretched  to  twice  its 
natural  size,  the  force  becomes  model-dependent. 


constant  force  applied  at  both  ends^  (see  Chapter  3,  Section  3.1).  These  relations 


(5.8) 
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Figure  5.8:  Attractive  force  between  walls,  vs.  separation  for  e = 1.0  ksT.  The 
scaling  prediction  is  based  on  the  use  of  the  appropriate  law  of  elasticity  (i.e.,  the 
Pincus  law  or  the  Hookean  relation). 

and  exclusion  of  conformations  that  penetrate  the  walls  diminish  the  attraction.  In 
multi-chain  systems  the  presence  of  steric  forces  can  go  as  far  as  to  create  an  overall 
repulsion  between  the  surfaces.  The  images  shown  in  Figure  5.2  demonstrate  the 
steric  effects  that  exist  in  the  cases  of  weak  and  strong  adsorption.  The  prediction  of 
force  per  bridge  that  scales  with  adsorption  energy  as  shown  in  Eqn.  5.1  is  probably 
limited  to  conditions  where  the  adsorption  strength  is  high  enough  that  almost  every 
bridge  is  strongly  stretched. 


5.3  Scaling  Relation  for  Pulling  a Chain  off  of  a Surface 
Testing  Eqn.  5.1  using  standard  experiments  is  difficult  because  of  the  diffi- 
culty of  isolating  the  contributions  from  strongly  stretched  bridges  from  the  overall 
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Figure  5.9:  Attractive  force  between  walls,  -/,  vs.  separation  for  e = 0.5  ksT.  The 
broken  line  represents  the  estimate  based  on  the  scaling  laws  {i.e.,  the  Pincus  or  the 
Hookean  relation). 

measurement.  Therefore,  it  is  useful  to  conceive  of  methods  to  test  Eqn.  5.1  more  di- 
rectly. Pulling  a chain  off  of  a surface  using  an  atomic  force  microscope  may  provide 
the  type  of  information  necessary  to  examine  Eqn.  5.1. 

A theoretical  analysis  of  an  adsorbed  (Gaussian)  chain  subjected  to  an  ex- 
ternal force  has  shown  that  the  competition  between  stretching  and  adsorption  leads 
to  a rich  behavior  with  first-  and  second-order  transitions  [136,  137,  138].  A simple 
scaling  analysis  for  this  problem  was  suggested  recently  for  a weakly  adsorbing,  ideal 
chain  [139].  In  what  follows  we  adapt  this  approach  to  obtain  a scaling  relation  be- 
tween the  force  and  the  adsorption  energy  for  the  case  of  strong  stretching  in  a good 
solvent.  First  one  determines  the  spatial  extension  (“thickness”)  D of  an  adsorbed 
chain  following  the  arguments  presented  by  de  Gennes  [6]  by  assuming  the  following 
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free  energy  functional: 

where  D = D/a  and  9 — e/ksT.  The  first  term  describes  the  free  energy  of  confine- 
ment for  an  adsorbed  layer  of  N segments  and  is  based  on  the  argument  that  the 
chain  is  a sequence  of  blobs  of  extension,  D,  each  with  an  energy  ksT.  The  second 
term  is  an  estimate  of  energy  due  to  the  surface  contacts  assuming  the  chain  segments 
are  evenly  distributed  in  the  thickness  D.  One  can  now  minimize  the  free  energy  with 
respect  to  D to  find  the  equilibrium  layer  thickness  (see  de  Gennes  [6],  Chapter  1), 

D*  ~ (510) 

For  a chain  being  pulled  by  one  of  its  segments,  we  modify  the  free  energy 
functional  in  Eqn.  5.9  as  follows  when  there  is  a bridge  with  p segments: 

\ 5/2 

where  H = H /a.  Eqn.  5.11  is  arrived  at  by  substituting  N —pfov  N and  D*  for  D in 
Eqn.  5.9  and  by  adding  the  second  term  on  the  right-hand  side  to  describe  the  bridge 
between  the  surface  and  the  AFM  tip  by  the  elastic  free  energy  expression  which 
corresponds  to  Eqn.  5.6.  The  minus  sign  on  the  first  term  in  Eqn.  5.11  is  necessary 
in  order  that  the  free  energy  of  the  adsorbed  part  of  the  chain  decreases,  rather  than 
increases,  with  increasing  chain  length,  temperature  and  adsorption  strength.  The 
minimization  of  Eqn.  5.11  with  respect  to  the  number  of  segments  in  the  ‘tether’ 
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yields  the  equilibrium  size  of  the  tether, 


p*  ~ 


6 ■ 


(5.12) 


The  Helmholtz  free  energy  J^{p)  at  p = 

ksT 


p*  is  then  used  to  find  the  equilibrium  force, 

^3/2 

~ . (5.13) 


It  is  not  surprising  that,  as  in  the  case  of  bridging  between  walls,  the  force  is  inde- 
pendent of  separation  and  depends  on  adsorption  energy  through  an  exponent  of  1.5 
since  the  same  type  of  balance  between  bridging  and  adsorbing  segments  exists.  In 
addition  to  establishing  a link  between  similar  problems,  this  correspondence  confirms 
that  the  physical  content  of  the  two  scaling  approaches  (namely,  the  one  presented 
here  and  the  one  by  Ji  et  al  [22])  is  the  same,  even  though  the  approach  by  Ji  et  al. 
[22]  is  based  on  an  extension  of  mean-field  calculations. 

This  scaling  analysis  suggests  that  an  experiment  may  be  designed  to  obtain 
direct  information  on  the  adsorption  energies  of  surfaces  from  the  measured  bridging 
forces.  Therefore,  it  is  useful  to  consider  how  a direct  comparison  of  the  simula- 
tion results  may  be  made  with  such  experimental  data.  The  correspondence  between 
real  chains  and  universal  polymer  models  used  in  simulations  depends  mainly  on 
the  relative  stiffness  of  the  type  of  polymer  being  considered  (see  also  Appendix  A). 
For  instance,  for  PEO  in  water  (a  good  solvent)  the  lattice  scale  corresponds  approxi- 
mately to  0.6  nm,  with  each  segment  in  the  lattice  containing  about  1.6  ethylene-oxide 
monomers.  Then,  our  data  for  the  force  and  adsorption  strength  of  (/,  e)  = (1.4 
ksT/a,  1.0  keT)  would  correspond  to  (/,  e)  = (9  pN,0.4  kcal/mol).  It  should  be 
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borne  in  mind  that  this  estimate  can  be  only  as  good  as  the  accuracy  of  the  cor- 
respondence between  coarse-grained  models  and  real  chains.  Therefore,  the  above 
comparison  should  be  viewed  mainly  as  an  instructive  example  of  the  possibility  for 
exploring  surface/polymer  characteristics. 

We  should  caution,  however,  that  Eqn.  5.13  was  developed  under  the  as- 
sumption that  the  system  remains  at  equilibrium,  and  hence  the  length  of  the  tether 
between  the  surface  and  the  AFM  tip  is  allowed  to  adjust  itself  at  each  separation. 
In  contrast,  this  is  not  usually  the  case  in  the  experiments  reported  in  the  litera- 
ture. That  is,  if  the  velocity  of  retraction  of  the  AFM  tip  is  large  compared  to  the 
relaxation  time  of  the  adsorbed  chain,  the  measured  force  would  be  determined  by 
the  appropriate  force-length  relationship  (which  follows  the  Pincus  law  for  moderate 
forces  and  becomes  model-dependent  for  large  stretching),  and  Eqn.  5.13  would  not 
be  valid  anymore. 


5.4  Toward  Multiple  Chains 

It  is  useful  at  this  stage  to  comment  on  the  case  of  multiple  chains  and 
contrast  it  with  single  chains.  A detailed  examination  of  multi-chain  systems  is 
presented  in  the  next  chapter;  here  we  present  preliminary  results  to  highlight  some 
of  the  similarities  and  differences  the  multi-chain  problem  has  with  the  present  study. 
In  particular,  we  present  results  for  the  forces  observed  between  two  partially  coated 
surfaces  in  a multi-chain  system.  The  system  consists  of  200-segment  chains  in  an 
athermal  solvent  between  two  adsorbing  walls  each  with  e = 1.0  ksT.  A total  of 
12  chains  is  considered  in  a 40  x 40  x lattice,  corresponding  to  3 monomers  for 
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Figure  5.10:  Attractive  force  per  unit  area,  -fjA,  vs.  separation  for  e = 1.0  ksT 
in  a multiple  chain  system  of  twelve  200-segment  chains  corresponding  to  a slightly 
undersaturated  layer.  The  broken  line  corresponds  to  an  estimate  based  on  the  Pincus 
and  Hookean  laws  of  elasticity. 

every  4 surface  sites.  Since  the  equivalent  adsorbed  layer  in  equilibrium  with  a dilute 
solution  has  roughly  1 monomer  for  every  surface  site,  this  simulation  corresponds  to 
a starved  surface.  The  results  obtained  for  a range  of  separation  distances  are  shown 
in  Figure  5.10.  Also  shown  in  the  figure  is  the  theoretical  curve  from  Eqn.  5.8  based 
on  summing  the  Pincus-law  prediction  for  single  bridges  of  a given  size  by  the  number 
of  bridges  of  that  size. 

First,  it  can  be  seen  that  the  net  force  observed  from  the  simulations  is 
attractive.  However,  the  force  observed  falls  well  below  the  prediction  based  on  Eqn. 
5.8,  indicating  that  steric  effects  play  a significant  role  in  multi-chain  systems.  This 
is  especially  clear  from  a comparison  with  Figure  5.8,  which  shows  that  the  Pincus 
law  works  well  for  an  identical,  but  single-chain,  system.  Dividing  the  total  force  by 
the  number  of  bridges  yields,  as  in  the  single-chain  cases,  a roughly  constant  force 
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Figure  5.11:  Force  per  bridge,  f^r,  vs.  H/a  for  a single  chain  (triangles)  and  a multi- 
chain system  (squares).  The  adsorption  energy  for  both  cases  is  e = 1.0  ksT. 

(see  Figure  5.11).  The  ‘constancy’  of  the  observed  force  per  bridge  might  mislead 
one  to  conclude  that  the  force  is  determined  largely  by  non-interacting  bridges  (as 
may  have  occurred  in  the  experimental  observations  by  Swenson  et  al.  [128,  129]). 
However,  it  should  be  borne  in  mind  that  the  observed  force  per  bridge  cannot  be 
related  unequivocally  to  the  adsorption  energy  of  the  surfaces,  nor  is  it  an  indication 
that  one  has  a collection  of  non-interacting  bridges.  The  fact  that  the  data  for 
multiple  chains  at  e = 1.0  ksT  are  well  below  the  corresponding  single-chain  case 
demonstrates  the  importance  of  steric  interactions  such  as  bridge-bridge  interactions, 
bridge-loop  interactions,  etc.  A detailed  analysis  of  the  competition  between  bridging 
attraction  and  steric  repulsion  is  presented  in  the  next  chapter,  where  we  examine 
the  interaction  between  two  undersaturated  polymer  layers. 
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CHAPTER  6 

STARVED  LAYERS:  TESTING  THE  NUMERICAL  MEAN-FIELD 

CALCULATIONS 

In  this  chapter^  we  extend  our  previous  work  to  examine  the  interactions 
between  two  physisorbed  layers  with  surface  coverages  F ranging  from  0 to  the  satu- 
ration coverage  Fq  (defined  near  the  end  of  the  next  section).  One  of  our  objectives 
is  to  complement  the  study  of  bridging  due  to  isolated  polymer  chains  presented  in 
Chapter  5 [21].  In  addition,  we  wish  to  examine  the  competition  between  bridging 
attraction  and  steric  repulsion  that  occur  between  undersaturated  polymer  layers.  As 
noted  in  Chapter  2,  the  mean-field  and  the  scaling  theories  differ  in  their  predictions 
of  the  forces  between  starved  layers.  While  mean-field  theories  indicate  that  attrac- 
tion between  the  surfaces  occurs  for  small  undersaturations  {i.e.,  for  a coverage  2% 
lower  than  the  saturation  value),  scaling/mean-field  hybrids  [9,  11,  12]  indicate  that, 
for  infinitely  long  chains,  the  onset  of  attraction  occurs  for  a rather  large  undersat- 
uration (i.e.,  for  a coverage  60%  lower  than  the  saturation  value).  Here  we  present 
unambiguous  results  for  the  force  of  interaction  between  two  physisorbed  polymer 
layers.  We  also  provide  a detailed  comparison  between  the  simulation  results  and 

^ This  chapter  is  a slightly  modified  and  updated  version  of  Jimenez  et  al.  [23]. 
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numerical  lattice  mean-field  predictions  based  on  first-order^  (standard  Scheutjens- 
Fleer  (SF)  theory)  and  second-order  Markov  chains  (hereinafter  referred  to  as  SFi 
and  SF2,  respectively).  A careful  analysis  of  the  conditions  for  which  the  numerical 
predictions  are  valid  and,  additionally,  an  understanding  of  the  deviations  of  these 
predictions  (in  those  situations  where  mean-field  is  not  appropriate)  are  of  utmost 
importance. 


6.1  Model  and  Simulation  Method 

The  computational  method  used  here  consists  of  three  steps:  (i)  generation 
of  polymer  configurations  using  Monte  Carlo  simulations,  (ii)  determination  of  force 
of  interaction  (i.e.,  change  in  the  Helmholtz  potential  as  the  two  polymer-bearing 
surfaces  approach  each  other),  and  (in)  calculation  of  the  same  quantities  using  self- 
consistent  mean-field  theory.  In  all  cases,  the  space  is  discretized  using  a simple 
cubic  lattice,  and  the  polymer  chains  are  modelled  as  appropriate  lattice  random 
walks.  Some  additional  details  follow  below. 

6.1.1  Lattice  Model 

The  model  used  here  is  identical  to  the  one  used  for  the  analysis  of  a single 
chain  between  two  adsorptive  walls  (see  Chapter  5).  We  model  the  polymer  chains 
as  self-avoiding  walks  (SAW’s)  inscribed  on  a simple  cubic  lattice  with  dimensions 
Lx  Lx  H.  Periodic  boundary  conditions  are  used  in  the  x and  y directions,  whereas 

^ A first-order  Markov  chain  admits  backfolding  of  segments,  although  such  back- 
folding  is  physically  forbidden. 
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two  impenetrable  walls  confine  the  layers  in  the  2:  direction.  The  two  walls  are 
adsorptive,  with  the  discrete  surface/segment  potential  given  as 


00  z < 1,  z > H 


Uext  — 


— e z = 1,  z = H ■ 


0 Kz<H 

K 


(6.1) 


The  chain  length  N used  in  most  of  the  simulations^  we  present  here  is  200. 

For  the  numerical  mean-field  calculations  shown  in  Section  6.3,  we  keep  the 
same  underlying  lattice  structure  {i.e.,  simple  cubic)  and  the  same  adsorption  energy  e 
and  chain  length  N used  in  the  simulations.  In  this  case,  as  noted  earlier,  the  polymers 
are  considered  as  first-  or  second-order  Markov  chains  (in  contrast  to  simulations, 
where  the  chains  are  essentially  {N  — l)^*-order  Markov  chains),  and  the  excluded- 
volume  interactions  are  taken  into  account  using  a mean-field  assumption. 


6.1.2  Simulation  Method  for  Structure  and  Forces 

For  generating  equilibrium  conformations  of  the  polymer  chains,  we  have 
employed  the  configurational  bias  method  [134].  As  described  in  Appendix  A,  in  this 
method,  sampling  of  the  SAW’s  consists  of  breaking  an  existing  chain  configuration 
into  two  parts,  erasing  one  of  these  subchains,  and  randomly  regrowing  it  using  an 
algorithm  based  on  the  well-known  Rosenbluth-Rosenbluth  method. 


^ This  value  of  N is  not  sufficient  to  test  the  scaling  limit  [14],  and  hence,  we  will 
not  be  able  to  make  definite  comparisons  with  the  scaling  free-energy  functional  used 
by  Rossi  and  Pincus  [11,  12]. 
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The  force  exerted  by  the  polymer  chains,  /,  is  calculated  using  the  contact- 
distribution  method  (CDM).  As  noted  in  Chapter  4,  Section  4.4.3,  the  calculation  of 
the  force  when  the  moving  wall  is  adsorptive  follows  a circuitous  path,  which  consists 
of  three  steps:  (i)  desorption  of  the  chains  from  one  of  the  walls  (taken  to  be  the 
upper  one  here),  (n)  moving  that  wall  to  the  new  position,  and  {in)  re-adsorption  of 
the  chains  back  on  that  wall. 

Although,  from  a fundamental  point  of  view,  the  desorption/re-adsorption 
approach  is  valid,  the  propagation  of  errors  in  such  a calculation  may  raise  questions 
concerning  the  accuracy  of  this  methodology.  Therefore,  in  order  to  test  the  accuracy 
of  the  proposed  method,  we  have  calculated  the  force  of  compression  between  a poly- 
mer layer  and  an  athermal  (neutral)  wall  following  two  different  paths.  In  the  first 
and  more  direct  path,  we  ‘move’  the  bare  surface.  In  this  case,  we  only  need  to  deter- 
mine the  probability  of  zero  surface  contacts'^  Pq.  In  the  second  and  more  tortuous 
path,  we  move  the  surface  with  the  adsorbed  chains  (referred  to  as  the  desorption/re- 
adsorption  path).  The  results  of  such  a calculation  are  shown  in  Figure  6.1.  It  is 
clear  from  the  figure  that  the  agreement  between  the  two  approaches  is  very  good 
and  that  the  methodology  used  to  estimate  forces  between  two  physisorbed  layers  is 
indeed  valid.  The  size  of  the  lattice  (L)  used  for  all  the  simulations  presented  in  this 
chapter  is  L/a  — 40,  where  a is  the  lattice  spacing.  The  number  of  polymer  chains 

^ In  most  of  the  cases  one  additional  simulation  with  a repulsive  potential 
Uext{H)  = 0.20  ksT  is  enough  to  get  a good  estimate  of  Pq.  However,  when  the 
chains  are  strongly  confined,  additional  simulations  with  stronger  repulsive  poten- 
tials may  be  required.  See,  for  example.  Chapter  4,  Section  4.2. 
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Figure  6.1:  Compression  of  a polymer  layer  by  an  ‘athermal’  (neutral)  wall,  illus- 
trating the  equivalence  of  two  different  ways  of  determining  the  force  of  confinement. 
The  surface  coverage  on  the  adsorptive  surface  (with  adsorption  energy  e = 1.0  ksT) 
is  r = 1.0  monolayers. 

was  varied  from  n = 1 to  n = 20  in  order  to  obtain  different  coverages  on  the  sur- 
faces. We  define  the  surface  coverage  T as  the  total  number  of  segments  belonging  to 
adsorbed  chains  divided  by  the  area  of  the  adsorptive  surfaces  [i.e.,  it  represents  the 
number  of  monolayers  adsorbed  on  each  surface).  For  instance,  the  results  presented 
in  Figure  6.1  correspond  to  a surface  coverage  F = 1.0  monolayer.  The  number  of 
polymer  chains  used  in  that  particular  simulation^  was  n = 8. 

It  is  also  useful  to  think  in  terms  of  the  surface  coverage  F relative  to  the 
equilibrium  coverage  F*  that  one  would  obtain  when  the  surface  is  in  contact  with  a 
polymer  solution.  Since  F*  is  a function  of  the  bulk  concentration  of  the  solution,  it 


® We  have  performed  simulations  for  the  same  surface  coverage  but  with  a larger 
simulation  box  {L  = 80  and  n — 32)  and  obtained  essentially  the  same  results  [107]. 
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is  common  practice  to  define  a saturation  coverage  Fq  as  the  coverage  obtained  when 
4>b  — > 0.  This  definition  presents  no  problem  when  the  chain  length  N approaches 
infinity,  since  the  adsorption  isotherm  {i.e.,  a plot  of  F*  as  a function  of  (j)b)  behaves 
essentially  like  a step  function  for  N oo.  For  short  chains,  however,  one  obtains 
Fo  from  the  extrapolation  of  F*  (beyond  the  initial  rise)  to  (j)b  = 0.  For  instance,  for 
the  system  under  study  {i.e.,  N = 200  and  e = 1.0  ksT),  the  saturation  coverage 
corresponds  to  1.2  monolayers,  as  determined  by  de  Joannis  and  Bitsanis  [14]. 

6.2  Interaction  between  Two  Undersaturated  Surfaces 
In  the  first  part  of  this  section  we  focus  on  the  number  of  bridges  and  the 
average  size  of  bridges  as  functions  of  the  surface  coverage  F and  the  separation 
between  walls  H/a.  The  chain  length  of  the  polymers  used  in  this  study  is  = 200, 
and  the  adsorption  energy  is  e = 1.0  ksT.  This  section  is  intended  to  complement  the 
work  presented  in  Chapter  5 on  bridging  by  a single  polymer  chain.  In  particular,  we 
want  to  determine  the  surface  coverage  at  which  the  steric  repulsion  sets  in.  Following 
this,  we  examine  the  competition  between  steric  repulsion  and  bridging  attraction 
between  two  undersaturated  surfaces  and  compare  our  results  with  numerical  self- 
consistent  calculations  using  first-order  Markov  chains  (standard  SF  theory,  denoted 
henceforth  as  SFi)  and  second-order  Markov  chains  (SF2),  the  latter  to  eliminate 
chain  backfolding. 


6.2.1  Structure  of  Bridges 
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Bridging  due  to  single  chains 

When  the  surface  coverage  is  very  small,  bridging  between  two  surfaces  may 
occur  due  to  isolated  polymer  chains.  This  limiting  case  has  been  examined  by  Ji  et 
al.  [22],  who  presented  analytical  results  based  on  scaling  arguments  for  the  average 
size  (i.e.,  average  number  of  segments  contained  in  the  bridge)  l^r  and  the  average 
number  of  bridges  ni^,  i.e., 

nfr,  ~ iV  f - j exp{-H/2D),  (6.2) 


and 


I 


br 


H 
D' 


(6.3) 


where  N is  the  chain  length,  a is  the  monomer  size,  and  D is  a length  comparable  to 
the  average  loop  extension,  which  depends  only  on  the  adsorption  energy.  In  order  to 
estimate  the  forces  due  to  bridging,  Ji  et  al.  assumed  each  bridge  to  be  individually 
governed  by  the  Pincus  law  of  elasticity  [70].  In  our  context,  the  Pincus  law  is 


f = C 


{H  - 1)3/2 


p5/2 

^br 


(6.4) 


where,  for  our  lattice  model,  H — 1 is  the  stretched  end-to-end  distance  of  the  bridge, 
Rbr  is  the  unperturbed  end-to-end  distance  of  a chain  containing  as  many  segments 
as  there  are  in  each  bridge,  and  C is  a proportionality  constant®  . This  assumption 
yields  a force  per  bridge  that  is  independent  of  the  separation  distance  H between 


® We  use  C « 3.2,  as  has  been  found  from  simulations  on  chains  with  hard-core 
excluded  volume  [76,  79]. 
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the  surfaces, 

ftxr  ~ D,  (6.5) 

or, 

U ~ (-)  . (6.6) 

since  D as  shown  by  de  Gennes  [6]. 

The  case  of  a semi-dilute  solution  confined  between  two  adsorbing  walls  has 
also  been  considered  by  Ji  et  al.  [22],  who,  in  order  to  make  predictions  for  the  force 
due  to  bridging,  have  assumed  that  each  bridge  is  still  governed  by  the  Pincus  law. 
The  discussion  they  present  in  this  context  is,  however,  speculative.  In  this  case,  one 
would  expect  the  steric  interactions  to  play  a major  role  and,  hence,  it  is  unlikely 
that  the  system  will  behave  as  a collection  of  elastic  tethers. 

Bridging  in  starved  layers 

In  Chapter  5 we  examined  the  case  of  bridging  by  a single  polymer  chain.  We 
presented  computer  simulation  results  for  the  number  and  size  distribution  of  bridge 
conformations  and  the  resulting  attractive  force  between  the  surfaces  as  a function  of 
inter-surface  distance  H/a  and  adsorption  energy  e.  These  results  showed  that  in  the 
case  of  strong  adsorption  (c  > 0.5  ksT),  the  force  per  bridge  f^r  is  indeed  independent 
of  separation  and  depends  only  on  the  adsorption  energy,  in  good  agreement  with  the 
scaling  predictions.  In  this  section  we  extend  our  previous  work  to  examine  the 
structure  of  the  bridges  and  the  forces  induced  by  undersaturated  polymer  layers. 
The  adsorption  energy  used  in  all  our  simulations  is  e = 1.0  HbT. 
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In  Figure  6.2  we  plot  the  average  number  of  bridges  as  a function  of  H/a 
for  different  values  of  the  surface  coverage  F (given  in  number  of  monolayers).  The 
broken  line  in  this  figure  represents  the  single-chain  results  presented  in  Chapter  5 
(obtained  for  a chain  of  1000  segments  and  for  the  same  adsorption  energy,  e = 1.0 
ksT).  As  evident  from  the  figure,  for  low  values  of  T (up  to  F = 0.5),  the  average 
number  of  bridges  shows  the  same  exponential  decay  with  H.  As  the  coverage  is 
increased  beyond  F = 0.5,  we  observe  a slower  (almost)  exponential  decay  {i.e.,  a 
larger  decay  length).  As  pointed  out  by  Ji  et  al.,  the  blob  size  (which  determines 
the  rate  of  exponential  decay  in  Eqn.  6.2)  cannot  be  determined  unambiguously  in 
systems  with  two  similar  characteristic  lengths.  When  the  chains  remain  isolated  the 
only  characteristic  length  in  the  system  is  the  layer  thickness  D (which  depends  only 
on  the  adsorption  energy).  As  the  surface  coverage  is  increased,  spatial  fluctuations 
in  the  system  (which  are  characterized  by  the  mesh  size  ^)  become  important  and 
cause  a deviation  from  Eqn.  6.2.  Ji  et  al.  [22]  speculate  that  when  ^ < D,  the  decay 
length  is  given  by  ^ (instead  of  D),  which  would  decrease  with  increase  in  surface 
coverage.  Figure  6.2  appears  to  be  inconsistent  with  the  suggestion  of  Ji  et  al,  but 
no  definitive  conclusions  may  be  drawn  without  data  for  larger  F’s. 

A similar  situation  can  be  observed  from  Figure  6.3,  which  shows  the  average 
size  of  the  bridges  as  a function  of  H/a.  The  results  show  a linear  relation  between  Ibr 
and  H/a.  For  low  coverages,  this  relation  is  characterized  only  by  the  layer  thickness 
D.  Similar  to  Figure  6.2,  the  deviations  in  the  slope  are  evident  for  F > 0.5. 

A more  interesting  situation  arises  when  one  compares  the  bridging  force 
directly  obtained  from  the  simulations  with  that  estimated  from  the  number  and  size 
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Figure  6.2:  Number  of  bridges  per  unit  area,  Ubr/a^,  as  a function  of  the  separation 
between  surfaces  H/a,  where  a is  the  lattice  spacing. 
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Figure  6.3:  Average  size  of  bridges  Ibr  {i.e.,  average  number  of  segments  within  a 
bridge)  as  a function  of  the  separation  between  surfaces,  H/a. 
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distribution  of  bridges  (i.e.,  following  the  assumption  that  each  bridge  is  individually 
governed  by  the  Pincus  law).  This  comparison  is  presented  in  Figure  6.4,  which  shows 
the  force  per  bridge  fbr  obtained  at  a fixed  inter-surface  distance  {H/a  = 4.5)  as  a 
function  of  the  surface  coverage  F.  The  force  per  bridge  is  defined  as  the  ratio  of  the 
total  force  and  the  number  of  bridges  obtained  from  the  simulation.  The  solid  line  in 
the  figure  corresponds  to  the  estimation  based  on  the  number  and  size  distribution  of 
the  bridges^  . The  deviation  between  the  two  curves  indicates  the  onset  of  the  steric 
repulsion  between  the  layers.  (A  similar  behavior  is  observed  in  Chapter  5 on  bridging 
by  a single  chain  for  low  adsorption  energies,  i.e.,  e < 0.5  ksT.)  The  fact  that  fbr 
obtained  directly  from  simulations  is  approximately  constant  as  F — > 0 indicates  that 
for  very  low  coverages  the  surface  can  be  considered  as  a collection  of  isolated  polymer 
chains.  In  this  regime  fbr  is  determined  by  the  segment-surface  adsorption  energy. 
As  the  coverage  is  increased,  steric  repulsion  sets  in  and  the  calculation  of  the  ratio 
between  the  total  force  and  the  number  of  bridges  loses  its  meaning. 

6.2.2  Competition  between  Bridging  Attraction  and  Steric  Repulsion 

The  force  profiles  obtained  for  a number  of  surface  coverages  are  shown  in 
Figures  6.5  and  6.6.  As  discussed  in  the  previous  section,  at  very  low  coverages  the 
surface  can  be  considered  as  a collection  of  isolated  polymer  chains.  An  increase  in 
the  coverage  gives  rise  to  an  increase  in  the  attractive  force,  as  seen  in  Figure  6.5, 

^ Calculations  based  on  the  average  number  and  average  size  of  the  bridges  are 
roughly  10%  lower  than  the  ones  shown  in  Figure  6.4,  which  are  calculated  using  the 
entire  number  and  size  distribution  of  bridges. 


142 


Figure  6.4:  Force  per  bridge,  f^r,  versus  F for  a fixed  separation  distance  between 
the  surfaces  {H/a  = 4.5).  The  solid  line  (scaling  prediction)  is  obtained  from  the 
appropriate  laws  of  elasticity  (i.e.,  the  Pincus  or  the  Hookean  laws). 

because  of  the  increase  in  the  number  of  bridges  (see,  for  example.  Figure  6.2).  As 
the  surface  coverage  increases  further,  the  polymer  chains  begin  to  overlap,  and  a 
competition  between  bridging  attraction  and  steric  repulsion  sets  in.  In  addition, 
the  average  size  of  the  bridges  also  increases  (see  Figure  6.3)  and  the  attractive  force 
therefore  extends  to  a longer  range  (see  Figure  6.5). 

For  further  increases  in  F,  the  steric  repulsion  becomes  dominant  and  a de- 
crease in  the  attractive  force  is  observed  (see  Figure  6.6).  For  a particular  coverage 
(r  ~ 1.12  in  our  case),  the  net  force  between  the  walls  becomes  roughly  zero  (except 
at  very  close  distances  where  a steep  repulsion  is  always  observed).  It  is  interesting  to 
note  that  the  value  of  F for  which  the  force  essentially  vanishes  is  below  the  saturation 
coverage  Fq,  and  corresponds,  approximately  to  F w 0.85  Fq.  (As  mentioned  earlier, 
Fo  was  obtained  from  the  extrapolation  of  the  adsorption  isotherm  to  4>b  = 0.  For 
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Figure  6.5:  Force  per  unit  area  / versus  H/a  for  a range  of  surface  coverage  F between 
0 and  0.75. 

N = 200,  e = 1.0  ksT,  Fo  = 1.2  monolayers.)  Any  increase  in  the  coverage  beyond 
the  saturation  value  results  in  a steric  repulsion  at  all  distances. 

To  better  illustrate  the  competition  between  bridging  attraction  and  steric 
repulsion  between  the  two  polymer  layers,  we  show,  in  Figure  6.7,  the  net  force 
at  a close  separation  {H/a  = 4.5)  as  a function  of  the  surface  coverage  F.  This 
figure  illustrates  several  of  our  previous  observations,  that  is,  the  linear  increase  of 
the  (attractive)  force  for  very  low  coverages  (where  the  surface  consists  of  isolated 
polymer  chains),  the  balance  between  bridging  and  steric  forces  and  the  onset  of 
purely  repulsive  forces  for  a coverage  of,  approximately,  85  % of  the  saturation  value 
(for  N = 200).  We  also  show  in  this  figure  a few,  additional  simulations  for  chain 
length  N = 100.  More  data  is  needed,  particularly  for  larger  chains,  before  one  can 
attempt  to  draw  any  conclusions  on  the  effect  of  chain  length  on  the  observed  force. 
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Figure  6.6;  Force  per  unit  area  / versus  H/a  for  a range  of  surface  coverage  F between 
0.75  and  1.25. 


Figure  6.7:  Force  per  unit  area  / at  a fixed  separation  {H/a  = 4.5)  as  a function  of 
surface  coverage,  F. 
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At  this  point  it  is  interesting  to  comment  on  the  predictions  based  on  the 
scaling/free-energy  functional  for  undersaturated  layers.  As  pointed  by  Rossi  and 
Pincus  [11,  12],  a coverage  lower  than  Pq  leads,  essentially,  to  a decrease  in  the  repul- 
sion between  the  surfaces.  Only  when  the  relative  surface  saturation  P/Po  becomes 
lower  than  0.4  does  the  net  force  become  attractive  (in  contrast  to  simulations,  where, 
for  N = 200,  the  forces  become  attractive  at  relative  saturations  P/Pq  < 0.85).  We 
note,  however,  that  the  chain  length  used  in  this  study  is  not  large  enough  to  test 
the  asymptotic  predictions  of  the  scaling  theory.  We  comment  more  on  this  issue 
in  the  next  section,  where  we  compare  our  simulation  results  with  lattice,  numerical 
mean-field  calculations  (Scheutjens-Fleer  theory). 

Before  we  proceed  to  the  comparison  with  the  lattice  mean-field  calculations, 
it  is  instructive  to  consider  a computer  ‘experiment’  designed  for  illustrative  purposes. 
In  this  ‘experiment’  we  compare  the  force  of  interaction  obtained  under  two  different 
conditions  (see  Figure  6.8).  In  the  first  case  (denoted  here  as  the  standard  case),  all 
polymer  segments  are  allowed  to  adsorb  on  both  surfaces.  This  situation  leads  to  the 
competition  between  bridging  attraction  and  steric  repulsion  we  have  shown  earlier 
in  this  section.  In  the  second  case  we  assume  that  the  polymer  chains  adsorbed  on 
the  upper  surface  do  not  adsorb  on  the  lower  surface.  (Similarly,  the  polymer  chains 
in  the  lower  surface  do  not  adsorb  on  the  upper  one.)  The  interactions  in  this  case 
are,  therefore,  purely  repulsive. 

The  interaction  force  resulting  from  this  experiment  is  shown  in  Figure  6.9. 
The  difference  between  these  profiles  serves  as  an  illustration  of  the  competition 
between  bridging  attraction  and  steric  repulsion  that  occurs  in  the  standard  case. 
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Figure  6.8:  Interaction  between  two  physisorbed  layers,  (a)  All  polymer  segments 
are  allowed  to  adsorb  on  both  surfaces.  This  case  leads  to  a competition  between 
bridging  attraction  and  steric  repulsion,  (b)  The  polymers  adsorbed  on  the  upper 
surface  do  not  adsorb  on  the  bottom  surface,  and  vice  versa.  This  case  leads  to  purely 
steric  repulsion. 
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Figure  6.9:  Force  per  unit  area  / versus  H/a  for  the  standard  case  (where  both  steric 
and  bridging  contributions  are  present),  and  for  a layer  that  leads  only  to  repulsive 
contribution. 
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One  can  also  argue  that  a situation  similar  to  that  presented  in  Figure  6.8b 
may  arise  in  non-equilibrium  conditions  {i.e.,  in  cases  where  the  characteristic  time  of 
the  compression  is  smaller  than  the  relaxation  time  of  the  polymers  within  the  gap). 
If  the  polymer  chains  cannot  rearrange  and  equilibrate  within  the  characteristic  time 
of  the  interaction,  it  is  possible  to  observe  repulsive  forces  even  if  the  surfaces  are 
undersaturated. 

At  this  point  it  is  important  to  note  that  although  one  can  use  simulations  to 
obtain  details  on  the  structure  and  interactions  under  a wide  variety  of  conditions, 
resorting  to  simulations  is  not  always  convenient  for  routine  guidance  in  practice. 
Therefore,  it  is  more  appropriate  and  useful  to  use  simulations  to  test  and  improve 
simpler,  analytical  and  numerical  approaches  for  predictive  purposes. 

6.3  Interaction  between  Two  Starved  Layers:  Comparison  with  Numerical 

Mean-Field  Theory 

It  is  important  to  compare  our  simulation  results  with  numerical  mean-field 
calculations  (based  on  the  Scheutjens-Fleer  formulation),  as  the  latter  are  the  only 
predictions  available  for  short  polymers.  Such  a comparison  can  be  used  to  identify 
the  most  useful  refinements  needed  for  improving  the  accuracy  of  the  numerical  lat- 
tice mean-field  theory.  The  original  theory  of  Scheutjens  and  Fleer  is  based  on  the 
conformations  of  a first-order  Markov  chain  subject  to  a mean-field  potential.  As 
well-known  [13],  this  formulation  contains  two  critical  assumptions,  namely,  (i)  ran- 
dom mixing  approximation  in  each  layer,  and  (u)  the  possibility  of  chain  backfolding. 
In  this  chapter,  in  addition  to  examining  the  original  formulation,  we  also  relax  the 
second  assumption  by  using  a second-order  Markov  chain  for  the  step-weighted  walk. 
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Figure  6.10:  Average  (attractive)  force  per  bridge,  f^r,  as  a function  of  adsorption 
energy  e.  The  results  shown  in  this  figure  are  for  single  chains. 

(For  additional  details,  see  Appendix  B.)  This  extension  of  the  standard  SF  theory 
was  introduced  by  Scheutjens  and  Leermakers  [140]  and  later  used  by  van  der  Linden 
et  al.  [44]  to  examine  the  adsorption  of  semi-flexible  copolymers.  It  is  interesting  to 
mention  that  for  a long  time  it  was  believed  that  the  use  of  a second-order  Markov 
chain  would  only  modify  the  length  of  the  statistical  segment  of  the  chain  (as  it  does 
for  polymers  in  solution)  [7].  However,  it  is  now  recognized  that  at  least  some  of 
the  structural  details  of  the  layer  (such  as  the  size  distribution  of  trains)  are  indeed 
sensitive  to  the  flexibility  of  the  chains.  A detailed  examination  of  the  layer  structure 
for  a second-order  Markov  chain  is  not  available  currently.  Some  preliminary  results 
are  shown  in  Appendix  B. 
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6.3.1  Bridging  by  Single  Polymer  Chains 

We  first  present  a brief  comparison  between  simulations  and  numerical  lattice 
calculations  for  the  case  of  a single  chain®  confined  between  two  adsorptive  surfaces. 
We  then  follow  with  a detailed  comparison  for  the  interaction  between  two  undersat- 
urated surfaces.  The  case  of  a single  chain  between  two  adsorptive  surfaces  can  also 
be  considered  as  that  of  a layer  with  a very  low  saturation  coverage.  As  discussed 
previously,  for  isolated  chains  (and  strong  adsorption)  one  can  take  the  ratio  of  the 
force  and  the  number  of  bridges  and  define  an  average  force  per  bridge,  fi^.  In  Figure 
6.10  we  show  the  average  (attractive)  force  per  bridge  as  a function  of  the  adsorption 
energy  e.  This  figure  shows  the  results  obtained  from  simulations  (see  Chapter  5, 
Figure  5.6)  and  the  ones  calculated  for  a first-order  (SFi)  and  a second-order  (SFa) 
Markov  chain.  As  evident  from  the  figure,  the  magnitude  of  fi^  increases  as  the  chain 
becomes  more  fiexible.  For  instance,  for  e = 1.0  ksT,  one  obtains  fbr,sFi  = 2.2  ksT/a 
for  a first-order  Markov  chain,  fbr,SF2  = 1-7  ksT/a  for  a second-order  Markov  chain 
and  fbr  ySim  — 1*4  ksT/a  for  a self-avoiding  walk  (from  simulations).  To  understand 
these  results  one  should  note  that  as  the  flexibility  of  a chain  increases,  its  unper- 
turbed size  decreases.  Therefore,  the  increase  of  with  the  flexibility  is  mainly  due 
to  the  corresponding  decrease  in  R^r  (see  Eqn.  6.4). 
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Figure  6.11:  Average  number  of  bridges  per  unit  area  as  a function  of  H/a  for 

r = 0.75. 
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Figure  6.12:  Average  size  of  bridges  /ftr  as  a function  of  H/a  for  F = 0.75. 
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6.3.2  Interactions  between  Two  Undersaturated  Layers 

We  next  consider  the  interaction  between  two  undersaturated  layers.  As  an 
illustrative  example,  we  show,  in  Figures  6.11  & 6.12,  the  number  and  size  of  the 
bridges  as  a function  of  H/a  for  a surface  coverage  F = 0.75.  The  improvement 
introduced  by  the  second-order  Markov  chain  in  the  prediction  of  the  structure  of 
the  bridges  is  evident  (especially  at  close  separations),  both  in  the  number  and  the 
size  of  the  bridges.  More  interesting,  however,  is  the  comparison  of  the  prediction 
for  the  force,  shown  in  Figure  6.13.  The  differences  in  this  case  are  drastic,  with  the 
simulations  giving  rise  to  a much  stronger  attractive  force  than  the  ones  predicted  by 
SFi  or  SF2.  Although  the  numerical  mean-field  calculations  based  on  the  second-order 
Markov  chain  (SF2)  yields  good  predictions  of  the  structural  properties  of  bridges, 
the  underestimation  of  the  saturation  coverage  Fq  affects  the  steric  contribution  to 
the  force. 

In  order  to  better  illustrate  the  differences  in  the  predictions,  we  show  in 
Figure  6.14  the  net  force  at  a close  separation  {H/a  = 4.5)  as  a function  of  the 
surface  coverage  F.  As  mentioned  earlier,  this  figure  illustrates  the  linear  increase 
of  the  (attractive)  force  with  a corresponding  increase  in  coverage.  In  this  regime, 
the  forces  are  mainly  due  to  intra-chain  interactions.  The  comparison  between  the 
predictions  of  SFi  and  SF2  also  confirms  our  previous  observations  that  the  use  of  a 

* The  restriction  of  the  numerical  lattice  calculations  for  single  chains  necessarily 
eliminates  the  use  of  a mean-field  potential.  Therefore,  such  calculations  are  based  on 
the  conformations  obtained  for  a first-order  {i.e.,  an  ideal  chain)  and  a second-order 
Markov  chain,  in  the  absence  of  any  mean  field.. 
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Figure  6.13:  Force  per  unit  area  / versus  separation  H/a  for  F = 0.75. 

second-order  Markov  chain  improves  the  standard  SF  theory  {i.e.,  SFi)  substantially. 
In  particular,  the  linearity  of  the  attractive  force  with  coverage  is  captured  by  SF2, 
although  SF2  still  falls  short  of  the  expected  magnitude  and  the  rate  of  increase, 
as  predicted  by  the  simulations.  As  the  coverage  is  increased  further,  one  observes 
a much  larger  deviation  between  mean-field  calculations  and  computer  simulations. 
The  forces  predicted  by  SF  1 and  SF 2 become  repulsive  at  lower  coverages  as  compared 
to  simulations,  which  is  something  one  would  expect  due  to  the  fact  that  the  mean- 
field  theories  underestimate  the  saturation  coverage.  For  instance,  for  N = 200 
and  e = 1.0  ksT,  the  saturation  coverage  obtained  from  mean-field  calculations  is 
To.sfi  = 0-7  for  the  first  order  Markov  chain  and  Fq, sfz  = 0.72  for  the  second-order 
Markov  chain,  in  contrast  to  Fo,sim  = 1.2  obtained  from  simulations. 

Due  to  the  differences  in  the  saturation  coverages  obtained  from  the  mean- 
field  theory  and  the  simulations,  one  would  expect  a better  agreement  if  the  forces 
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Figure  6.14:  Force  per  unit  area  at  a fixed  separation  between  surfaces  {H/a  = 4.5) 
as  a function  of  surface  coverage,  F. 

are  compared  for  the  same  relative  surface  saturation  F/Fq.  In  order  to  examine  this, 
we  have  replotted  the  data  in  Figure  6.14  in  terms  of  F scaled  with  respect  to  Fq;  see 
Figure  6.15.  Such  a comparison,  in  fact,  reveals  discrepancies  not  evident  in  Figure 
6.14.  For  example,  we  observe  from  our  simulations  that  (for  N = 200)  a crossover 
from  attraction  to  steric  repulsion  occurs  at  F/Fq  0.85.  According  to  the  numerical 
mean-field  calculations,  this  crossover  occurs  at  a value  of  F/Fq  slightly  larger  than 
1.0.  In  general,  one  would  expect  better  agreement  between  the  simulations  and  the 
mean-field  theories  for  larger  coverages  (i.e.,  oversaturated  surfaces),  where  spatial 
fluctuations  in  segment  density  are  expected  to  be  less  pervasive. 

An  observation  concerning  for  what  values  of  F/Fq  the  crossover  from  attrac- 
tion to  repulsion  occurs  is  in  order.  As  N oo,  the  mean-field  calculations  (based 
on  both  the  two-order-parameter  theory  of  Semenov  and  Joanny  [38]  and  the  SF 
theory)  predict  a steric  repulsion  for  F/Fq  = 1.0  for  all  values  of  H.  However,  this 


154 


Figure  6.15:  Force  per  unit  area  at  a fixed  separation  between  surfaces  {H/a  = 4.5) 
as  a function  of  the  relative  surface  coverage,  F/Fq. 

steric  repulsion  turns  into  attraction  even  for  very  low  undersaturations  {i.e.,  even  for 
small  deviations  below  1.0).  For  instance,  the  two-order-parameter  theory  predicts 
attraction  for  any  F/Fq  ~ 0.98,  as  noted  in  Chapter  2.  It  is  interesting  to  contrast 
the  above  mean-field  prediction  with  that  of  the  scaling/free-energy  functional  the- 
ory, which  shows  a sharp  crossover  from  steric  repulsion  to  attraction  at  a very  low 
value  of  r/Fo  0.4).  As  we  pointed  out  earlier,  in  our  simulations  we  observe  the 
crossover  at  F/Fq  « 0.85  for  N = 200.  It  would,  therefore,  be  interesting  to  extend 
the  simulations  to  longer  chain  lengths  to  examine  whether  the  asymptotic  prediction 
of  the  scaling/free  energy  functional  theory  is  reached®  . 


® It  is  not  clear  what  the  limitations  of  the  ‘scaling-based’  free-energy  functional 
proposed  by  de  Gennes[27]  are.  Examination  of  the  extensions  of  this  theory  for  the 
case  of  undersaturated  layers  is,  therefore,  highly  relevant. 
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It  is  also  important  to  mention  that  the  discrepancies  between  the  simula- 
tions and  the  mean-field  predictions  are  not  surprising  because  of  the  large  density 
fluctuations  that  exist  in  undersaturated  polymer  layers  in  good  solvents.  Neverthe- 
less, an  analysis  of  the  differences  between  SF  theory  and  the  simulations  (for  both 
undersaturated  and  saturated  surfaces),  and  efforts  to  improve  the  numerical  mean- 
field  predictions  are  important  since,  as  mentioned  earlier,  the  numerical  mean-field 
theory  predicts  the  only  easily  accessible  guidelines  for  the  structure  and  forces  in- 
duced by  short  physisorbed  polymers.  In  this  context,  it  is  interesting  to  make  a final 
comparison  between  SFi  and  SF2.  In  Figure  6.16  we  show  the  overall  density  profile 
obtained  for  F = 0.75  and  H/a  = 5.  The  agreement  observed  here  between  SFi  and 
simulations  is  fortuitous.  Although  the  two  density  profiles  (those  obtained  from  SFi 
and  simulations)  are  almost  identical,  the  forces  obtained  from  them  are  drastically 
different  from  each  other  (see  Figure  6.13).  This  illustrates  the  fact  that  agreements 
in  segment  density  profiles  or  surface  coverage  alone  are  not  sufficient  as  a measure 
of  the  acceptability  of  the  theoretical  predictions.  Finally,  SF2  predicts  a larger  con- 
centration at  the  midpoint  between  the  surfaces  (partly  due  to  the  prediction  of  a 
larger  number  of  bridges).  However,  it  underestimates  the  segment  concentration  at 
the  surface.  It  would  be  interesting  to  extend  the  numerical  mean-field  calculations 
to  also  incorporate  bond  correlations.  This  correction,  which  leads  to  the  so-called 
self-consistent  anisotropic  theory  [141,  142],  attempts  to  improve  over  the  random 
mixing  approximation  in  each  layer  (and,  hence,  improve  the  numerical  predictions). 
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Figure  6.16:  Segment  density  profile  between  two  adsorptive  walls.  The  coverage  in 
each  surface  is  F = 0.75  and  the  separation  between  walls  is  H/a  = 5. 


CHAPTER  7 

CONCLUDING  REMARKS 

In  the  previous  chapters  we  have  focused  on  the  calculation  of  polymer- 
induced  forces  using  specially  developed  computer  simulations.  Such  calculations 
pave  the  way  for  testing  the  scaling,  mean-field  and  density  functional  theories  of 
polymer  adsorption.  Moreover,  our  results  also  demonstrate  that  the  forces  of  inter- 
action provide  a more  discriminating  measure  of  the  acceptability  of  the  theoretical 
approximations  than  mere  comparisons  of  structural  properties.  Indeed,  both  are 
often  necessaxy.  Therefore,  the  approach  we  have  used  constitutes  an  essential  step 
for  improving  the  theoretical  predictions  for  routine  guidance  in  practice. 

A summary  of  the  key  contributions  of  this  work  is  already  presented  in 
Chapter  1,  Section  1.3.  In  this  chapter  we  present  some  concluding  remarks  and  final 
thoughts.  Our  goal  here  is  to  go  beyond  the  results  we  have  presented  in  the  previous 
chapters  and  to  (i)  emphasize  the  motivation  for  the  method  proposed  here,  and  (ii) 
elaborate  on  some  of  the  directions  we  recommend  for  the  future. 

7.1  Some  Comments  on  the  Rationale  for  CDM 
7.1.1  Forces  from  Lattice  Simulations 

It  is  important  to  emphasize  here  that  there  are  a large  number  of  problems 
where  lattice  simulations  axe  preferred  over  off-lattice  simulations.  In  addition  to 
the  fact  that  simple,  efficient  methods  for  sampling  polymer  conformations  are  more 
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easily  developed  for  lattice  models,  there  are  a large  number  of  lattice-based  theories 
{e.g.,  Flory-Huggins  theory,  Scheutjens-Fleer  theory)  that  require  examination  with 
a corresponding  lattice  simulation  [143].  While  in  off-lattice  simulations  the  forces 
can  be  obtained  easily  from  the  calculation  of  the  pressure  from  the  virial,  in  lattice 
models  the  forces  need  to  be  obtained  from  differences  in  the  Helmholtz  potential. 
The  contact-distribution  method  (CDM)  constitutes  a new,  efficient  technique  to  cal- 
culate such  forces  and  Helmholtz  potential  differences.  The  agreement  between  the 
results  obtained  from  lattice  and  off-lattice  approaches  (and,  hence,  the  correspon- 
dence between  CDM  and  results  obtained  from  the  virial)  has  been  established  in  the 
analysis  of  ‘simple’  single-chain  problems,  namely,  the  compression  of  a single  end- 
grafted  chain  with  a flat  bottom  object  [19],  and  the  confinement  of  a single  chain 
between  two  ‘athermal’  (neutral)  walls  [101].  In  the  case  of  polymers  confined  to  very 
small  spaces,  the  underlying  lattice  structure  may  have  some  effect  on  the  observed 
results^  . Simulations  with  different  lattice  models  {e.g.,  bond-fluctuation  model  and 
simple  cubic  lattice)  can  be  used  to  shed  light  on  that  issue. 


^ For  instance,  in  the  particular  case  of  compression  of  a single  end-grafted  chain 
with  a flat  bottom  object,  we  observe  a transition  region  {i.e.,  a noise  dominated  re- 
gion, where  the  chain  is  neither  confined,  nor  escaped)  that  extends,  approximately, 
over  three  lattice  units  (see  Chapter  4).  We  have  also  observed  that,  in  the  com- 
pression of  a polymer  brush,  the  transition  region  becomes  almost  imperceptible  [15]. 
These  results  were  obtained  with  the  bond-fluctuation  model  (BFM),  a model  in 
which  the  average  bond  length  is,  approximately,  2.5  lattice  units.  It  would  be  inter- 
esting to  examine  whether  such  a transition  can  be  observed  with  less  flexible  models 
{i.e.,  a simple  cubic  lattice). 
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7.1.2  Need  for  CDM 

The  need  for  efficiency  in  simulations  cannot  be  overemphasized  even  in  these 
days  of  fast  supercomputers.  For  instance,  the  time  needed  to  obtain  the  results 
presented  on  Chapters  5 and  6 is  of  the  order  of  months  (based  on  a single  pro- 
cessor). Extensions  to  more  complicated  systems,  containing  additional  interactions 
and  longer  chains  cannot  be  performed  without  efficient  simulation  procedures.  The 
contact-distribution  method  is  not  the  first  attempt  designed  to  obtain  forces  from 
lattice  Monte  Carlo  simulations.  As  discussed  in  Chapter  4,  Dickman  [24]  has  used  a 
thermodynamic  integration  approach  to  examine  the  interaction  between  two  polymer 
brushes.  The  CDM  is,  however,  more  efficient  than  Dickman’s  approach.  As  evident 
from  its  name,  CDM  uses  the  entire  distribution  of  contacts  between  the  polymer 
segments  and  the  compressing  object,  as  opposed  to  Dickman’s  method  where  only 
the  average  number  of  contacts  is  used.  Moreover,  Dickman’s  procedure  requires  an 
extrapolation  of  the  average  number  of  contacts  for  the  case  of  an  infinite  repulsive 
potential  applied  at  the  surface  layer.  In  CDM  such  an  extrapolation  is  not  needed. 
Overall,  the  approach  we  have  introduced  allows  for  a two  to  three-fold  reduction  in 
computer  time  at  the  expense  of  a slight  increase  in  the  allocated  memory. 

7.2  Immediate  Extensions 

We  plan  to  use  the  contact-distribution  method  for  scrutinizing  some  of  the 
basic  approximations  in  the  current  theories  of  polymer  layers.  In  this  context,  we 
have  initiated  a focussed  research  program  on  the  examination  and  development  of 
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improved  numerical  lattice  mean-field  theories.  In  the  short  term  our  efforts  will 
focus  on  the  following: 

7.2.1  Anisotropic  Mean-Field  Models 

We  have  already  shown  how  the  use  of  a second-order  Markov  chain  improves 
the  predictions  of  some  of  the  structural  features  of  polymer  layers  and  the  interac- 
tions induced  by  such  layers.  The  addition  of  the  anisotropic  potential  is  expected  to 
improve  the  calculation  of  the  excluded-volume  interaction.  In  the  standard  numeri- 
cal lattice  theory,  the  interactions  are  based  on  a random  mixing  approximation.  The 
anisotropic  model  improves  this  approximation  by  taking  into  account  correlations 
between  the  bonds  of  the  chain.  For  instance,  let  us  assume  that  in  one  particular 
conformation  of  a polymer  chain,  segment  s is  in  layer  i,  and  segment  s — 1 is  in 
layer  i — 1.  (The  bond  connecting  segments  s — 1 and  s is,  in  this  case,  a vector 
in  the  i direction.)  If  we  take  into  account  bond  correlations,  any  other  segment 
in  layer  z with  a similar  bond  vector  would  not  be  used  to  calculate  the  excluded- 
volume  interaction  of  segment  s.  This  anisotropic  model,  which  was  first  introduced 
by  Leermakers  and  Scheutjens  [141],  was  used  by  Fleer  and  coworkers  [44]  to  calculate 
the  surface  coverage  and  thickness  of  a polymer  layer.  Such  study,  however,  cannot 
offer  any  information  on  the  improvement  an  anisotropic  field  may  introduce.  It  is 
through  comparisons  with  the  type  of  simulations  we  have  shown  here  that  a proper 
assessment  of  this  extension  of  the  numerical  mean-field  theory  can  be  made. 
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7.2.2  Two-Dimensional  Numerical  Lattice  Mean-Field  Theory 

In  the  standard  numerical  lattice  theory,  a random  mixing  approximation  is 
used  in  each  layer  parallel  to  the  surface.  Therefore,  this  theory  cannot  be  used  to 
examine  interactions  between  polymer  layers  and  finite-sized  objects.  To  this  order, 
a two-dimensional  numerical  lattice  mean-field  theory  has  been  formulated  based  on 
the  use  of  a cylindrical  lattice  [144].  (In  this  case,  the  random  mixing  approximation 
is  applied  at  concentric  rings.)  This  2D  theory  may  have  important  implications  for 
developing  guidelines  for  interactions  between  polymer  layers  and  finite-sized  objects. 
A rigorous  examination  of  the  forces  induced  by  both  end-grafted  and  physisorbed 
chains  can  shed  light  on  the  potentials  and  limitations  of  this  theory. 

In  the  following  sections  we  present  some  additional  thoughts,  remarks  and 
recommendations.  These  comments  are  meant  to  draw  attention  to  the  fact  that  the 
investigation  of  the  fundamentals  of  the  theories  of  polymer  adsorption  has  just  begun. 
They  are  also  intended  to  inspire  some  additional  directions  for  future  research. 

7.3  An  Overview  of  Some  Outstanding  Issues 
The  scaling,  mean-field  and  density  functional  theories  have  focused  on  the 
study  of  a simple,  ideal  system:  adsorption  of  homopolymers  on  flat  homogeneous 
surfaces.  In  this  section  we  present  some  final  comments  that  complement  the  assess- 
ment of  the  theories  provided  in  Chapter  2. 

• Approximate  and  Exact  Solutions  to  the  Mean-Field  Approximation. 
The  ground  state  dominance  approximation  and  the  two-order-parameter  the- 
ory contain,  in  addition  to  the  mean-field  assumption,  a number  of  additional 
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approximations  which  are  needed  to  obtain  analytical  solutions  to  the  Edwards 
equation.  Although  the  technical  improvements  introduced  by  the  two-order- 
parameter  theory  are,  in  some  cases,  significant,  the  predictions  are  still  limited 
to  polymer  chains  of  infinite  molecular  weight.  In  contrast,  the  numerical  lat- 
tice theory  of  Scheutjens  and  Fleer  provides  a complete  solution  to  a discretized 
version  of  the  Edwards  equation.  This  theory  constitutes  the  only  approach 
that  has  the  potential  to  provide  guidelines  for  polymers  of  arbitrary  molecular 
weight.  Furthermore,  it  can  be  extended  to  study  the  behavior  of  systems  with 
greater  complexity,  such  as  the  adsorption  of  copolymers  or  the  adsorption  of 
branched  polymers.  Much  needs  to  be  done  in  this  regard. 

• Approximate  and  Exact  Solutions  to  Lattice  Models.  The  numerical 
lattice  mean-field  theory  and  the  lattice  Monte  Carlo  simulations  presented  in 
Chapters  5 and  6 are  based  on  the  same  model,  i.e.,  a polymer  chain  embedded 
in  a simple  cubic  lattice.  However,  while  the  results  of  the  simulations  are  based 
on  the  conformational  statistics  of  self-avoiding  walks,  the  theoretical  predic- 
tions are  based  on  the  conformations  of  a single  (first-  or  second-order)  Markov 
chain  in  a mean-field.  Although  the  mean-field  assumption  is  not  expected  to 
provide  correct  guidelines  in  systems  with  concentration  fluctuations  {e.g.,  in 
semi-dilute  solutions  or  starved  layers),  comparisons  between  simulations  and 
the  theoretical  calculations  allow  one  to  examine  the  extent  of  the  deviations 
in  the  theoretical  predictions  and  the  refinements  that  may  improve  the  accu- 
racy of  such  predictions.  Therefore,  a systematic  examination  of  the  lattice 


mean-field  theory  is  in  order. 
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• Does  Marginal  Solvent  Conditions  Exist  in  Lattice  Models?  Fleer  et  al. 
[7]  have  argued  that  a mean-field  approach  is  valid  under  the  so-called  marginal 
solvent  conditions.  In  good  solvents,  marginal  conditions  occur  somewhere  be- 
tween the  semidilute  and  the  concentrated  regimes  [7].  As  well  known,  the 
characteristic  length-scale  associated  with  the  density  fluctuations  in  semidi- 
lute polymer  solutions  is  known  as  the  correlation  length^  , This  length  is 
used  to  define  the  notion  of  ‘blobs’,  which  are  sections  of  a chain  where  inter- 
chain interactions  are  negligible  {i.e.,  the  segments  inside  the  blob  behave  as  an 
unperturbed  chain).  As  the  concentration  of  the  polymer  solution  increases,  the 
correlation  length  (or  the  blob  size)  decreases.  Fleer  et  al.  [7]  argue  that  there 
is  a critical  concentration  beyond  which  the  number  of  segments  inside  the  blob 
is  so  small  that  they  behave  as  an  ideal  chain.  It  is  above  this  concentration 
where  marginal  conditions  exist.  The  above  argument  is  indeed  correct,  but 
only  for  a special  class  of  models.  Recent  simulations  with  off-lattice  models 
have  shown  that  ‘short’  polymers  may  behave  as  either  ‘ideal’  or  as  ‘rigid-rods’, 
depending  on  the  magnitude  of  the  excluded- volume  interaction  [145].  For  in- 
stance, simulations  of  a pearl- necklace  model  [145]  {i.e.,  hard-spheres  connected 
by  rigid  bonds)  have  shown  that  the  alluded  ‘ideal’  behavior  of  short  chains  is 
observed  only  if  the  size  of  the  hard  sphere  is  small  compared  to  the  length  of 
the  bond.  These  results  imply  that  there  are  a large  number  of  models  (lattice 

^ The  correlation  length,  i^,  is  associated  with  the  average  distance  between  entan- 
glements in  a polymer  solution. 
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models  included)  where  the  existence  of  marginal  conditions  (in  good  solvents) 
is  highly  questionable. 

• Scaling  Theories  Provide  Asymptotic  Behavior  in  a Few,  Simple  Cases. 

Polymer  solutions  in  good  solvents  have  been  described  successfully  by  theo- 
ries based  on  simple  scaling  laws.  These  scaling  arguments  have  been  used 
to  describe  the  structure  and  interactions  of  saturated  polymer  layers  contain- 
ing chains  of  infinite  molecular  weight.  It  has  been  shown  recently  [62]  that 
the  scaling  prediction  of  the  power  law  decay  of  the  segment  density  profile 
is  in  agreement  with  results  obtained  from  simulations  of  long  self-avoiding 
walks  (containing  more  than  1000  segments).  The  asymptotic  prediction,  how- 
ever, differs  markedly  from  the  power  law  observed  from  simulations  of  shorter 
chains.  A closer  examination  of  this  is  warranted. 

• Scaling/Mean-Field  Hybrids:  What  Can  They  Do  for  Us?  A combi- 
nation of  scaling  arguments  and  the  free  energy  functional  approach  has  been 
used  to  explore  the  behavior  of  undersaturated  (starved)  layers  [9,  11,  12].  The 
formulation  of  such  scaling-based  functionals  is,  however,  somewhat  arbitrary. 
It  is  not  clear  how  well  such  an  approach  performs.  As  shown  in  Chapter  6,  the 
calculation  of  the  interaction  between  two  undersaturated  layers  can  provide  a 
rigorous  test  of  these  theories.  The  results  we  have  presented,  however,  were 
obtained  for  chains  containing  only  200  segments.  Therefore,  one  cannot  draw 
any  conclusions  regarding  the  accuracy  of  the  scaling/FEF  approach  until  one 
examines  the  behavior  of  longer  chains. 
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7.4  Nonequilibrium  Forces 

It  is  important  to  emphasize  that  all  the  interaction  forces  calculated  in  our 
study  are  based  on  Helmholtz  energies  of  interaction,  and  hence,  represent  the  forces 
one  would  obtain  in  a quasistatic  compression  (or  retraction)  of  the  surfaces.  If  the 
characteristic  time  of  compression  or  retraction  is  small  compared  to  the  relaxation 
time  of  the  polymers,  the  resulting  forces  may  differ  drastically  from  the  forces  we  have 
reported  here.  For  instance,  the  non-equilibrium  retraction  of  a single  adsorbed  chain 
may  not  give  rise  to  a constant  force  but  to  a rather  distinctive  profile  characterized  by 
both  entropic  and  enthalpic  contributions.  Similarly,  the  non-equilibrium  interaction 
between  two  starved  polymer  layers  may  be  repulsive  rather  than  attractive.  The 
possibility  of  using  computer  simulations  to  get  an  insight  on  non-equilibrium  forces 


needs  further  examination. 


APPENDIX  A 

LATTICE  MODELS  OF  POLYMER  SYSTEMS 
The  structure  of  polymer  systems  is  characterized  by  a wide  range  of  length 
scales  that  vary  from  m,  for  details  on  the  chemical  structure  of  the  chain,  to 
at  least  10~*-10~®  m,  for  details  on  the  meso-  or  microstructure  (e.p.,  coil  size,  layer 
thickness,  etc.).  Similarly,  the  dynamics  of  these  systems  are  characterized  by  time 
scales  that  range  from  s (for  bond  vibrations)  to  at  least  10“®  s for  the  relaxation 
of  a polymer  chain  or  other  collective  motions.  In  order  to  simulate  these  systems,  it 
is  necessary  to  omit  fast  motions  (and  structural  details)  on  very  small  length  scales 
and  focus  on  the  slow  motions  occurring  on  larger  scales.  Although  this  simplified 
(often  called  ‘coarse-grained’)  model  of  a polymer  chain  cannot  capture  differences 
introduced  by  the  chemical  architecture  of  the  chain,  it  does  capture  the  universal 
behavior  of  flexible  chains  and,  therefore,  it  is  very  useful  for  testing  the  validity  of 
theoretical  concepts  or  for  providing  guidelines  on  the  behavior  of  a particular  class 
of  polymers. 

A large  number  of  ‘coarse-grained’  models  have  been  implemented  in  the 
continuum  space.  Some  of  the  models  that  have  been  used  include  the  freely  jointed 
chain,  the  pearl  necklace  model,  and  the  bead-spring  model.  In  general  terms,  these 
models  use  interaction  potentials  that  depend  on  the  bond  length  (or  effective  bond 
length),  bond  angles,  torsional  angles  and  also  an  effective  potential  between  non- 
bonded  monomers.  A detailed  description  of  these  models  is  available  in  the  literatme 
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[146,  147].  Here  we  are  interested  in  models  that  contain  an  additional  level  of 
simplification,  namely,  the  discretization  of  the  space  by  a lattice  structure  [146]. 
These  models  are  very  popular  because  fast  algorithms  are  simple  to  implement.  In 
what  follows  we  describe  two  lattice  models  we  have  used  in  the  work  presented  in 
Chapters  4,  5 and  6. 


A.l  Simple-Cubic  Lattice  Model 

In  this  model  the  segments  of  the  polymer  are  located  at  lattice  sites,  and 
the  bonds  have  the  length  of  the  lattice  spacing,  a.  Each  lattice  site  can  be  occupied 
at  most  by  one  bead,  thereby  automatically  taking  into  account  the  excluded  volume 
interaction.  The  polymer  chains  in  this  model  are  therefore  considered  as  self-avoiding 
walks  (SAWs). 

The  conformations  of  self-avoiding  walks  can  be  sampled  using  a variety  of 
methods,  such  as  kink-jump,  crankshaft  motion  or  reptation  [146,  148].  These  meth- 
ods can  be  used  not  only  to  sample  the  conformations  of  the  polymer  chains,  but 
also  to  examine  the  dynamic  behavior  of  the  system  as  it  approaches  equilibrium^  . 
However,  these  methods  are  not  adequate  to  overcome  large  energetic  barriers  (i.e., 
the  system  may  get  trapped  in  a metastable  state),  and  hence,  they  are  not  adequate 


^ The  dynamic  information  obtained  from  such  an  approach  seems  to  capture  the 
‘real’  dynamics  provided  the  characteristic  time  associated  with  an  elemental  move 
in  the  simulation  is  much  smaller  than  the  characteristic  time  one  wants  to  examine. 
For  instance,  the  reptation  move  may  not  give  appropriate  ‘dynamic’  information 
for  the  relaxation  of  a polymer  chain  in  a solution,  but  it  seems  to  capture  the 
essential  features  of  phenomena  occuring  at  longer  time  scales,  such  as  the  kinetics  of 
adsorption  and  desorption  of  ‘short’  polymers  near  an  attractive  wall  [148,  149,  150]. 
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to  sample  the  conformation  of  long  polymers  confined  by  adsorbing  surfaces.  A more 
appropriate  simulation  method  for  such  a purpose  is  the  ‘configurational-bias’  method 
of  Siepmann  and  Frenkel  [134], 

A.  1.1  Configurational-Bias  Method 

The  configurational-bias  method  consists  of  ‘breaking’  an  existing  chain  into 
two  parts,  erasing  one  of  these  subchains  and  randomly  regrowing  it  using  an  algo- 
rithm based  on  the  Rosenbluth-Rosenbluth  method  [135].  The  detailed  procedure 
proceeds  as  follows; 

1.  Select  a chain  from  the  system,  randomly. 

2.  Select  the  end-segment  to  be  removed  (denoted  here  as  segment  A^),  randomly. 

3.  Select  the  segment  at  which  the  chain  is  broken  (say  segment  k)'^  , randomly. 

4.  Remove  segments  A:  -f- 1 to  A^. 

5.  Regrow  segments  starting  from  segment  1. 

6.  Accept  the  move  with  probability  Pacc  given  by 

Pace  ~ min[l,  Wold/^^trial]f  (■^•1) 

where  Wgid  and  Wtriai  are  the  Rosenbluth  weights  of  the  old  subchain  and  the 
trial  subchain,  respectively. 

^ To  improve  the  ratio  of  accepted  moves,  it  is  convenient  to  place  a limit  on  the 
size  of  the  subchain  to  be  regrown  (say  40-100  segments).  For  additional  details,  see 
Ref.  [62]. 
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In  the  original  version  of  the  configurational  bias  method  the  segments  are 
regrown  starting  from  segment  k.  However,  de  Joannis  [62]  has  shown  that  regrowing 
the  segments  on  the  opposite  side  results  in  a more  efficient  algorithm. 

As  mentioned  earlier,  the  simulations  on  a simple  cubic  lattice  are  very  efficient 
in  the  use  of  computer  time.  One,  however,  has  to  be  careful  with  the  influence  of  the 
underlying  lattice  structure  on  the  observed  results,  especially  in  the  case  of  polymers 
confined  in  small  spaces.  One  attempt  to  increase  the  flexibility  of  polymer  chains  in 
a lattice  is  the  so-called  bond-fluctuation  model.  The  basic  characteristics  of  such  a 
model  are  presented  in  the  next  section. 

A. 2 Bond-Fluctuation  Model 

The  bond-fluctuation  model  (BFM)  constitutes  an  intermediate  description 
between  a highly  flexible  continuum  model  and  the  traditional  lattice  model  [151]. 
Each  segment  (or  effective  monomer)  occupies  eight  sites  in  the  lattice,  as  illustrated 
in  Figure  A. 2.  The  vector  that  connects  two  segments  can  take  108  values  (which 
result  from  a combination  of  5 bond  lengths  and  87  bond  angles),  as  opposed  to  6 
values  in  the  simple  cubic  model.  The  fluctuation  in  the  bond  length  in  this  model  is 
physically  reasonable  since  it  corresponds  to  fluctuations  in  the  length  of  the  Kuhn 
or  statistical  segment^  of  a polymer  chain. 

^ A statistical  segment  is  defined  as  length  of  the  chain  whose  orientation  is  uncor- 
related to  its  neighboring  statistical  segments. 
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In  order  to  maintain  the  effective  performance  of  lattice  models,  the  set  of 
bond  lengths  in  BFM  is  restricted  to  those  bonds  where  the  excluded-volume  con- 
straint simultaneously  guarantees  that  bonds  cannot  cross  each  other.  This  set  of 
bond  lengths  depends  on  the  elementary  motion  of  the  segments  during  the  simu- 
lation. In  BFM,  the  elementary  motion  is  similar  to  the  kink-jump  motion  in  the 
simple  cubic  lattice,  i.e.,  the  position  of  one  segment  is  changed  from  r to  r -t-  Ar, 
where  Ar  = P{1,0, 0},  P{a,  ft,  c}  being  the  set  of  vectors  obtained  from  all  possible 
permutations  and  sign  combinations  of  ±a,  ±b  and  ±c. 

It  is  instructive  to  follow  the  derivation  of  the  set  of  bond  lengths  allowed 
in  the  model  [151].  As  evident  from  Figure  A. 2,  bond  intersection  is  possible  if  any 
coordinate  difference  between  segments  is  equal  or  larger  than  4.  After  elimination 
of  all  bond  vectors  with  Aa;  < 4 (or  Ay  < 4,  A^;  < 4),  the  remaining  set  can  be 
separated  in  two  groups  or  subsets  Bi  and  B2: 

Bi  = P(2, 0, 0)  U P(2, 1, 0)  U P(2, 1, 1)  U P(3, 0, 0) 

The  bonds  in  Bi  are  such  that  for  any  conformations  where  bonds  could  intersect 
the  segments  overlap. 

B2  = P(2, 2, 0)  U P(2, 2, 1)  U P(2, 2, 2)  U P(3, 1, 0)  U P(3, 1, 1)  U P(3, 2, 0)  U P(3, 2, 1) 
UP(3, 2, 2)  U P(3, 3, 0)  U P(3, 3, 1)  U P(3, 3, 2)  U P(3, 3, 3) 

The  bond  vectors  in  B2  are  such  that  bonds  might  cross  without  violating  the 


excluded-volume  constraint. 
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One  can  still  include  some  of  the  bond  vectors  from  the  subset  B2  provided 
one  eliminates  all  bond  vectors  that  are  not  connected  with  one  elementary  move 
to  any  vector  in  Bi  and,  at  the  same  time,  one  does  not  include  vectors  from  B2 
that  are  connected  by  one  elementary  move  to  any  other  bond  vector  of  B2.  These 
restrictions  lead  to  four  subsets  of  B2  that  can  be  incorporated  into  Bi: 

• P(2,2,0)UP(3,1,0)  or 

• P(2,2,0)UP(3,1,1)  or 

• P(2,2,1)UP(3,1,1)  or 

• P(2,2,1)UP(3,1,0). 

Prom  these  options,  the  subsets  P(2, 2, 1)  U P(3, 1, 0)  are  used  because  it  is 
the  option  that  contributes  with  more  bond  vectors  and  keeps  the  interval  of  bond 
lengths  as  small  as  possible.  In  this  manner,  the  set  of  bond  vectors  that  is  used  in 
the  3D  bond-fluctuation  model  is: 

B = P(2, 0, 0)  U P(2, 1, 0)  U P(2, 1, 1)  U P(2, 2, 1)  U P(3, 0, 0)  U P(3, 1, 0) 

Note  that  if  the  elemental  moves  of  the  segments  include  the  vectors  P{1, 1, 0} 
or  P{1, 1, 1},  z.e.,  if  the  segments  or  effective  monomers  could  jump  a larger  distance 
in  a single  step,  the  set  of  bond  vectors  in  the  model  must  be  restricted  to  Bi. 

The  bond-fluctuation  model  can  be  considered  as  a complement  to  the  simple- 
cubic  lattice  model.  While  the  latter  may  be  used  for  systems  with  large  energetic 
barriers  (provided  an  efficient  technique  such  as  configurational- bias  is  used),  BFM 
is  more  appropriate  to  examine  the  behavior  of  branched  or  tethered  polymers.  It 
is  also  more  appropriate  for  studying  the  dynamic  properties  of  polymer  systems, 
especially  when  confined  in  small  spaces. 
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In  the  next  section  we  present  some  comments  on  the  relation  that  ‘coarse- 
grained’ models  such  as  the  ones  presented  here  may  have  with  ‘real’  systems. 

A. 3 Mapping  Real  Systems  onto  ‘Coarse-Grained’  Models 

Although  the  basic  idea  of  using  ‘coarse-grained’  models  in  a simulation  is  to 
get  a better  understanding  of  the  universal  aspects  or  physics  underlying  a particular 
system,  it  is  a good  practice  to  get  a rough  idea  of  the  type  of  ‘real’  systems  that 
may  have  a behavior  similar  to  the  one  observed  in  simulations  [152,  153].  At  the 
same  time,  it  is  important  to  understand  how  a real  system  can  be  transformed  to  a 
coarse-grained  model. 

A simple  way  to  map  real  systems  onto  coarse-grained  models  is  to  establish 
the  correspondence  between  the  ‘flexibility’  of  a particular  polymer  and  the  flexibility 
of  the  lattice  model.  The  bond  length  in  a simple  cubic  lattice  is  1 lattice  unit,  while 
in  the  bond-fluctuation  model  it  is,  approximately,  2.5  lattice  units.  In  principle,  a 
bond  in  a lattice  model  represents  not  an  interatomic  bond  but  rather  one  statistical 
segment  in  a polymer  chain.  For  example,  for  the  case  of  poly  (styrene),  a statis- 
tical segment  has  a length  of,  roughly,  1 nm  and  consists  of  four  to  five  repeating 
(monomeric)  units.  One  can  therefore  say  that,  for  a poly(styrene)  system,  each  lat- 
tice unit  in  the  bond-fluctuation  model  represents  0.4  nm  and  that  each  unit  of  force 
( ksT/a  ),  therefore,  represents  about  10  pN. 

It  is  also  important  to  point  out  that  the  adsorption  energy  used  in  the  model 
is  actually  the  net  energy  of  adsorption  of  a segment,  which  depends  on  both  the 


characteristics  of  the  substrate  and  the  solvent. 
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One  should  emphasize,  however,  that  the  relations  presented  above  are  merely 
guidelines  used  to  obtain  a better  idea  of  the  magnitudes  of  the  forces  and  distances 
involved  in  the  calculations  and  to  facilitate  the  interpretation  and  analysis  of  both 
experiments  and  simulations. 


APPENDIX  B 

SCHEUTJENS-FLEER  THEORY 


The  lattice-based,  numerical,  mean-field  theory  of  Scheutjens  and  Fleer  pro- 
vides guidelines  for  the  structure  and  forces  induced  by  polymers  at  interfaces  or 
polymers  confined  in  small  spaces.  The  fundamental  ideeis  behind  this  theory  are 
introduced  in  Chapter  2,  Section  2.1.3.  In  this  appendix,  we  present  specific  de- 
tails on  the  implementation  of  the  numerical  mean-field  calculations  for  first-order 
and  second-order  Markov  chains.  (A  first-order  Markov  chain  admits  backfolding 
of  segments,  whereas  a second-order  Markov  chain  avoids  such  backfolding.)  The 
Scheutjens-Fleer  theory  considers  lattice  chains  confined  by  two  walls  and  filled  with 
solvent  molecules,  as  depicted  in  Figure  B.l.  We  are  particularly  interested  in  the 
so-called  restricted  equilibrium  case  (see  Chapter  2,  Section  2.3),  in  which  the  number 
of  polymer  chains  between  two  surfaces  is  fixed  and  independent  of  the  intersurface 
separation. 

The  lattice  system  depicted  in  Figure  B.l  has  an  area  of  L sites  (in  the  plane 
parallel  to  the  surface)  and  M layers  between  the  walls.  The  underlying  structure  of 
the  lattice  is  characterized  by  two  parameters,  Aq  and  Ai,  which  are  related  to  the 
fraction  of  nearest  neighbors  that  lie  in  the  same  and  adjacent  layers,  respectively.  For 
instance,  for  a simple  cubic  lattice,  Aq  = 4/6  and  Ai  = 1/6.  We  denote  the  number 
of  polymer  chains  in  our  system  as  n and  the  length  of  the  chains  as  r,  following  the 
convention  used  by  Scheutjens  and  coworkers  [7]. 
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Figure  B.l:  Lattice  system  confined  between  two  adsorbing  walls.  All  sites  are  filled 
with  either  polymer  segments  or  solvent  molecules. 

The  statistical  weight,  Gz,  of  a monomer  located  in  layer  z is  determined  by 
the  effective  potential  felt  by  a polymer  segment.  For  an  athermal  system  Gz  is  given 
as 

ln(?z  = + Xs2^M,i  + ln(l  — 0i),  (B.l) 

where  Xai  denotes  the  free-energy  of  interaction  between  a segment  and  the  surface 
next  to  layer  1 (bottom  surface  in  Figure  B.l)  and  Xs2  the  corresponding  free-energy 
of  interaction  between  a segment  and  the  surface  next  to  layer  M.  The  Kronecker 
delta  guarantees  that  the  surface  potential  is  only  felt  at  the  surface  layer  {i.e.,  at 
z = 1 and  z = M).  The  last  term  in  Eqn.  B.l  represents  the  excluded- volume 
interaction  within  a random-mixing  approximation. 

In  what  follows  we  first  focus,  in  Section  B.l,  on  the  calculations  for  a first- 
order  Markov  chain  (SFi),  as  presented  in  the  original  theory  of  Scheutjens  and  Fleer 
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[40].  In  Section  B.2  we  focus  on  the  calculations  for  a second-order  Markov  chain 
(SF2),  i.e.,  a chain  that  avoids  backfolding.  Preliminary  results  on  the  structure  of 
polymer  layers  obtained  from  these  theories  follow  in  Section  B.3. 

B.l  First-Order  Markov  Chain 

In  a first-order  Markov  chain,  the  position  of  segment  s -t- 1 depends  only  on 
the  position  of  segment  s.  The  statistical  weight  of  an  s-mer  having  the  last  segment 
in  layer  z,  denoted  here  as  G{z,s),  is  then  written  in  terms  of  a recurrence  relation: 


G{z,  s)  — [AqG(2:,  s — 1)  -|-  \]^G{^z  -f- 1,  s — 1)  -f-  XiG{z  — 1,  s — 1)]  Gz-  (®-2) 

The  use  of  this  relation  allows  one  to  obtain  the  statistical  weight  of  a chain  of  r 
segments  from  the  statistical  weight  of  a monomer  (defined  in  Eqn.  B.l): 
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Once  the  statistical  weights  are  known,  one  can  obtain  the  segment  density  (pi  in  layer 
i,  from  the  so-called  composition  rule: 

^ G(i,  s)G{i,  r - s -I- 1), 

‘ S=1 


(B.4) 
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Equation  B.4  represents  the  probability  of  joining  two  different  subchains,  both  having 
the  end-segment  in  layer  i.  (Note  that  the  prefactor  1/Gj  arises  due  to  the  double 
counting  of  the  end-segment  of  each  subchain).  The  other  term  in  the  prefactor  is  a 
normalization  constant  C,  which  is  given  by^ 


C = 


n 

IWY 


(B.5) 


where  G{r)  is  the  total  weight  of  an  r-segment  chain  in  the  lattice  system,  z.e., 

M 

G(r)  = 53G(i,r)  (B.6) 

i=l 

It  is  interesting  to  notice  that  if  we  drop  the  excluded-volume  term  from 
Eqn.  B.l  {i.e.,  if  we  consider  a system  of  ideal  chains),  one  can  obtain  the  density 
profile  directly  from  Eqns.  B.2  and  B.4.  For  the  case  of  non-ideal  chains,  a recursive 
or  iterative  procedure  is  required.  In  this  case  the  solution  of  Eqn.  B.4  is  used  to 
recalculate  the  weights  of  a single  segment,  Gz  (which  are  then  used  to  recalculate  the 
segment  density).  We  note,  however,  that  when  simple  iteration  is  used,  the  method 
does  not  converge  to  a self-consistent  solution.  Therefore,  the  segment  density,  <l>i,k+i, 
to  be  used  for  recalculating  the  weights  (in  iteration  A;  -f  1)  is  taken  as 


(B.7) 


^ The  normalization  constant  given  in  Eqn.  B.5  corresponds  to  the  restricted- 
equilibrium  case,  where  the  number  of  polymer  chains  between  the  surfaces  is  fixed 
[41].  Details  on  the  equilibrium  case  (where  the  polymer  system  is  in  equilibrium 
with  a bulk  solution)  are  given  in  Refs.  [7,  40]. 
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where  (j>i^new  represents  the  results  obtained  from  Eqn.  B.4  and  uj  is  an  adjustable 
parameter  0 < a;  < 1.  In  practice,  we  find  that  uj  in  the  range  between  0.001  and  0.1 
gives  appropriate  results.  The  longer  the  chain  length,  however,  the  more  difficult  it  is 
to  obtain  a self-consistent  solution.  For  very  long  chain  lengths  {N  > 1000),  the  iter- 
ation procedure  described  above  becomes  inefficient  and  other  numerical  procedures 
are  preferred  [40]. 

For  non-athermal  solvents,  the  procedure  to  calculate  the  segment  density  is 
identical  to  the  one  described  above.  In  this  case,  however,  the  weight  of  a single 
segment  is  defined  as 


InGz  — + Xs2^M,i  + ln(l  — 4>i)  + X {{4>i)  ~ {4>Tj)  (B-8) 


where  = 1 - 0^,  and  (0i)  = Ao0i  -1-  + XAi-x- 

Once  a self-consistent  solution  has  been  obtained,  one  can  calculate  the 
polymer-induced  forces  from  the  difference  in  Helmholtz  potential  between  two  sys- 
tems containing  M and  M — 1 layers,  i.e., 

/IfcM-i  = - 1)  - T(M),  (B.9) 


where  d is  a free  variable  that  indicates  the  separation  between  walls.  The  Helmholtz 
potential  of  the  system  is  given  by 
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It  is  important  to  emphasize  that  the  formulation  presented  in  this  section 
corresponds  to  a chain  that  admits  backfolding  (which  is  physically  forbidden).  In 
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order  to  examine  the  behavior  of  semi-flexible  chains  or  the  behavior  of  chains  where 
backfolding  is  prohibited,  one  should  consider  the  use  of  a second-order  Markov  chain. 
The  calculations  for  this  case  are  described  in  the  next  section. 


B.2  Second-Order  Markov  Chain 

In  a second-order  Markov  chain,  the  position  of  segment  s -I- 1 depends  not 
only  on  the  position  of  segment  s,  but  also  on  the  position  of  segment  s-1.  Therefore, 
the  statistical  weight  of  an  s-mer,  with  the  end-segment  in  position  z,  is  represented 
as  G{z,s,/3),  where  /?  indicates  the  position  of  segment  s — 1.  It  is  convenient  to 
define  /3  relative  to  the  position  of  the  end  segment.  For  instance,  for  a simple  cubic 
lattice,  /3  can  take  six  different  values,  as  illustrated  in  Figure  B.2.  The  weight  of  an 
s-mer  can  therefore  be  represented  as  a six-component  vector.  Each  component  of 
the  vector  represents  one  of  the  six  possible  positions  of  segment  s — 1,  i.e., 


G(z,»)=  G(z,s,l)  G{z,s,2)  G(z,s,Z)  G{z,s,4)  G{z,s,5)  G(z,s,6) 

(B.ll) 


Note  that  for  a monomer  (i.e.,  a single  segment),  the  different  values  of  /5  have  no 
physical  meaning.  One  can  still  define  the  weight  of  a single  segment  in  a vector  form 


as: 


G(^,l) 


G,/6  G,/6  GJ6  GJ%  GJ6  GJQ 


(B.12) 


An  important  consequence  of  the  vector  representation  used  for  the  statistical 
weights,  G(z,  s),  is  the  need  to  use  6x6  matrices  to  characterize  the  lattice  structure 
(in  contrast  to  the  scalars  Aq  and  Ai  used  in  Section  A.l).  Furthermore,  there  are 
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Figure  B.2:  Second-order  Markov  chain  formulation.  This  figure  illustrates  the  con- 
vention used  to  denote  the  position  of  segment  s — 1 relative  to  segment  s. 

two  sets  of  matrices  that  should  be  used  in  the  recurrence  relation.  One  set  for 
the  calculation  of  the  statistical  weight  of  a dimer,  which  does  not  suffer  from  the 
restrictions  due  to  backfolding,  and  the  other  for  the  calculation  of  the  weights  for 
an  s-mer  (with  s > 2).  The  statistical  weight  of  a dimer  is 


G{z,  2)  = [\oG{z,  1)  XiG(z  + 1,1)  + A_iG(^  - 1, 1)]  G(z,  1), 


(B.13) 
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where  the  A matrices  are  given  by 
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The  corresponding  weight  for  an  s-mer  (with  s > 2)  is  given  as 

G{z,  s)  = [AS^G(2:,  s — 1)  + A"^G(z  + 1,  s — 1)  + A"iG(^  — 1,  s — 1)]  G{z,  1), 

(B.14) 

with  the  set  of  matrices  A"*’ (the  nb  superscript  stands  for  no-backfolding)  given  by 

1/5  1/5  0 1/5  1/5  1/5 

0 0 0 0 0 0 

_ 0 1/5  1/5  1/5  1/5  1/5 

^0  - , 

0 0 0 0 0 0 

1/5  1/5  1/5  1/5  1/5  0 

1/5  1/5  1/5  1/5  0 1/5 

0 0 0 0 0 0 

1/5  1/5  1/5  0 1/5  1/5 

^ 0 0 0 0 0 0 
0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 
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Af 


0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

1/5  0 1/5  1/5  1/5  1/5 


0 0 0 0 0 0 


0 0 0 0 0 0 


Once  the  weights  of  the  r-segment  chains  have  been  calculated,  the  segment 
densities,  (j>i,  are  obtained  from  the  composition  rule,  which  is  rewritten  for  a second- 
order  Markov  chain  as 

2 G{z,l,k)  G{z,r,k) 


<i>i  = C 


A;=l 


k=l 

gT~ 


(B.15) 


E 


^ G{z,  s,  k)  Y!,  G{z,  r - s + l,k)  - Y G{z,  s,  k)G{z,  r - s + l,k) 

k=l fe=l fc=l 


Note  that  the  first  term  in  Eqn.  B.15  is  identical  to  the  terms  used  in  the  first- 
order  Markov  chain  (see  Eqn.  B.4).  The  second  term,  in  contrast,  excludes  the 
conformations  of  the  subchains  that  cause  backfolding.  The  normalization  constant 
is  again  obtained  from  Eqn.  B.5,  with  the  total  weight  of  the  r-segment  chain,  G{r), 
given  by 

M 6 

G{r)  = '^'^G{i,r,k).  (B.16) 

t=i  jfe=i 

The  iteration  procedure  used  to  obtain  a self-consistent  solution  is  the  same 
as  the  one  described  in  the  previous  section.  Once  this  self-consistent  solution  is  ob- 
tained, all  the  structural  details,  such  as  the  contribution  to  the  segment  density  due 


184 


to  loops,  tails,  trains  and  bridges,  can  be  calculated.  The  details  of  these  calculations 
are  presented  in  Ref.  [7]. 

In  what  follows  we  provide  a preliminary  comparison  of  the  structure  of  a 
polymer  layer  predicted  by  the  first-order  (SFi)  and  the  second-order  (SF2)  formula- 
tions. 


B.3  Structure  of  a Polymer  Layer:  Comparison  between  SFi  and  SF2 
As  mentioned  earlier,  the  original  theory  of  Scheutjens  and  Fleer  is  based 
on  the  conformations  of  a first-order  Markov  chain  subject  to  a mean-field  potential. 
This  formulation  contains  two  critical  assumptions,  namely,  (z)  random  mixing  ap- 
proximation in  each  layer,  and  (ii)  the  possibility  of  chain  backfolding.  Here  we  relax 
the  second  assumption  by  using  a second-order  Markov  chain  for  the  step-weighted 
walk.  This  extension  of  the  standard  SF  theory  was  introduced  by  Scheutjens  and 
Leermakers  [140]  and  later  used  by  van  der  Linden  et  al.  [44]  to  examine  the  ad- 
sorption of  semi-flexible  copolymers.  It  is  interesting  to  mention  that  for  a long  time 
it  was  believed  that  the  use  of  a second-order  Markov  chain  would  only  modify  the 
length  of  the  statistical  segment  of  the  chains  (as  it  does  for  polymers  in  solution) 
[7].  However,  it  is  now  recognized  that  the  structure  of  a polymer  layer  is  indeed 
sensitive  to  such  type  of  assumptions  [44].  A detailed  analysis  of  the  layer  structure 
calculated  for  a second-order  Markov  chain  has  not  been  presented  yet.  In  what  fol- 
lows we  present  some  preliminary  results  on  the  structure  of  a layer  obtained  when 
the  separation  between  the  walls  is  large  and  a ‘bulk’  region  develops  in  the  middle 
of  the  gap. 
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Figure  B.3:  Surface  coverage,  F,  as  a function  of  the  bulk  concentration,  4>i,. 

We  first  present,  in  Figure  B.3,  the  surface  coverage,  F,  obtained  for  different 
values  of  the  bulk  concentration,  (f)b-  We  define  the  surface  coverage  as  the  number  of 
segments  that  belong  to  adsorbed  chains  divided  by  the  area  of  the  surface.  We  also 
present,  in  Figure  B.4,  the  segment  density  profile  as  well  as  the  contribution  of  loops 
and  tails  to  the  density  profile.  It  is  interesting  to  note  that,  although  the  segment 
density  on  the  first  layer  is  larger  for  SFi,  both  the  surface  coverage  and  the  thickness 
of  the  layer  are  larger  for  the  SF2  calculation.  The  ability  of  the  segments  to  back- 
fold  in  the  SFi  formulation  produces  an  overestimation  of  the  fraction  of  trains.  In 
particular,  the  SFi  calculation  allows  for  the  existence  of  trains  of  just  one  segment, 
contrary  to  SF2,  where  such  trains  cannot  exist.  This  observation  is  further  illus- 
trated in  Figure  B.5,  where  we  show  the  calculation  of  the  total  fraction  of  segments 
belonging  to  loops,  trains  and  tails  for  different  values  of  the  bulk  concentration. 
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Figure  B.4:  Segment  density  profile  and  contribution  of  loops  and  tails  obtained  from 
SFi  and  SF2  for  a polymer  layer  in  equilibrium  with  a bulk  concentration  06  ~ 10~^. 
The  continuous  and  broken  lines  correspond  to  the  predictions  from  SFi.  The  rms 
thickness  of  the  layer  (in  lattice  units)  is  2.2  for  SFi  and  2.7  for  SF2. 


Figure  B.5:  Fraction  of  segments  that  belong  to  trains,  tails  and  loops  as  a function 
of  the  bulk  concentration,  06-  For  instance,  for  06  ~ 10“^  (see  Figure  B4),  SFi 
predicts  that  69%  of  the  segments  belong  to  trains,  23%  to  loops  and  8%  to  tails. 
The  corresponding  values  obtained  from  SF2  are  61%,  30%  and  9%,  respectively. 
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The  results  shown  above  are  just  meant  to  illustrate  some  of  the  structural 
differences  one  can  obtain  when  the  conformations  where  the  chain  backfolds  are 
excluded  from  the  ensemble  average.  Additional  calculations  for  a wide  range  of  bulk 
concentrations,  chain  length  and  solvency  conditions  should  be  obtained  to  properly 
asses  the  differences  between  the  two  formulations  (SFi  and  SF2).  Furthermore,  it 
would  be  interesting  to  improve  the  assumption  of  random-mixing  on  each  layer  by 
accounting  for  bond  correlations  (ie.,  using  the  so-called  self-consistent  anisotropic 
theory  [141,  142]).  We  expect  the  segment  density  in  the  first  layer  to  increase 
(particularly  for  the  SF2  calculation)  when  bond  correlations  are  taken  into  account. 
This  in  turn  will  increase  the  prediction  of  the  surface  coverage  and  of  the  thickness  of 
the  layer.  Comparison  between  these  calculations  and  computer  simulations  can  then 
be  used  to  identify  the  most  useful  refinements  needed  for  improving  the  accuracy  of 


the  numerical  mean-field  calculations. 
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